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Abstract 


Experiments  were  conducted  on  lifted,  turbulent  jet  diffusion  flames.  An  au¬ 
tomated  technique  using  a  linear  photodiode  array  was  implemented  to  measure 
the  temporal  history  of  the  liftoff  height  h.  The  measurements  enabled  accurate 
determination  of  the  mean  liftoff  height  h  under  a  wide  range  of  flow  conditions, 
including  severed  fuels,  nozzle  diameters,  and  exit  velocities  Us.  The  results  showed 
an  approximately  linear  relationship  between  h  and  Us,  with  a  slight  dependence 
on  Reynolds  number.  A  strain-rate  model  for  liftoff,  based  on  far-field  scaling  of 
turbulent  jets,  provides  an  explanation  for  the  linear  dependence  of  h  on  Us.  Mea¬ 
surements  were  also  made  in  which  the  nozzle  fluid  contained  varying  amounts  of 
air,  where  it  was  found  that  the  slope  of  the  h  vs.  Us  line  increases  faster  than 
predicted  by  far-field  scaling  of  turbulent  jets.  The  discrepancy  is  attributed  to 
near-field  effects. 

The  amplitudes  of  the  fluctuations  in  h  were  found  to  be  of  the  order  of  the 
local  large  scale  of  the  jet.  There  is  a  slight  increase  in  normalized  fluctuation  level 
h' /h  with  h,  and  there  is  some  variation  of  h' /h  with  fuel  type.  The  time  scales 
of  the  fluctuations  of  h  were  found  to  be  considerably  longer  than  the  local  large- 
scale  time  of  the  turbulence  Tg.  By  using  fuels  of  different  chemical  times  to  vary 
T(,  the  measured  correlation  time  rj/2  normalized  by  was  found  to  collapse  with 
Richardson  number  £h-  Experiments  in  which  the  nozzles  were  oriented  horizontally 
showed  no  change  in  rj/2,  however.  Additional  experiments  were  conducted  to 
investigate  alternative  explanations  for  the  variation  of  ti/2/t«  with  £&.  These 
experiments  included  measuring  the  flame  length  L  simultaneously  with  h ,  and 
measuring  the  visible  radiation  I  simultaneously  with  h.  L(t)  was  found  to  be 
nearly  uncorrelated  with  h(t),  dismissing  the  possibility  that  a  feedback  mechanism 
from  L  to  h  controls  the  fluctuations  of  h.  Although  I(t )  is  highly  correlated  with 
h(t)  for  the  most  sooting  fuel,  acetylene,  it  is  not  deemed  responsible  for  the  longer 
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correlation  times  of  that  fuel.  This  was  deduced  from  experiments  using  mixtures 
of  hydrogen  with  other  fuels,  which  produce  very  little  radiation,  but  which  have 
values  of  ti/2/t6  comparable  to  those  of  acetylene  flames. 

Another  experiment  was  conducted  in  which  two-dimensional  images  of  fuel 
concentration  (CH4)  and  reaction  zones  (indicated  by  CH)  were  obtained.  The 
images  showed  a  wide  variety  of  structure  types,  indicating  that  there  is  no  universal 
description  of  the  flow  field  at  the  flame  base.  The  flame  stabilization  position 
showed  large  fluctuations  in  both  the  axial  and  radial  directions.  The  shot  to  shot 
variation  in  methane  number  density  at  the  flame  base  was  also  large. 
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CHAPTER  1 


Introduction 


Many  combustion  devices  employ  turbulent  flow  to  achieve  efficient  mixing  of 
fuel  and  oxidizer.  In  particular,  turbulent  jets  are  used  extensively.  To  reach  a  high 
combustion  efficiency,  conditions  are  sometimes  chosen  such  that  the  flame  stability 
is  a  concern.  There  are  many  factors  which  determine  the  stability  of  the  flame, 
including  the  rate  of  chemical  reactions,  strain  rates  in  the  flow,  and  the  mixing  of 
products  (which  serve  as  the  ignition  source)  with  the  reactants. 

The  lifted,  turbulent  jet  flame  is  a  useful  flow  to  investigate  these  complex 
issues.  Consisting  of  a  jet  of  fuel  issuing  from  a  round  nozzle  into  an  oxidizer  (air), 
it  has  a  relatively  simple  geometry.  Under  certain  flow  conditions,  the  flame  is 
detached  from  the  nozzle  (lifted).  With  such  a  flame,  there  is  a  region  from  the 
nozzle  to  the  flame  base  in  which  the  fuel  mixes  with  the  oxidizer,  but  no  burning 
ccurs.  Further  downstream,  where  the  velocities  and  strain  rates  are  lower,  the 
flame  is  stabilized.  A  shadowgraph  image  of  a  lifted,  turbulent  propane  flame  is 
shown  in  Fig.  1.1.  The  mixing  of  the  turbulent  jet  below  the  flame  base  is  visible, 
as  is  the  abrupt  transition  to  the  flame  base,  where  there  is  a  substantial  change 
in  the  structure  of  the  turbulence,  as  evidenced  by  the  initial  change  in  spreading 
angle  of  the  flow.  Because  the  flow  below  the  flame  base  is  largely  unaffected  by 
the  combustion  occurring  above,  it  is  the  much-studied  turbulent  jet. 

The  goal  of  this  research  is  to  gain  an  understanding  of  the  stabilization  mecha¬ 
nism  of  lifted  flames.  In  particular,  the  temporal  behavior  is  of  interest.  It  is  hoped 
that  knowledge  gained  from  studying  lifted  flames  can  be  applied  to  other  problems 
in  flame  stability. 


Fk;.  1.1  Shadowgraph  image  of  a  liftod  turbulent  propane  flame.  The  nozzle  di¬ 
ameter  is  3.S  111m.  the  exit  velocity  is  GG  111/s.  and  the  liftoff'  height  is  ~ 
lo()  mm. 

copy  A". ; r.r.z  :r.: '  >?4  m  •  • 
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1.1  Background 

Lifted  diffusion  flames  have  been  studied  by  many  investigators  in  the  past  ( e.g., 
Vanquickenborne  and  Vein  Tiggelen  1966,  Kalghatgi  1984,  Peters  and  Williams  1983, 
Eickhoff  et  al.  1985,  Lockwood  and  Moneib  1982).  Some  studies  have  focused  on 
flames  stabilized  in  the  near-field  region  of  the  jet  in  which  the  shear  layer  begins 
to  roll  up  (e.g.,  Sava s  and  Gollahalli  1986a,  Gollahalli  et  al.  1987).  In  particular, 
the  processes  of  detachment  from  the  nozzle  and  reattachment  to  the  nozzle  have 
received  considerable  attention.  Other  studies  have  examined  the  behavior  of  flames 
stabilized  further  downstream  (up  to  ~  50  nozzle  diameters,  depending  on  the  fuel). 
In  addition,  Savas  and  Gollahalli  (1986b)  and  Chung  and  Lee  (1991)  have  investi¬ 
gated  lifted  laminar  flames,  in  part  to  provide  a  starting  point  for  understanding 
the  turbulent  flames. 

Experimental  results  of  Kalghatgi  (1984)  have  shown  that  the  mean  (time- 
averaged)  liftoff  distance  h  of  turbulent  flames  increases  approximately  linearly  with 
nozzle  exit  velocity  U0,  independent  of  the  nozzle  diameter  d,  once  h  is  sufficiently 
far  downstream  (approximately  20 d).  The  slope  of  the  h  vs.  Ua  line  depends  on 
the  fuel  used.  At  a  maximum  velocity  Ub,  proportional  to  d,  the  flame  blows  out. 
Because  h  ~  I70,  blowout  occurs  when  h  reaches  a  maximum  distance  proportional 
to  d,  for  a  given  fuel. 

Several  models  have  been  proposed  to  explain  the  liftoff  behavior.  These  mod¬ 
els  can  generally  be  divided  into  two  categories,  those  which  propose  that  the  sta¬ 
bilization  distance  is  governed  by  a  balance  between  the  propagation  rate  of  the 
flame  and  the  velocity  of  the  approaching  flow,  and  those  which  propose  that  the 
liftoff  location  is  the  most  upstream  position  in  which  the  turbulent  mixing  does 
not  extinguish  the  combustion.  The  former  type  of  model  is  described  by  Kalghatgi 
(1984),  in  which  an  expression  for  the  turbulent  flame  speed  in  terms  of  the  lam¬ 
inar  flame  speed  and  the  turbulence  properties  is  used.  Models  of  the  later  type 
include  the  mixing-rate  model  of  Broadwell  et  al.  (1985),  the  scalar-dissipation-rate 
model  of  Peters  and  Williams  (1983),  and  our  own  strain-rate  model  (Miake-Lye 
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and  Hammer  1989).  Additional  classifications  of  the  models  in  terms  of  the  impor¬ 
tance  given  to  the  large-scale  structure  of  the  turbulence  can  also  be  made.  Another 
model  which  collapses  the  data  well  but  has  no  theoretical  basis  has  been  proposed 
by  Pitts  (1989). 

Although  most  of  the  experiments  on  flames  stabilized  in  the  far  field  have 
focused  on  the  mean  behavior,  there  have  been  some  measurements  of  fluctuating 
properties.  Lockwood  and  Moneib  (1982)  measured  fluctuating  temperatures  with 
thermocouples,  using  a  compensation  technique  to  extend  the  frequency  range  of 
the  measurements.  Chen  et  al.  (1989)  measured  temperatures  along  a  line  using 
thin  film  pyrometry,  allowing  them  to  track  the  radial  position  of  the  reaction  zone. 
Using  a  video  camera,  Birch  and  Hargrave  (1989)  measured  fluctuating  liftoff  heights 
of  choked,  underexpanded  natural  gas  flames. 

There  are  many  unresolved  issues  concerning  the  mechanism  by  which  lifted 
flames  are  stabilized,  and  the  present  experiments  were  undertaken  to  help  elucidate 
the  stabilization  process. 


1.2  Current  Experiments 

To  improve  upon  previous  measurements  of  the  liftoff  height,  a  new  measure¬ 
ment  system  was  designed  in  which  the  flame  light  was  focused  directly  onto  a  linear 
photodiode  array  aligned  with  the  axis  of  the  flame.  The  liftoff  height  could  then  be 
found  by  detecting  the  visible  radiation  at  the  flame  base.  This  enabled  temporal 
information  to  be  obtained,  in  addition  to  the  time-averaged  liftoff  height. 

Measurements  of  the  time  history  of  the  liftoff  height  were  made  for  moderate 
durations,  for  the  purpose  of  determining  the  mean  value.  The  fuels  used  in  these 
experiments  were  natural  gets  (abbreviated  as  n.g.;  93.5%  CH4,  3.1%  C2H6,  1.9% 
N2,  0.6%  C3H8,  0.6%  C02,  0.2%  C4H]0,  0.1%  C5H12),  ethane  (C2H«),  ethylene 
(C2H4),  and  acetylene  (C2H2).  For  each  of  the  four  fuels,  two  or  three  different 
nozzles  of  varying  diameter  were  used.  Liftoff  heights  ranging  from  the  maximum 
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attainable  for  each  nozzle  down  to  approximately  half  of  that  height  were  measured. 
Limiting  the  conditions  to  this  range  ensured  that  the  flames  studied  were  stabilized 
far  enough  downstream  of  the  nozzle  exit  that  self  similar  behavior  of  the  jet  below 
the  flame  base  was  achieved,  at  least  to  a  reasonable  approximation.  In  addition  to 
jets  of  pure  fuel  exiting  the  nozzle,  measurements  were  made  in  which  the  fuel  was 
premixed  with  air,  but  still  fuel  rich  and  thus  requiring  ambient  air  for  complete 
combustion. 

More  extensive  measurements  were  also  made  with  all  four  fuels  in  which, 
for  each  fuel,  two  different  nozzle  diameters  were  used  and  four  different  liftoff 
heights  studied  for  each  nozzle.  Using  two  nozzles  for  each  fuel  allowed  the  ef¬ 
fects  of  Reynolds  number  to  be  investigated,  while  using  different  fuels  provided  a 
mechanism  for  varying  the  local  fluid-mechanical  time  scale  as  well  as  a  buoyancy 
parameter  relevant  to  the  dynamics  of  the  flame  base.  From  the  longer  data  sets 
statistical  analyses  of  the  temporal  fluctuations  were  obtained.  In  addition  to  the 
above  listed  fuels,  experiments  were  done  with  fuel  mixtures.  A  mixture  of  50% 
C2H2  /  50%  C2H6  was  used,  primarily  to  compare  with  C2H4  flames.  Mixtures 
of  hydrogen  (H2)  with  each  of  the  fuels  C2H6,  C2H4,  and  C2H2  were  also  used  to 
achieve  greater  control  over  the  chemical  rate  of  the  fuel. 

As  an  additional  test  of  the  importance  of  buoyancy  to  the  liftoff  height  dy¬ 
namics,  experiments  were  conducted  with  the  nozzles  oriented  horizontally.  The 
fuels  used  in  these  experiments  were  natural  gas,  C2H4,  and  C2H2. 

To  investigate  whether  the  behavior  at  the  flame  tip  affects  the  liftoff  height, 
additional  experiments  were  conducted  in  which  the  liftoff  height  and  the  flame 
length  were  measured  simultaneously.  The  fuels  investigated  in  these  tests  were 
natural  gas  and  C2H6. 

Similarly,  simultaneous  measurements  were  made  of  the  liftoff  height  and  the 
flame  radiation  at  a  location  downstream  of  the  liftoff  height.  This  was  done  to  see 
whether  the  flame  radiation  fluctuations  are  correlated  with,  and  might  affect,  the 
liftoff  height  of  the  more  radiative  fuels  (C2H2  in  particular). 
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To  better  understand  the  flame  structure  and  the  stabilization  process,  experi¬ 
ments  were  conducted  in  which  planar,  two-dimensional  images  of  both  the  fuel  and 
reaction  zones  of  a  methane  (CH4)  flame  were  obtained  by  simultaneous  measure¬ 
ments  of  Raman  scattering  from  CH4  and  fluorescence  from  CH,  which  marks  the 
reaction  zones.  These  experiments  were  done  in  collaboration  with  Sandia  National 
Labs  in  Livermore,  CA.  In  the  flow  condition  reported  here,  the  fuel  was  premixed 
with  air,  and  the  flame  was  lifted  to  an  intermediate  liftoff  height,  ~  0.6  times  the 
maximum  height  before  blowout. 
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CHAPTER  2 

Experimental  Facilities 


2.1  Liftoff  Height  Measurements 

2.1.1  Nozzles  and  Flow  Control 

The  nozzles  used  in  the  liftoff  height  experiments  performed  at  Caltech  were 
designed  to  produce  fully  developed  turbulent  pipe  flow  at  the  exit.  All  nozzles  were 
made  from  brass  tubes,  except  the  smallest  nozzle  which  was  made  from  stainless 
steel.  The  nozzles  had  length  to  diameter  ratios  of  ~  60.  The  outside  of  each 
nozzle  was  tapered  at  a  15  degree  half  angle  to  reduce  disturbances  to  the  entrained 
flow  near  the  exit.  If  the  jet  is  modeled  as  a  semi-infinite  line  sink,  the  streamline 
approaching  the  nozzle  exit  is  a  parabola.  Although  a  parabola  is  then  the  ideal 
choice  for  the  external  shape  of  the  nozzle,  a  straight  taper  was  chosen  because  it 
is  simple  to  construct,  and  it  does  not  directly  obstruct  the  natural  streamline.  To 
facilitate  connection  to  a  plastic  tube  supp.ying  the  fuel,  the  nozzles  were  soldered  to 
successively  larger  brass  tubes  until  reaching  an  outer  diameter  of  9.5  mm  (3/8  in.). 
The  entrance  to  the  final  (smallest)  diameter  tube  of  each  nozzle  was  sharp  edged 
to  produce  turbulent  flow  right  from  the  beginning  of  the  tube.  Figure  2.1  shows 
the  design  of  the  2.16  mm  nozzle.  Several  different  nozzle  diameters  were  used,  to 
have  flexibility  in  adjusting  the  Reynolds  number  and  liftoff  height  for  the  various 
fuels.  The  nozzle  diameters  used  in  these  experiments  were  1.02,  1.55,  2.16,  3.12, 
3.88,  4.61,  4.92,  6.25,  and  7.73  mm.  For  all  of  the  nozzles,  the  exit  was  located 
53  cm  above  the  floor. 


8 


Fig.  2.1  Design  of  the  2.16  mm  nozzle.  Brass  tubes  with  outer  diameters  of  9.5  mm, 
6.4  mm,  and  3.2  mm  axe  soldered  together  as  shown,  reaching  a  final  inner 
diameter  of  2.16  mm. 

Combustion  gases  were  carried  out  of  the  lab  through  a  fume  hood  using  a 
blower.  The  hood  was  1.5  m  square  and  located  2.4  m  above  the  laboratory  floor. 
The  blower  had  a  flow  rate  of  1.2  m3/s.  Screens  were  placed  around  the  flame  to 
minimize  the  disturbances  from  room  drafts.  A  double  layer  of  screens  was  used 
around  the  perimeter  of  the  hood  and  an  additional  cylindrical  screen  of  diameter 
68  cm  and  height  120  cm  was  placed  around  the  flame.  In  the  studies  of  horizontal 
flames,  the  nozzles  protruded  through  the  double  layer  of  screens  around  the  hood, 
and  the  cylindrical  screen  was  suspended  horizontally.  All  screens  had  an  open  area 
of  70%,  and  were  painted  black  to  reduce  reflections  of  light. 

Flow  rates  of  the  fuels  were  controlled  with  valves  and  measured  with  laminar 
flow  elements  (LFE’s),  Meriam  models  50MW20-1  and  50MJ10-9.  The  pressure  drop 
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Ap  across  the  LFE  was  measured  with  a  Baxocel  pressure  transducer  type  590-D, 
with  a  full  range  of  10  torr.  The  upstream  pressure  p(m  (the  subscript  fm  stands  for 
conditions  at  the  LFE)  was  measured  with  an  Ashcroft  test  gauge,  with  a  range  of 
0  to  15  psi.  The  gas  temperature  upstream  of  the  flow  element  was  measured  with 
a  type  K  thermocouple  to  determine  the  viscosity  /ifm  and  density  ptm.  Values  of 
viscosity,  required  for  converting  pressure  drop  into  flow  rate,  were  taken  from  Reid 
et  al.  (1977). 

The  volumetric  flow  rate  at  the  LFE  was  obtained  from  the  solution  of  the 
equation 

Ap  =  Ci  p{m  Qim  +  C2  Plm  Qfm  i  (2-1) 

where  Q  is  volumetric  flow  rate,  and  Ci  and  C2  are  constants  from  a  calibration  of 
the  LFE.  The  quadratic  term  in  Eq.  2.1  was  always  small. 

The  mean  velocity  Ua  at  the  nozzle  exit  was  then  found  by  assuming  that  the 
gas  reached  atmospheric  pressure  and  temperature,  i.e., 

jj  _  Qfm  (Pfm/Poo)  (Too/T{m ) 

°  ”  (tt/4)  cP 

where  the  subscript  00  stands  for  ambient  conditions  present  at  the  nozzle  exit. 

The  natural  gas  used  in  the  experiments  was  gas  was  taken  from  a  commercial 
supply  line,  with  its  composition  determined  from  the  gas  company  sample  analysis. 
The  purities  of  the  bottled  fuels  were:  C2H6  99%,  C2H4  99.5%,  C2H2  99.6%,  and 
H2  99.95%.  Warm  water  was  used  to  keep  the  regulators  from  icing,  and  a  water 
bath  was  used  upstream  of  the  LFE  to  bring  the  gases  close  to  room  temperature. 
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2.1.2  Instrumentation  and  Measurement  Technique 

For  measuring  the  liftoff  height,  a  Reticon  linear  photodiode  array  was  used. 
The  array  was  a  type  T  with  512  pixels  spaced  at  50  pm;  each  pixel  was  2.5  mm 
wide.  The  visible  radiative  light  emission  was  imaged  directly  onto  the  array  using 
various  camera  lenses.  The  lenses  used  were  a  50  mm  f/1.4  Nikon  lens,  a  90  mm 
f/2.5  Pentax  lens,  and  a  75 — 300  mm  f/5.6  Pentax  zoom  lens.  The  particular  lens 
used  in  a  given  experiment  was  chosen  based  on  the  size  of  the  flame  and  its  light 
intensity. 

Data  from  the  linear  array  were  collected  using  an  Everex  1800  personal  com¬ 
puter  (having  an  Intel  286  microprocessor)  with  an  R.C.  Electronics  data  acquisition 
board  installed.  The  array  was  driven  by  a  clock  from  an  evaluation  board  supplied 
by  Reticon,  which  was  also  used  to  synchronize  data  acquisition. 


FlG.  2.2  Sample  scan  of  the  linear  array  imaging  a  C2H6  flame  to  illustrate  the 
method  for  determining  the  liftoff  height. 

FYom  each  seem  of  the  array,  the  liftoff  height  was  found  by  looking  for  the  first 
pixel  with  a  voltage  (light  intensity)  above  a  threshold  value.  This  is  illustrated 
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in  Fig.  2.2,  which  shows  one  scan  of  the  linear  array  imaging  a  C2H6  flame.  The 
background  level  of  the  array,  the  threshold  set  some  level  above  the  background, 
and  the  actual  array  scan  are  all  shown.  There  is  a  rather  sudden  increase  in  light 
intensity  at  the  flame  base. 

For  fuels  producing  large  amounts  of  soot  (C2H4  and  C2H2),  the  radiation  from 
which  results  in  background  light  in  the  images,  a  slight  modification  was  made  to 
the  algorithm  for  finding  the  liftoff  height.  Because  this  background  noise  was  found 
to  be  relatively  uniform  along  the  array  extent,  it  could  be  subtracted.  In  practice, 
before  comparing  a  pixel’s  intensity  with  the  threshold  for  determining  liftoff  height, 
the  intensity  32  pixels  below  the  current  pixel  was  subtracted. 

Because  the  computer  memory  was  insufficient  to  store  enough  array  scans  to 
obtain  meaningful  statistics  on  the  liftoff  height,  processing  had  to  be  done  in  real 
time  with  only  the  liftoff  height  being  stored.  As  a  sample  of  the  raw  data,  the  last 
64  scans  of  the  array  were  saved,  however.  This  technique  allowed  up  to  128K  seems 
to  be  acquired  for  each  flow  condition.  At  one  scan  every  4.4  ms,  this  amounted  to 
~  9j/2  minutes  of  data.  For  the  faster  flames,  switches  were  selected  on  the  circuit 
board  controlling  the  array  to  cause  only  half  of  the  pixels  to  be  output,  thereby 
increasing  the  data  acquisition  rate  to  one  scan  every  2.2  ms. 


2.2  Imaging  of  Fuel  and  Flame  Structures 

Experiments  to  obtain  planar  images  of  CH4  and  CH  were  conducted  at  the 
Combustion  Research  Facility,  Sandia  National  Laboratories  in  Livermore,  CA. 
Both  Raman  scattering  from  CH4  (fuel)  and  fluorescence  from  CH,  which  is  present 
only  in  the  reacting  zones,  were  measured  simultaneously  using  the  DIANA  pulsed 
dye  laser.  The  simultaneous  images  allowed  some  insight  to  be  gained  into  the 
relation  between  the  fuel  structures  and  the  reaction  zones. 

As  with  the  liftoff  height  measurements,  the  nozzle  used  in  the  imaging  exper¬ 
iments  was  designed  to  produce  fully  developed,  turbulent  pipe  flow  at  the  exit.  In 
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the  imaging  experiments,  however,  the  jet  emerged  from  a  50  mm  bluff  body,  rather 
than  being  tapered  on  the  outside.  To  form  the  laser  light  into  a  sheet,  a  multipass 
cell  was  used  (Namazian  et  al.  1988),  which  reflects  the  laser  beam  back  and  forth 
through  the  flame,  while  keeping  it  in  focus.  More  details  on  the  burner  and  mea¬ 
surement  technique  are  provided  in  Ch.  5,  where  these  experiments  are  discussed  in 
detail. 
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CHAPTER  3 

Mean  Liftoff  Heights 


Before  the  temporal  behavior  of  h  is  discussed  in  Ch.  4,  the  time-averaged  liftoff 
heights  axe  discussed  in  this  chapter.  Previous  results  of  other  researchers  are  first 
reviewed  in  Sec.  3.1.  These  findings  led  to  the  development  of  a  model  by  Dr.  R.  C. 
Miake-Lye  and  myself  which  implicates  the  large-scale  strain  rate  as  the  criterion 
which  determines  the  liftoff  height;  this  model  is  discussed  in  Sec.  3.2.  Results  on 
mean  liftoff  heights  from  the  current  work  are  presented  in  Sec.  3.3. 


3.1  Previous  Measurements 

Vanquickenbome  and  Van  Tiggelen  (1966)  measured  mean  liftoff  heights  of 
methane  using  nozzle  diameters  of  1.3,  1.8,  and  2.4  mm.  Using  such  small  nozzles 
limited  the  maximum  Reynolds  number  attainable  before  blowout  to  ~  1,200.  The 
most  extensive  set  of  measurements  of  mean  liftoff  heights  was  made  by  Kalghatgi 
(1984)  who  used  as  fuels  CH4,  C2H4,  H2,  and  propane  (CsHg).  He  used  tubes 
for  nozzles,  similar  to  the  current  experiments.  In  Table  3.1  are  listed  the  nozzle 
diameters  used  with  each  fuel  along  with  a  value  of  h/U0  obtained  by  fitting  a  line 
to  the  data  in  his  figures. 

An  important  experimental  finding  of  Kalghatgi  was  that  h  is  proportional  to 
U0 ■  The  slope  of  the  line  h/U0  depends  on  the  fuel,  but  is  essentially  independent 
of  nozzle  diameter,  once  h  is  sufficiently  far  downstream. 

Many  of  the  data  for  H2  were  obtained  with  choked  and  underexpanded  jets. 
For  these  data,  Kalghatgi  assumed  an  isentropic  expansion  to  atmospheric  pressure 
to  determine  an  effective  value  of  U0. 
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Table  3.1:  Liftoff  Height  Experiments  of  Kalghatgi  (1984) 

fuel 

nozzle  diameters 

h/Ua 

(mm) 

(ms) 

h2 

1.08,  1.74,  2.95,  4.06,  5.03,  6.1 

0.030 

ch4 

4.06,  5.03,  6.1,  8.3,  10.1 

2.4 

c3h8 

4.06,  5.03,  6.1,  8.3,  10.1 

1.9 

c2h4 

4.06,  6.1,  8.3 

0.67 

Kalghatgi  interpreted  his  experimental  results  by  arguing  that  the  speed  of  the 
gas  just  upstream  of  the  stabilization  position  is  equal  to  a  turbulent  flame  speed, 
which  depends  on  both  the  fuel’s  laminar  flame  speed  and  the  turbulence  properties. 
An  alternative  explanation  of  the  results  is  presented  in  the  following  section. 


3.2  Strain-Rate  Model 

The  strain-rate  model  provides  a  plausible  explanation  for  the  experimental 
observations  that,  for  a  given  fuel,  the  mean  liftoff  height  is  proportional  to  the 
fuel  velocity  exiting  the  nozzle,  independent  of  the  nozzle  diameter.  The  model 
was  initially  developed  by  Dr.  R.  C.  Miake-Lye.  He  and  I  worked  out  the  details 
and  conducted  experiments  to  test  one  aspect  of  the  model — the  effect  on  the  liftoff 
height  of  adding  some  air  to  the  fuel  (Miake-Lye  and  Hammer  1989). 

The  premise  of  the  model  is  that  the  flame  stabilizes  where  the  strain  rate  is 
Just  low  enough  to  allow  burning  without  extinction,  t.e.,  at  the  critical  strain  rate  of 
the  fuel,  <rcr.  Considering  the  large  structure  observed  in  non-reacting  turbulent  jets 
( e.g .,  Dimotakis  et  al.  1983),  the  relevant  strain  rate  is  taken  to  be  that  between 
large  turbulent  structures.  This  strain  rate  must  scale  with  the  local  large-scale 
variables  Uc i  and  S,  where  Uc\  is  the  local  centerline  velocity  and  6  is  the  local  jet 
width,  taken  to  be  the  full  visual  width  0.44  x. 


Because  hydrocarbon  fuels  have  reaction  rates  that  are  very  sensitive  to  com¬ 
position,  with  the  highest  reaction  rates  near  stoichiometric  conditions,  the  highest 
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strain  rates  that  the  flame  can  sustain  also  occur  near  stoichiometric  conditions. 
Hence,  this  model  also  requires  that  the  local  mass  fraction  of  fuel  Y  be  at  the 
stoichiometric  value  Y^t. 


In  the  derivation  that  follows,  as  well  as  throughout  this  thesis,  two  often-used 
quantities  are  the  equivalent  source  velocity  Us  and  source  diameter  <4-  They  are 
defined  in  terms  of  the  momentum  flux  J0  and  mass  flux  m0  exiting  the  nozzle  by 
the  solution  of  the  equations 


m0 


Jo 


(3.1) 


The  quantitative  predictions  of  this  model  are  derived  by  considering  the  simi¬ 
larity  laws  of  non-reacting,  turbulent  jets.  As  a  result,  the  model  is  applicable  only 
to  flames  stabilized  a  sufficient  distance  from  the  nozzle  exit.  In  the  jet  far  field, 
the  mean  velocity  is  described  by 

a  =  C.  V,  ^  /.(,)  ,  (3.2) 

X 

and  the  mean  concentration  by 

Y  =  CyY^Mi)  ■  (3-3) 

X 

The  large  scale  strain  rate  oi  is  given  by 

°L  ~  •  (3-4) 

Noting  that  6  ~  x  and  using  Eq.  3.2  results  in  the  similarity  law 

»i  =  c,  ^  u-n) .  (3.5) 

In  these  equations,  C„,  Cy,  and  C„  are  constants  determined  from  experiments,  Y0 
is  the  initial  mass  fraction  of  fuel  in  the  nozzle  fluid  with,  the  balance  being  air, 
t?  is  the  non-dimensional  radial  coordinate,  r/x,  and  /u(r?),  /y(r/),  and  fa{rj)  are 
similarity  functions,  also  determined  from  experiments,  where  possible. 
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The  radial  flame  location  will  fluctuate,  approximately  tracking  the  position 
where  instantaneously  Y  =  Y'si;  this  will  be  distributed  around  the  point  where 

Y  —  V;t.  For  simplicity  in  describing  the  strain- rate  model,  it  is  assumed  that 

Y  =  Yst.  Additionally  it  is  assumed  that  o i  is  equal  to  the  critical  strain  rate  for 
the  fuel,  <7Cr-  A  more  detailed  consideration  of  the  statistics  of  the  strain-rate  and 
concentration  fields  would  not  alter  the  qualitative  predictions  of  the  model. 

Solving  for  ds/x  from  Eq.  3.3  with  Y  =  Fst, 


ds  =  1  Y,t  1 

x  Cy  Y0  /y(f?) 


(3.6) 


putting  this  into  Eq.  3.5  and  setting  x  =  h  results  in  the  expression  for  the  mean 
liftoff  height 


-  ca  y;t  us  Uv) 
cY  y0  act  /y(r?)  ’ 


(3.7) 


where  r/  is  given  by  Eq.  3.3  as 


n 


(3.8) 


If  the  the  strain-rate  profile  and  concentration  profile  are  the  same,  j.  e.,  if  fai’n)  = 
/y(^)  for  T) min.  <  rj  <  r/max  ,  then  the  strain-rate  model  predicts  a  linear  relation 
between  h  and  Us, 


(3-9) 


Without  this  assumption  about  the  strain  rate  and  concentration  profiles,  the 
model’s  predictions  can  be  written  as 


Yst  h_  =  C*_ 

Y0  (4  Cy 


»!))• 


where 


/srm(^)  =  fa(r))/fY{V )  ■ 


(3.11) 


This  equation  implicitly  expresses  the  non-dimensional  liftoff  height  ( Yst/Y0)h/ds 
as  a  function  of  a  non-dimensional  velocity  (Yst/Y0)2(l/aCT)Us/ds. 
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FlG.  3.1  Non-dimensional  liftoff  height  vs.  velocity  as  predicted  by  the  strain-rate 
model,  when  ctl  is  approximated  by  the  largest  principal  compressive 
(—dill)  or  extensive  (dj)  strain  rate  of  the  mean  flow  field. 

To  solve  for  (Yst/Y0)h/ds  requires  knowledge  of  the  two  functions  fyiv)  and 
fa(  17).  Chen  and  Rodi  (1980)  recommend 

frW  =  e-c'y'  (3.12) 

with  C'Y  =  57.  The  relevant  strain  rate  for  flame  stabilization  is  that  produced  by 
the  large  structures  of  the  turbulence.  This  is  difficult  to  measure,  and  therefore 
only  an  approximate  treatment  is  possible.  One  choice  is  to  use  the  maximum 
principal  strain  rate  obtained  from  the  mean  velocity  field  for  ai,  even  though  this 
may  not  be  a  very  good  approximation  to  the  relevant  strain  rate  produced  by 
the  turbulence.  Using  the  expression  from  Chen  and  Rodi  (1980)  for  the  velocity 
function, 

f.M  =  e" ° (3.13) 

with  C'u  =  94,  the  principal  strain  rates  can  be  calculated,  as  discussed  in  Sec.  A.2. 
Choosing  either  the  maximum  compressive  strain  rate  (— ^m)  or  the  maximum  ex- 
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tensive  strain  rate  ( d\ )  allows  Eq.  3.10  to  be  solved  for  ( Yst/Y0)h/ds  as  a  function 
of  ( Yst/Yo)2Us/ds .  The  solutions  for  both  choices  are  shown  in  Fig.  3.1.  Below 
the  bottom  dotted  line,  ( Yst/Y'0)h/ds  <  1.2,  which  corresponds  to  h/ds  <  20  for 
Yat  /Y0  =  0.06.  This  may  be  too  close  to  the  nozzle  origin  for  the  fax- field  equa¬ 
tions  to  apply.  (The  kink  in  the  a l  =  —dm  curve  is  a  result  of  a  change  in  the 
strain-rate  component  of  maximum  compression  to  the  azimuthal  component  for 
low  ( Yst/Y0)h/ds ,  or  large  r?.)  Above  the  top  dotted  line,  (Ysl/Y0)h/ds  >  2.6,  where 
flames  cannot  be  stabilized.  Between  the  two  dotted  lines,  both  curves  show  an 
approximately  linear  relationship,  with  a  negative  virtual  origin  on  the  liftoff  height 
axis.  This  is  an  indication  that  the  strain-rate  model  does  a  reasonable  job  of  pre¬ 
dicting  a  linear  relation  between  h  and  Us.  Because  the  mean  strain- rate  field  was 
used  for  approximating  fa(v)  while  the  turbulence  will  cause  the  true  variation  of 
a i  with  T]  will  be  different,  the  model’s  actual  predictions  may  differ  from  Fig.  3.1. 
As  long  as  the  effective  strain  rates  increase  as  rj  is  reduced  (moving  toward  the  jet 
centerline)  in  approximately  the  same  way  as  the  concentration  profile,  the  strain 
rate  model  predicts  a  linear  h  vs.  Us  relation. 


Table  3.2:  Strain-Rate  Model  Predictions 

fuel 

’Ho  —  h/Yst  Us 
(ms) 

l/<7ext 

(ms) 

'Ho  ^ext 

k/S2 

(ms) 

Tlo  S2/k 

ch4 

43.2 

3.13 

13.8 

0.698 

61.9 

c2h4 

10.4 

- 

- 

0.174 

59.8 

c3h8 

31.2 

2.86 

10.9 

0.503 

62.0 

h2 

1.04 

0.0769 

13.5 

0.012 

86.0 

In  addition  to  predicting  a  linear  relation  between  h  and  Us,  a  correct  model 
for  liftoff  must  predict  the  observed  slope  of  the  h/Us  line  for  different  fuels.  A 
comparison  of  experimental  results  with  model  predictions  is  made  in  Table  3.2. 
An  effective  chemical  time  T|0  can  be  defined  for  each  fuel  with  the  equation 

h  —  Ysi  7i0  Us  , 


(3.14) 
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where  h/Ua  is  taken  from  experimental  data.  In  the  second  column  of  Table  3.2 
Kalghatgi’s  experimental  liftoff  results  have  been  used  to  estimate  Tj0.  These  values 
can  be  compared  to  properties  of  the  fuels  which  provide  estimates  of  the  critical 
strain  rates.  In  the  third  column,  the  extinction  limits  of  laminar,  counterflow 
diffusion  flames  are  listed.  Values  for  CH4  and  C3H8  were  taken  from  Tsuji  and 
Yamaoka  (1969),  and  the  value  for  H2  was  taken  from  Dixon- Lewis  et  al.  (1986). 
Although  the  combustion  in  a  lifted  turbulent  jet  flame  does  not  take  place  at  the 
interface  between  pure  fuel  and  pure  air  because  of  premixing  upstream  of  the  flame 
base,  <7ext  is  used  as  an  approximate  measure  of  crcr.  Flames  with  partial  premixing 
have  higher  extinction  limits  than  pure  diffusion  flames.  For  instance,  Mastorakos 
et  al.  (1992)  found  that  the  extinction  limit  of  counterflow  diffusion  flames  of  20 
%  CH4  /  80  %  air  is  <rext  =  600  s-1,  compared  with  350  s-1  for  pure  CH4.  The 
fourth  column  of  Table  3.2  is  the  ratio  of  the  second  and  third  columns,  aiding 
comparison  of  the  model’s  predictions  for  the  different  fuels.  That  the  numbers 
in  the  fourth  column  are  reasonably  close  to  one  another  provides  support  for  the 
strain-rate  model.  Using  the  definition  of  ri0  from  Eq.  3.14  along  with  Eqs.  3.2, 
3.13,  and  3.3,  it  can  be  shown  that  air\0  k  7.5  for  Ysth/ds  =  1.5  to  2.5,  where 
d\  is  calculated  on  the  contour  Y  =  Yst.  That  the  values  in  the  fourth  column  of 
Table  3.2  are  in  the  neighborhood  of  7.5  indicates  that  the  strain  rates  at  the  flame 
base,  as  estimated  by  <7j,  are  close  to  (but  slightly  less  them)  <rext.  It  is  likely  that  the 
turbulence  increases  the  strain  rates  above  that  obtained  from  the  mean  velocity 
field.  However,  premixing  of  fuel  and  air  will  increase  the  maximum  sustainable 
strain  rates. 

As  an  alternative  estimate  for  the  critical  strain  rate,  52/k  is  tabulated  in  the 
fifth  column  of  Table  3.2,  where  k  is  the  thermal  diffusivity  and  5  is  the  laminar 
flame  speed  Because  k/S2  is  a  relevant  chemical  time  in  premixed  flame  propaga¬ 
tion  ( e.g .,  Glassman  1977),  its  inverse,  which  has  dimensions  of  a  strain  rate,  can 
be  thought  of  as  a  characteristic  maximum  strain  rate  of  the  fuel.  In  the  table,  k  is 
defined  as  fc//pst  CPiSt  where  kf  is  the  thermal  conductivity  of  the  reacted  mixture, 
here  estimated  as  that  of  nitrogen  at  the  adiabatic  flame  temperature,  pst  is  the 
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density  of  the  unburned,  stoichiometric  mixture,  and  CP)St  is  the  specific  heat  of  the 
unbumed  mixture.  As  with  the  estimates  of  <rcr  based  on  <7ext,  k/S 2  provides  a  rea¬ 
sonable  collapse  of  the  liftoff  data  for  the  fuels,  as  seen  by  the  reasonable  agreement 
between  the  numbers  in  the  sixth  column  of  Table  3.2. 


3.3  Current  Measurements 

Measurements  of  the  time- varying  liftoff  height  were  made  as  described  in  Ch.  2 
using  a  linear  photodiode  array.  From  these  measurements  mean  values  were  easily 
obtained.  Two  sets  of  measurements  are  reported  in  this  chapter.  First,  liftoff 
heights  of  pure  fuels  for  various  nozzle  diameters  and  exit  velocities  are  reported, 
then  measurements  of  h  with  premixing  of  air  with  the  fuel  are  presented. 


3.3.1  Details  of  the  Measurement  Technique 

In  making  the  measurements  reported  in  this  chapter,  for  each  nozzle  used, 
the  array  and  camera  lens  were  positioned  at  a  vertical  height  equal  to  the  liftoff 
height  in  the  middle  of  the  range  of  h  to  be  measured,  specifically,  h  for  which 
Ysih/ds  =  2.2.  The  horizontal  distance  from  the  nozzle  to  the  array  was  chosen 
such  that  the  flame  width  just  filled  the  lateral  extent  of  the  array  when  the  flame 
was  at  this  middle  height,  as  illustrated  in  Fig.  3.2.  Because  the  optics  were  not 
moved  as  the  liftoff  height  was  changed,  flames  with  h  greater  than  the  middle  value 
had  images  somewhat  wider  than  the  array,  and  therefore  the  left  and  right  sides  of 
the  image  were  chopped  off.  This  was  found  to  have  a  negligible  effect  on  the  mean 
value  measured,  however,  as  documented  in  Sec.  D.6.  For  each  condition,  8192  data 
points  were  acquired  at  a  repetition  time  of  4.4  ms. 
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FlG.  3.2  Schematic  showing  the  image  of  the  flame  on  the  linear  array.  Pixels  are 
50  times  as  wide  as  their  vertical  spacing. 

3.3.2  Pure  Fuels 

The  four  fuels  used  in  measurements  of  mean  liftoff  heights  were  natural  gas, 
C2H6,  C2H4,  and  C2H2.  Properties  of  these  fuels  (molecular  weight  MW,  Y^t,  and 
adiabatic  flame  temperature  rise  AT/)  are  listed  in  Table  3.3.  The  properties  of 
H2  are  also  listed  because  that  fuel  was  used  in  experiments  described  in  Ch.  4. 
The  properties  of  natural  gas  were  calculated  from  a  sample  analysis  supplied  by 
the  gas  company;  the  composition  is  given  on  p.  4.  The  flame  temperatures  and 
stoichiometric  mass  fractions  for  all  the  fuels  were  taken  from  Becker  and  Liang 
(1978).  The  three  fuels  C2H6,  C2H4,  and  C2H2  were  chosen  in  part  because  their 
molecular  weights  are  close  to  that  of  air,  thereby  minimizing  the  effect  of  different 
jet  and  ambient  densities;  natural  gas  was  used  because  of  its  convenience. 

The  flow  conditions  for  the  measurement  of  mean  liftoff  heights  are  listed  in 
Table  3.4.  For  each  fuel,  three  nozzle  diameters  were  used,  except  for  C2H2  for 
which  only  two  nozzles  were  used.  The  corresponding  minimum  and  maximum 


Table  3.2:  Properties  of  Fuels  Used 

fuel 

symbol 

MW 

1st 

ATf 

(K) 

natural  gas 

see  p.  4 

17.1 

0.0582 

1925 

ethane 

c2h6 

30.0 

0.0584 

1959 

ethylene 

c2h4 

28.0 

0.0633 

2069 

acetylene 

c2h2 

26.0 

0.0700 

2239 

hydrogen 

h2 

2.0 

0.0284 

2117 

velocities  for  each  fuel/nozzle  combination  are  listed,  along  with  the  corresponding 
minimum  and  maximum  liftoff  heights.  All  of  the  measurements  discussed  in  this 
chapter  are  tabulated  in  Appendix  B. 


Table  3.4:  Conditions  of  h  Measurements  Without  Air  Premixing 


fuel 

d 

min.  Us 

max.  Us 

min.  h 

max.  h 

(mm) 

(m/s) 

(m/s) 

(mm) 

(mm) 

n.g. 

4.61 

33 

56 

89 

164 

6.25 

42 

74 

133 

229 

7.73 

39 

90 

137 

272 

c2h6 

3.12 

46 

71 

74 

124 

3.88 

48 

93 

89 

180 

4.61 

58 

111 

115 

210 

c2h4 

2.16 

73 

127 

44 

91 

3.12 

95 

181 

66 

135 

3.88 

111 

200 

82 

143 

c2h2 

1.02 

173 

240 

20 

37 

1.55 

226 

266 

32 

40 

Mean  liftoff  heights  vs.  source  velocity  Us  are  shown  in  Fig.  3.3  for  natural 
gas  using  three  different  nozzle  diameters.  Details  of  the  method  used  to  calculate 
Ua  and  ds  can  be  found  in  Sec.  A.l.  From  the  figure,  it  can  be  seen  that  h  is 
approximately  proportional  to  £/s,  although  there  is  some  variation  with  nozzle 
diameter. 
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FlG.  3.3  Mean  liftoff  heights  for  natural  gas,  d  —  4.61,  6.25,  and  7.73  mm. 

To  see  more  clearly  the  effect  of  nozzle  diameter  it  is  useful  to  plot  the  data 
non-dimensionally.  Recalling  the  discussion  in  Sec.  3.2,  the  proper  non-dimensional 
variable  for  the  liftoff  height  is  Ysth/ds  (for  Y0  =  1).  For  a  given  fuel  with  known  Fst 
and  MW,  h/ds  and  h/d ,  which  are  relevant  to  determining  whether  or  not  h  is  in  the 
far  field  of  the  jet,  are  linearly  related  to  Ysth/ds.  The  proper  variable  to  use  for  the 
velocity  is  less  obvious.  Following  the  strain  rate  model,  a  possibility  is  Y//tTxUs/d9, 
where  tx  is  a  chemical  time  of  the  fuel.  In  the  strain-rate  model  rx  is  proportional 
to  l/<rcr.  Other  models  use  Y^rxUs/ds  as  a  non-dimensional  velocity,  with  different 
choices  of  rx,  and  possibly  with  Yst  replaced  with  another  concentration,  such  as 
Y  at  the  maximum  laminar  flame  speed.  To  avoid  choosing  a  chemical  time  rx, 
Yith/ds  is  plotted  against  Y9tUs/ds  in  Fig.  3.4,  making  the  slope  of  each  line  on  the 
plot  have  dimensions  of  ms.  The  slope  of  a  line  is  then  the  effective  value  of  tig, 
defined  as  h/Usd9  in  Eq.  3.14.  Values  of  h/d9  and  h/d  are  labeled  along  the  right 
side  of  the  plot. 

FYom  Fig.  3.4  it  can  be  seen  that  changing  the  nozzle  diameter  does  have  a 
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FlG.  3.4  Normalized  mean  liftoff  heights  vs.  normalized  velocity  for  natural  gas. 
Values  of  h/ds  and  h/d  are  shown  along  the  right  side  of  the  graph. 


slight  effect  on  the  normalized  liftoff  height.  As  d  is  increased  from  4.61  mm  to 
6.25  mm,  Ysth/ds  at  a  constant  Y^Us/ds  increases.  However,  there  is  no  further 
increase  in  going  to  the  7.73  mm  nozzle.  It  is  possible  that  the  changing  Reynolds 
number  is  responsible  for  this  behavior  in  which  the  normalized  liftoff  height  first 
increases,  and  then  remains  constant  with  further  increase  in  the  Reynolds  number. 
The  Reynolds  number  in  the  fax  field  of  the  jet  is  defined  as 

Re oo  =  ,  (3.15) 

Voo 

and  is  approximately  equal  to  the  Reynolds  number  just  below  the  flame  base  Re*, 
which  is  defined  in  Sec.  A. 5.  The  ranges  of  Re o©  corresponding  to  the  data  points 
in  Figs.  3.3  and  3.4  are:  for  the  4.61  mm  nozzle,  7,300  to  12,000;  for  the  6.25  mm 
nozzle,  13,000  to  22,000;  and  for  the  7.73  mm  nozzle,  14,000  to  34,000.  The  density 
in  the  flame  is  lower  than  p <*,  and  the  viscosity  is  higher  them  p<x>,  making  the 
effective  Reynolds  number  in  the  flame  considerably  lower  than  Re Therefore,  it 
is  plausible  that  variation  of  Reynolds  number  could  have  an  effect  on  these  flames 
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even  with  Re^  >  15,000. 

As  larger  nozzle  diameters  are  used,  the  required  velocity  for  a  given  non- 
dimensional  liftoff  height  increases,  and  eventually,  compressibility  effects  become 
important,  at  least  in  some  region  of  the  jet  near  the  nozzle  exit.  Because  the  speed 
of  sound  in  the  fuel  a0  is  in  general  not  equal  to  that  in  the  air  a^.,  there  is  some 
ambiguity  in  the  exact  definition  of  a  Mach  number.  Here,  the  Mach  number  of 
the  jet  will  be  defined  as  M„  =  f/s/«oo-  Because  most  of  the  fuels  studied  (C2H6, 
C2H4,  and  C2H2)  have  speeds  of  sound  close  to  that  of  air,  the  precise  definition  of 
Mg  will  not  qualitatively  change  the  arguments  presented  here. 

For  the  natural  gas  results  shown  in  Fig.  3.4  it  is  unlikely  that  the  change  in 
Mg  is  responsible  for  the  variation  with  nozzle  diameter,  since  the  maximum  value 
of  Mg  achieved  in  natural  gas  experiments  was  0.26.  However,  this  may  not  be  true 
for  all  of  the  fuels,  as  is  discussed  below. 


h/ds 


Fig.  3.5  Normalized  liftoff  height  vs.  normalized  velocity  for  C2H6,  d  =  3.12,  3.88, 
and  4.92  mm. 
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The  second  fuel  used  in  these  experiments  was  C2H6.  The  mean  liftoff  heights 
are  plotted  in  Fig.  3.5,  where  Us  and  h  have  been  normalized  as  was  done  with 
natural  gas  in  Fig.  3.4.  Because  ds  ~  d  for  C2H6,  only  values  of  h/ds  are  labeled 
along  the  right  side  of  the  plot.  The  variation  with  nozzle  diameter  (or  Reynolds 
number)  is  slight;  the  ranges  of  Re oo  for  the  C2H6  data  plotted  are:  for  the  3.12  mm 
nozzle,  9,000  to  14,000;  for  the  3.88  mm  nozzle,  12,000  to  23,000;  and  for  the 
4.92  mm  nozzle,  18,000  to  35,000.  Although  the  data  from  the  last  two  nozzles 
do  not  collapse  completely,  there  is  less  variation  in  liftoff  height  in  changing  from 
the  3.88  mm  nozzle  to  the  4.92  mm  nozzle,  than  in  changing  from  the  3.12  mm 
nozzle  to  the  3.88  mm  nozzle.  As  with  the  natural  gas  experiments,  compressibility 
effects  should  be  insignificant,  as  the  highest  Mach  number  realized  for  C2H6  was 
Ma  =  0.32. 


h/ds 


Fig.  3.6  Normalized  liftoff  height  vs.  normalized  velocity  for  C2H4,  d  =  2.16,  3.12, 
and  3.88  mm. 


C2H4  reacts  faster  than  C2H6  and  consequently  higher  velocities  are  needed  to 
achieve  the  same  liftoff  height.  Figure  3.6  is  a  plot  of  normalized  liftoff  height  vs. 
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velocity  for  nozzle  diameters  of  2.16,  3.12,  and  3.88  mm.  Again,  as  d  is  increased, 
Ysth/ds  increases  slightly  at  a  fixed  value  of  Y^Us/ds,  but  the  effect  is  smellier 
in  going  from  the  3.12  mm  to  the  3.88  mm  nozzle  than  going  from  the  2.16  mm 
to  the  3.12  mm  nozzle.  The  ranges  of  Re for  the  various  nozzles  are:  for  the 
2.16  mm  nozzle,  9,700  to  17,000;  for  the  3.12  mm  nozzle,  18,000  to  35,000;  and 
for  the  3.88  mm  nozzle,  27,000  to  48,000.  This  range  of  Rex  spans  the  range 
explored  with  the  natural  gas  and  C2H6  flames  in  which  no  Reynolds  number  effect 
was  found  above  ~  13,000  for  natural  gas,  and  little  Reynolds  number  effect  found 
above  20,000  for  C2H6.  Therefore,  the  continued  increase  in  Ysth/ds  with  Re^  for 
C2H4  requires  another  explanation. 

Possibly  compressibility  effects  are  beginning  to  affect  the  data.  The  Mach 
number  ranges  for  the  nozzles  are:  for  the  2.16  mm  nozzle,  Ms  =  0.21  to  0.37;  for 
the  3.12  mm  nozzle,  0.28  to  0.52;  and  for  the  3.88  mm  nozzle,  0.32  to  0.58.  It  is 
possible  that  close  to  the  nozzle  where  velocities  are  high,  reduced  entrainment  in 
the  mixing  layer  (Papamoschou  1988),  caused  by  the  high  Mach  numbers,  effectively 
acts  to  create  a  virtual  origin  downstream  of  the  nozzle  exit.  Assuming  that  this 
virtual  origin  is  the  same  for  all  properties  of  the  jet  (concentration,  velocity,  strain 
rates,  etc.),  this  would  cause  Ysth/ds ,  for  which  no  virtual  origin  effect  has  been 
considered,  to  increase. 

Finally,  mean  liftoff  height  data  for  C2H2  are  shown  in  Fig.  3.7.  The  achievable 
conditions  for  C2H2  were  limited  by  the  allowable  delivery  pressure  for  that  gas 
of  15  psig.  When  the  data  for  the  1.02  mm  nozzle  are  extrapolated  to  Us  =  0,  a 
fine  drawn  through  the  data  intersects  the  vertical  axis  at  a  large,  negative  non- 
dimensional  distance.  Although  still  having  a  negative  intercept  on  the  vertical 
axis,  the  intercept  is  smaller  for  the  1.55  mm  nozzle.  As  with  ethylene,  an  increase 
in  diameter  results  in  an  increase  in  Ysth/da.  The  ranges  of  Re in  Fig.  3.7  Eire:  for 
the  1.02  mm  nozzle,  11,000  to  15,000;  and  for  the  1.55  mm  nozzle,  20,000  to  24,000. 
The  Mach  numbers  are  the  highest  of  all  the  fuels  used  in  this  study:  0.51  to  0.69 
for  the  1.02  mm  nozzle  and  0.64  to  0.76  for  the  1.55  mm  nozzle. 
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h/ds 


FlG.  3.7  Normalized  liftoff  height  vs.  normalized  velocity  for  C2H2,  d  =  1.02  and 
1.55  mm. 

To  better  illustrate  the  effects  of  Re <*,  and  Ms  on  mean  liftoff  heights,  data 
from  all  four  fuels  have  been  plotted  on  a  single  graph,  Fig.  3.8.  This  was  done 
by  selecting  a  single  value  of  Ysth/ds  =  2.0  and  plotting  the  corresponding  value 
of  Ys\Us/ds  vs.  Re  oo.  Because  the  four  fuels  have  vastly  different  effective  chemical 
times  ti0,  each  curve  has  been  normalized  by  the  value  of  Y£Us/ds  at  Re 00  =  15,000, 
to  bring  the  curves  onto  the  same  scale  and  force  them  all  to  pass  through  the  point 
(15,000,  1).  With  this  normalization,  the  value  of  Fst  for  each  fuel  does  not  affect 
the  vertical  axis,  and  has  been  left  off  the  vertical  axis  label.  One  point  for  C2H2 
had  to  be  extrapolated  from  the  data  because  the  highest  liftoff  height  reached  for 
the  1.55  mm  nozzle  was  Ysth/ds  =  1.9. 

As  the  Reynolds  number  is  increased,  the  velocity  needed  to  reach  Ysth/ds  =  2.0 
decreases.  All  the  fuels  do  not  lie  on  a  single  curve,  however,  indicating  that  Re  <*,  is 
not  the  only  parameter  causing  this  decrease.  The  curves  for  natural  gas  and  C2H6 
level  out  at  Re =  15,000  and  18,000,  respectively,  while  those  of  C2H4  and  C2H2 
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FlG.  3.8  Normalized  velocity  YatUa/da  at  Yath/da  =  2.0  vs.  Reynolds  number  Re <», 
for  all  four  fuels.  The  curves  have  been  normalized  to  pass  through  the 
point  (15,000,  1)  to  allow  all  fuels  to  be  displayed  on  one  plot. 

continue  to  fall.  It  is  possible  that  this  continued  decrease  for  the  latter  two  fuels 
is  actually  caused  by  an  increase  in  Mach  number,  which  presumably  decreases  the 
mixing  in  the  region  near  the  nozzle,  creating  a  virtual  origin  downstream  of  the 
nozzle  exit.  It  is  also  possible,  however,  that  these  two  fuels  simply  have  a  different 
dependence  on  Re <»,  which  results  from  details  of  the  chemical  reaction  mechanism. 

As  discussed  in  Dimotakis  (1993),  there  is  evidence  of  a  transition  in  the  be¬ 
havior  of  turbulent  jets  in  tne  Reynolds  number  range  of  the  current  experiments. 
Gilbrech  (1991)  measured  flame  lengths  of  low-heat-release  gas  flames.  He  fit  his 
results  to  an  equation  of  the  form 

^  =  A<f>+B,  (3.16) 

<h 

where  L  is  the  mean  flame  length  and  <f>  is  the  stoichiometric  air  to  fuel  ratio, 
1/Ht  —  1.  The  stoichiometric  coefficient  A  and  the  normalized  virtual  origin  B 
were  both  found  to  depend  on  Reynolds  number.  In  particular,  A  wa s  found  to 
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decrease  as  Rea 0  increases  from  10,000  to  20,000,  above  which  the  Reynolds  number 
dependence  was  found  to  be  less,  if  any,  at  least  up  to  150,000. 

For  natural  gas  and  C2H6,  the  change  in  liftoff  height  behavior  shown  in  Fig.  3.8 
occurs  at  approximately  the  same  value  of  Re^  as  the  change  in  A  shown  in 
Gilbrech’s  Fig.  4.8  (Fig.  7  in  Dimotakis  1993).  The  other  two  fuels  plotted  in 
Fig.  3.8  may  have  compressibility  effects,  as  discussed  above.  (In  the  flame-length 
study,  the  Mach  number  did  no  change  with  Re0 a  because  Re oo  was  controlled  by 
adjusting  the  pressure  in  the  reaction  chamber.)  Dimotakis  notes  that  many  other 
flows,  in  particular  the  shear  layer,  have  a  transition  in  this  same  range  of  Re^. 

A  property  of  the  jet  which  changes  with  Re oo  is  the  strain  rate  at  the  smallest 
scale  of  the  flow.  The  strain  rate  at  the  Kolmogorov  scale,  ok,  is  proportional  to 
01  RelJ? .  This  increase  in  a ^  with  Re^  may  have  the  effect  of  destabilizing  the 
flame,  causing  it  to  move  to  a  slightly  higher  position,  where  both  ai  and  (tk  are 
lower.  Why  the  effect  would  diminish  above  Re  oo  «  20,000  is  not  clear,  however. 

To  get  an  idea  of  the  liftoff  behavior  for  even  larger  Reynolds  numbers,  the 
results  of  Birch  and  Hargrave  (1989)  can  be  examined.  They  studied  very  large 
natural  gas  flames  in  which  the  jet  was  choked  and  underexpanded,  using  pressures 
up  to  28  bar.  By  assuming  an  isentropic  expansion  to  ambient  pressure  from  sonic 
conditions  at  the  nozzle  exit,  an  equivalent  source  velocity  can  be  calculated,  and 
an  equivalent  source  diameter  can  be  similarly  calculated.  (For  isentropic  flow 
equations,  see,  e.g .,  Liepmann  and  Roshko  1957,  pp.  53,  126.)  In  Fig.  3.9a  their 
measured  mean  liftoff  heights  are  plotted  vs.  the  calculated  source  velocity  Us  for 
nozzle  diameters  of  50.6  mm  and  111.2  mm,  along  with  data  from  the  current 
experiments  using  a  6.25  mm  nozzle.  h/Us  is  ~  3.2  ms  for  the  50.6  mm  nozzle,  and 
ranges  from  ~  3.4  to  4.4  ms  for  the  111.2  mm  nozzle,  all  reasonably  close  to  the 
value  of  3.1  ms  measured  in  the  current  experiments  (Fig.  3.3).  Re ranges  from 
1.2  x  106  to  6.2  x  106  for  the  50.6  mm  nozzle,  and  from  1.6  x  106  to  7.5  x  106  for 
the  111.2  mm  nozzle,  more  than  200  times  as  large  as  in  the  current  experiments. 

The  data  are  plotted  again  in  Fig.  3.9b,  using  the  same  normalization  as  in 
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Fig.  3.9  Mean  liftoff  heights  for  natural  gas,  from  Birch  and  Hargrave  (1989),  along 
with  the  current  results  for  a  6.25  mm  nozzle,  (a)  dimensional,  (b)  nor¬ 
malized. 

Fig.  3.4.  Because  the  calculated  value  of  ds  changes  as  the  stagnation  pressure  is 
changed,  the  shapes  of  the  curves  are  different  from  Fig.  3.9a.  There  are  multiple 
values  of  the  abscissa,  because  Ua/ck  reaches  a  maximum  at  an  intermediate  pres¬ 
sure  and  then  decreases.  The  left-most  point  on  Fig.  3.9b  for  the  50.6  mm  nozzle 
corresponds  to  the  right-most  point  in  Fig.  3.9a.  For  the  111.2  mm  nozzle,  however, 
the  left-most  points  of  each  plot  correspond.  As  long  as  the  flames  are  out  of  the 
near  field  of  the  jet,  the  value  of  h/Us  is  approximately  the  same  as  for  the  much 
smaller  nozzles  of  the  current  experiments. 


32 


3.3.3  Fuels  Premixed  with  Air 

Experiments  were  conducted  in  which  the  fuel  was  premixed  with  air  prior  to 
exiting  the  nozzle.  The  strain-rate  model  described  in  Sec.  3.2  predicts  that  when 
air  is  added  to  the  nozzle  fluid  reducing  the  initial  mass  fraction  of  fuel  in  the  nozzle 
fluid  1'0,  h  oc  ( 1  /  Y0 )  for  a  given  fuel  and  source  velocity.  These  experiments  were 
conducted  largely  to  test  this  prediction.  Several  other  models  of  liftoff  which  use 
far-field  scaling  of  jets  make  this  same  prediction,  as  summarized  in  Hammer  (1990). 

Miake-Lye  and  Hammer  (1989)  conducted  experiments  in  which  the  nozzle 
diameter  was  kept  constant  and  Ya  was  varied.  In  that  study,  as  Y0  was  reduced, 
the  slope  of  the  liftoff  curve  exceeded  predictions.  Because  those  experiments  were 
conducted  with  the  nozzle  diameter  fixed  while  Y0  was  varied,  the  Reynolds  number 
dropped  as  Ya  was  decreased. 

In  the  current  experiments,  however,  the  nozzle  diameter  was  changed  system¬ 
atically  to  keep  the  Reynolds  number  range  of  the  experiments  the  same  at  each 
value  of  Y0.  For  each  fuel,  this  meant  keeping  Y0ds  constant  so  that  at  a  given 
value  of  (Yst/Y0)2Us/ds  (at  which  ( Yst/Y0)h/ds  is  predicted  to  be  fixed)  Re^  will 
be  fixed,  for  any  value  of  Y„.  This  is  more  clearly  seen  by  examining  the  equation 

Re<x,  =  (Usdj/uoo  =  (i/y9t)2[(nt/n)2^s/d9][F0ds]2/^oo .  (3.17) 

The  specific  conditions  that  were  studied  are  listed  in  Table  3.5.  Experiments 
were  done  with  natural  gas,  C2H6,  and  C2H4,  using  four  nozzles  and  four  corre¬ 
sponding  values  of  Y0  for  each  fuel. 

Data  from  measurements  using  natural  gas  are  plotted  in  Fig.  3.10  for  several 
values  of  Ya.  For  all  the  curves  plotted,  the  range  of  Re0 0  is  ~  6,000  to  12,000, 
which  is  too  low  to  rule  out  effects  of  Re o©  on  the  data.  The  Mach  numbers  Ms  are 
all  less  than  0.2,  which  is  sufficiently  low  to  have  negligible  compressibility  effects. 
The  strain-rate  model  predicts  that  all  of  the  curves  would  lie  on  top  of  one  another. 
However,  as  more  air  is  added  to  the  nozzle  fluid,  the  slope  of  the  curve  increases, 
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Table  3.5:  Conditions  of  h  Measurements  with  Air  Premixing 


fuel 

d 

(mm) 

Y0 

min.  Us 
(m/s) 

max.  Us 
(m/s) 

min.  h 

(mm) 

max.  h 

(mm) 

n.g. 

4.61 

1 

33 

56 

89 

164 

4.92 

0.92 

25 

51 

75 

169 

6.25 

0.69 

19 

38 

73 

181 

7.73 

0.53 

15 

30 

68 

189 

c2h6 

3.88 

1 

47 

89 

91 

170 

4.92 

0.79 

37 

68 

89 

168 

6.25 

0.63 

30 

55 

88 

186 

7.73 

0.50 

24 

41 

77 

168 

c2h4 

2.16 

1 

74 

127 

44 

93 

3.12 

0.69 

48 

83 

39 

91 

3.88 

0.55 

38 

63 

34 

91 

4.92 

0.43 

30 

49 

27 

92 
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creating  a  fan  pattern.  When  V'0  is  small,  at  low  velocities  the  liftoff  heights  axe 
lower  than  predicted,  and  at  high  velocities  liftoff  heights  are  higher  than  predicted. 

One  possible  reason  for  the  disagreement  between  experimental  results  and  the 
model  is  that  as  Y0  is  decreased,  the  flames  axe  stabilized  closer  to  the  nozzle  ( h/ds  is 
reduced)  thereby  moving  out  of  the  fax  field  of  the  jet  where  self-similar  behavior  is 
most  completely  achieved.  Additionally,  there  may  be  some  complicated  Reynolds 
number  effects  in  this  range  of  Re^.  These  issues  are  discussed  later  in  this  section, 
following  the  presentation  of  all  of  the  experimental  data. 


FIG.  3.11  Normalized  liftoff  height  vs.  normalized  velocity  for  C2H6,  for  varying  fuel 
mass  fraction  Y0  in  the  nozzle  fluid. 

Results  similar  to  those  using  natural  gas  were  obtained  for  C2H6  as  shown 
in  Fig.  3.11.  Again,  as  Ya  is  reduced,  the  slope  of  the  line  increases.  The  range  of 
Re 00  for  each  of  the  curves  in  Fig.  3.11  is  ~  12,000  to  20,000,  and  the  maximum 
Mach  numbers  Ms  are  0.26,  0.20,  0.16,  and  0.12  for  Y0  =  1,  0.79,  0.63,  and  0.51, 
respectively.  This  is  still  in  the  range  where  Re oo  could  have  an  effect,  although 
these  values  of  Re ^  sure  roughly  twice  sis  high  as  in  the  natural  gas  flames  (Fig.  3.10). 
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FlG.  3.12  Normalized  liftoff  height  vs.  normalized  velocity  for  C2H4,  for  varying  fuel 
mass  fraction  Ya  in  the  nozzle  fluid. 

Experiments  with  air  premixing  were  also  performed  using  C2H4.  The  results 
are  plotted  in  Fig.  3.12  and  show  behavior  similar  to  the  previous  fuels.  The  slope 
of  the  line  increases  as  Ya  is  reduced.  For  each  of  the  curves  in  Fig.  3.12  the  range 
of  Re oo  is  ~  9,000  to  16,000,  falling  in  between  those  of  the  natural  gas  and  C2H6 
flames.  The  maximum  values  of  Ma  are  0.37,  0.24,  0.18,  and  0.14  for  Y0  =  1,  0.69, 
0.55,  and  0.43,  respectively. 


Figure  3.13  shows  more  clearly  how  the  slope  of  the  liftoff  curve  changes  with 


Y0.  Here,  the  slope  of  each  of  the  lines  of  Figs.  3.10  through  3.12  are  plotted  against 
1  /Y0.  Each  slope  was  found  by  fitting  a  parabola  through  the  data  and  taking  the 


derivative  at  the  point  where  ( Yst/Y0)h/ds  =  2.0,  *.e., 


d[(Yst/Y0)h/ds) 

d[{Ysi/Y0)2Us/ds] 


at  ( Ysi/Y0)h/ds  =  2.0  . 


(3.18) 


Each  slope  was  then  normalized  by  the  slope  for  Y0  =  1. 


As  increased  amounts  of  air  Eire  added  to  the  nozzle  fluid,  the  slope  of  the  curve 
increases,  in  disagreement  with  predictions  baised  on  far-field  scaling  of  turbulent 


1.0  1.4  1.8  2.2  2.6 

1/Y0 

Fig.  3.13  Variation  of  the  slope  of  the  curve  of  ( Yai/Y0)h/ds  vs.  ( YsX./Y0)2Us/ds  at 
(Ysi/Y0)h/ds  =  2.0  with  Y0.  The  slopes  have  been  normalized  by  the 
slope  for  Y0  =  1. 

jets,  as  was  used  in  the  strain-rate  model  described  in  Sec.  3.2.  Although  there 
is  some  scatter  in  this  plot,  there  is  a  reasonable  collapse  for  the  different  fuels. 
Such  a  collapse  was  not  obtained  in  the  results  of  Miake-Lye  and  Hammer  (1989), 
presumably  because  of  the  changing  Reynolds  number  as  Y0  was  changed. 

As  mentioned  above,  near-field  effects  and  Reynolds  number  effects  may  be 
responsible  for  the  discrepancy  between  experimental  results  and  predictions  based 
on  far-field  scaling.  In  particular,  the  mixing  virtual  origin  found  by  Gilbrech  (1991) 
in  his  study  of  flame  lengths  may  be  important.  The  normalized  mixing  virtual 
origin,  £,  (see  Eq.  3.16  and  the  discussion  in  Sec.  3.3.2)  was  found  to  have  a  strong 
dependence  on  Reynolds  number  in  the  range  10,000  <  Re^  <  150,000. 

The  non-dimensional  liftoff  height  can  be  defined  to  include  a  virtual  origin  x0 
(in  general,  different  from  the  flame  length  virtual  origin),  i.e., 
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Considering  first  the  results  for  natural  gas,  if  this  quantity  is  plotted  on  the  vertical 
axis  of  Fig.  3.10,  the  effect  would  be  to  shift  the  curves  down  for  Ya  <  1,  assuming 
xQ/ds  is  positive  and  independent  of  Re oo.  This  would  not  collapse  the  data  for 
different  Y0.  To  collapse  the  data,  xa/ds  would  need  to  be  ~  —4  at  Re0 0  =  6,000  and 
+6  at  Reoo  =  12,000.  This  can  be  compared  with  the  flame-length  results,  in  which 
B  was  found  to  increase  from  17  to  28  as  Re^  increases  from  10,000  to  20,000.  The 
required  change  in  x0/ds  with  Reoo  is  therefore  comparable  to  that  observed  in  the 
flame-length  studies  (although  the  range  of  Re  oo  is  different).  However,  a  negative 
virtual  origin  would  be  required  at  the  lower  values  of  Reoo,  and  this  is  in  contrast 
to  the  flame-length  results.  The  flame  lengths  measured  were  60  to  130  diameters 
from  the  nozzle  exit,  and  so  it  would  not  be  expected  that  the  extrapolated  origin 
B  would  be  the  same  as  xa/ds  of  the  current  measurements  which  were  made  in  the 
range  of  h/da  =  10  to  40. 

The  C2H6  results  plotted  in  Fig.  3.11  show  the  same  behavior  as  the  natural 
gas  results.  However,  the  values  of  Re  oo  are  roughly  twice  as  high  as  in  the  natural 
gas  flames  (Fig.  3.10).  As  a  result,  Re  oo  at  the  lower  liftoff  heights  of  Fig.  3.11  is 
approximately  equal  to  Re  oo  at  the  upper  liftoff  heights  of  Fig.  3.10.  The  required 
dependence  of  x0/ds  on  Re  oo  is  then  different  for  the  C2H6  and  natural  gas  results. 
This  rules  out  the  possibility  that  it  is  a  change  in  virtual  origin  with  Re  oo  that 
is  responsible  for  the  results  for  premixing  with  air.  Although  there  are  Reynolds 
number  effects,  they  are  not  the  dominant  reason  for  the  lack  of  a  collapse  at  the 
different  values  of  Y0. 

A  possible  explanation  for  the  systematic  variation  with  Y0  seen  in  Figs.  3.10 
to  3.13  can  be  obtained  by  modifying  the  strain-rate  model  to  allow  two  different 
virtual  origins,  one  for  the  concentration  field  and  one  for  the  strain-rate  field. 
Equation  3.3  for  the  concentration  field  is  modified  to  have  a  virtual  origin  iy, 


y  =  Cyn(j.^*3,y)A-(y/(j  -  ly))  , 


(3.19) 


38 


FlG.  3.14  Normalized  liftoff  height  vs.  normalized  velocity  using  Eqs.  3.19,  3.20,  and 
3.21  in  the  strain-rate  model.  Fst  =  0.06,  xy/ds  =  0.5,  and  x„/ds  =  1.75. 

and  Eq.  3.5  is  allowed  to  have  a  different  virtual  origin  x„, 

°L  =  C.  ,  /„(»/(*  -  x„»  .  (3.20) 

tx  Xa) 

To  ensure  a  linear  relationship  when  iy  =  xa  —  0,  it  is  assumed  that 

Ml)  =  Ml)  =  .  (3.21) 

Making  the  substitution  rio0  =  Cff/(Cy  aCT),  (Yst/Y0)h/ds  can  be  found  as  a  function 
of  (lst/Fo)27io017s/^s,  for  choices  of  the  virtual  origins,  xy  and  x9.  (With  this 
definition  of  7io0,  when  xy  —  xa  =  0,  ( Yst/Y0)h/ds  =  {Yst/Yo)2T\o0Us/ds.) 

It  is  assumed  that  xy/ds  and  xa/ds  are  not  changed  by  varying  Y0.  For  choices 
of  xy/ds  =  0.5  and  xa/d9  =  1.75,  and  with  Kst  =  0.06,  the  results  of  the  calculation 
for  Y0  =  1  and  Y0  =  0.5  are  shown  in  Fig.  3.14.  Here,  the  same  qualitative  behavior 
is  obtained  as  in  the  data.  As  Ya  is  reduced,  the  slope  of  the  curve  increases. 
However,  the  magnitude  of  the  change  is  less  than  observed.  The  ratio  of  the  slope 
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of  the  Y0  =  0.5  curve  to  the  Y0  =  1.0  curve  at  ( Yst/Y0)h/ds  =  2.0  is  only  1.14,  less 
than  the  experimental  value  of  1.5  obtained  from  Fig.  3.13. 

There  is  not  much  room  for  adjusting  the  values  of  the  virtual  origins.  xa  must 
be  greater  than  iy  to  have  the  slope  increase  as  Y0  is  decreased.  Increasing  xy/ds 
and  Xtr/ds  by  the  same  amount  has  the  effect  of  moving  the  Ya  =  0.5  curve  up  more 
than  the  Y0  —  1  curve,  without  changing  the  slope.  If  larger  values  of  {xa  —  xy  )/ds 
are  tried,  no  solution  can  be  found  for  the  lower  values  of  the  non-dimensional 
velocity,  because  the  ai  =  aCT  curve  is  moved  entirely  inside  the  Y  =  l'st  curve, 
t.e.,  the  flames  should  begin  at  the  nozzle  exit. 

Although  the  details  of  the  actual  physical  process  of  the  flame  stabilization 
certainly  deviate  from  the  above  model,  this  does  provide  a  possible  qualitative 
explanation  for  the  observed  variation  of  liftoff  height  with  Ya  and  Us. 
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CHAPTER  4 

Fluctuations  of  the  Liftoff  Height 


Measurements  of  liftoff  height  vs.  time  can  reveal  information  concerning  the 
dynamics  of  lifted  flames.  In  this  chapter  aspects  of  the  temporal  fluctuations 
in  liftoff  height  axe  examined.  The  motivation  for  making  these  measurements 
was  to  gain  insight  into  the  stabilization  process.  Particular  issues  of  interest  are 
whether  the  magnitude  of  the  fluctuations  scales  with  the  local  large  scale  of  the  jet, 
whether  the  time  scale  of  the  fluctuations  scales  with  the  local  passage  time  of  large 
structures,  and  whether  the  liftoff  height  fluctuates  primarily  at  one  frequency.  It 
was  hoped  that  answering  these  questions  would  help  elucidate  the  importance  of 
large  turbulent  structures  to  the  flame  stabilization  mechanism. 

To  enable  accurate  analyses  of  the  data,  the  measurements  reported  in  this 
chapter  were  made  for  longer  times  than  those  of  Ch.  3.  For  most  of  the  conditions, 
131072  (217)  data  points  were  acquired.  Exceptions  to  this  were  the  measurements 
of  Sec.  4.4  where  h  and  the  flame  length  were  measured  simultaneously  and  those 
of  Sec.  4.5,  where  h  and  the  flame  radiation  were  measured  simultaneously,  in  both 
of  which  65536  measurements  of  h  were  acquired. 

One  difference  in  the  experimented  procedure  for  the  data  of  this  chapter  from 
that  of  Ch.  3  is  that  the  array  was  repositioned  for  each  flow  condition.  The  vertical 
height  of  the  array  was  adjusted  to  be  at  the  mean  liftoff  height,  and  the  horizontal 
distance  from  the  flame  to  the  array  was  chosen  such  that  the  flame  image  just 
filled  the  lateral  extent  of  the  linear  array,  as  illustrated  in  Fig.  3.2.  This  was  done 
to  prevent  any  systematic  error  related  to  either,  the  fraction  of  the  flame  width 
imaged,  or,  the  relative  spatial  resolution,  from  contaminating  the  results;  both  of 
these  factors  were  kept  constant  for  all  the  conditions  reported  in  this  chapter. 
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Setting  a  sufficiently  low  threshold  for  determining  h  would  allow  the  lowest 
point  anywhere  along  the  flame  base  to  be  recorded  as  the  liftoff  height.  To  distin¬ 
guish  the  flame  light  from  noise  in  the  data,  a  sufficiently  large  threshold  had  to  be 
set;  the  technique  therefore  required  that  a  large  enough  amount  of  light  be  present 
to  call  that  position  ft.  To  a  reasonable  approximation,  however,  the  liftoff  heights 
measured  in  these  experiments  represent  the  lowest  position  along  the  lateral  extent 
of  the  flame  base  where  burning  exists. 


4.1  Fluctuation  Amplitudes  of  the  Liftoff  Height 

Before  discussing  the  temporal  behavior  of  /i(f),  results  on  the  spatial  scales 
of  the  fluctuations  are  presented.  This  includes  measures  such  as  the  probability 
density  function  (pdf)  and  root-mean-squaxe  (rms)  of  the  liftoff  height. 


h  /  h 


FlG.  4.1  Probability  density  function  (pdf)  of  liftoff  height.  C2H6,  ds  =  3.88  mm, 
Us  =  88  m/s,  h  =  151  mm. 
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Table  4.1:  Flow  Conditions  Studied,  Group  A  Nozzles 


fuel 

d 

(mm) 

us 

(m/s) 

h 

(mm) 

Reh 

rl/2 

(ms) 

n.g. 

4.61 

36 

88 

8,500 

58.6 

47 

116 

10,900 

63.7 

57 

144 

13,000 

76.9 

65 

172 

14,600 

110.1 

c2h6 

3.88 

56 

98 

15,700 

39.6 

74 

124 

20,000 

41.8 

88 

151 

23,500 

51.8 

96 

166 

25,300 

66.9 

c2h4 

2.16 

83 

49 

12,400 

21.0 

98 

59 

14,300 

20.5 

116 

74 

16,900 

22.7 

127 

83 

18,200 

25.5 

c2h2 

1.02 

173 

20 

11,800 

10.9 

193 

24 

12,900 

10.1 

219 

30 

14,300 

10.1 

235 

34 

15,200 

10.6 

The  histogram  of  h  from  one  of  the  flow  conditions  is  shown  in  Fig.  4.1,  where 
the  rectangular  bins  each  have  a  width  corresponding  to  two  pixels  on  the  linear 
array.  In  the  limit  of  infinitesimal  bin  size,  the  histogram  would  converge  to  the 
true  pdf.  The  curve  has  been  normalized  such  that  its  area  is  one,  t.e., 

OO 

pd{(h/h)d(h/h)  =  1  .  (4.1) 

The  pdf  is  mono-modal,  centered  approximately  at  the  mean  and  reasonably  sym¬ 
metric.  The  maximum  excursions  of  h  from  the  mean  are  ~  ±0.25  h,  making  the 
range  of  the  fluctuations  just  larger  them  the  local  jet  diameter,  6  =  0.44  h.  For 
these  data  the  rms  level  h'  is  0.066  h,  which  is  slightly  less  than  our  earlier  measure¬ 
ments  (Miake-Lye  and  Hammer  1989)  in  which  h' /h  was  found  to  be  ~  0.08  for  a 
wide  range  of  liftoff  heights.  That  the  spatial  scale  of  the  fluctuations  is  compara¬ 
ble  to  the  jet  large  scale  is  consistent  with  the  idea  of  large-scale  structures  of  the 
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FlG.  4.2  Normalized  rms  level  h' /h  vs.  normalized  liftoff  height  for  the  conditions 
listed  in  Table  4.1  (group  A  nozzles). 

turbulence  governing  the  dynamics  of  the  liftoff  height. 

For  each  of  the  first  four  fuels  listed  in  Table  3.3,  experiments  were  conducted 
with  two  nozzles.  This  allowed  two  different  Reynolds  numbers  to  be  investigated 
at  each  non-dimensional  liftoff  height.  The  conditions  achieved  with  the  smaller  set 
of  nozzles,  denoted  group  A,  a r?  listed  in  Table  4.1  and  those  of  the  larger  group 
of  nozzles,  group  B,  are  listed  in  Table  4.2.  The  Reynolds  number  listed  Re*  is 
the  local  Reynolds  number  just  below  the  flame  base,  which  is  approximately  equal 
to  the  far-field  Reynolds  number  Re =  Usds/u^.  Its  precise  definition  is  given 
in  Sec.  A. 5.  For  each  fuel,  four  different  liftoff  heights  were  studied,  spanning  the 
range  of  h  from  ~  1.5  to  2.5,  where  h  is  the  non-dimensional  liftoff  height,  which  is 
approximately  equal  to  Ysih/ds.  Its  precise  definition,  which  accounts  for  the  fact 
that  the  density  just  below  the  liftoff  height  is  not  equal  to  p^,  is  given  in  Secs.  4.3 
and  A.4.  Having  h  >  1.5  ensured  that  all  flames  were  lifted  to  at  least  20 <4-  For 
the  three  pure  fuels,  ds  ~  d,  and  so  all  flames  were  lifted  to  at  least  20  nozzle 
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diameters.  For  natural  gas,  all  flames  were  lifted  to  at  least  16  nozzle  diameters. 

A  quantity  derivable  from  the  pdf  is  the  rms  of  the  fluctuations,  h' .  By  making 
accurate  measurements  of  h  over  a  wide  range  of  conditions,  the  current  experiments 
are  able  to  provide  more  information  on  the  fluctuation  levels  of  lifted  flames  than 
previous  investigations.  In  Fig.  4.2  the  normalized  rms  level  h' /h  is  plotted  against 
normalized  liftoff  height  for  all  the  conditions  listed  in  Table  4.1.  These  rms  values 
have  been  corrected  for  the  finite  sampling  time  of  the  linear  array  as  described  in 
Appendix  C. 


Table  4.2:  Flow  Conditions  Studied,  Group  B  Nozzles 


fuel 

d 

(mm) 

Us 

(m/s) 

h 

(mm) 

Reh 

rl/2 

(ms) 

n.g. 

7.73 

46 

148 

17,900 

60.9 

61 

183 

23,800 

62.0 

79 

221 

30,300 

69.7 

90 

245 

34,300 

81.9 

c2h6 

4.92 

68 

125 

24,000 

38.7 

90 

157 

31,000 

43.5 

109 

187 

36,700 

53.9 

120 

214 

39,800 

80.4 

C2H4 

3.12 

105 

71 

22,300 

21.0 

128 

87 

26,700 

22.7 

154 

106 

31,600 

25.7 

171 

121 

34,700 

31.0 

C2H2 

1.55 

215 

30 

22,100 

8.2 

254 

38 

25,500 

8.7 

The  most  apparent  observation  from  Fig.  4.2  is  that  as  h  increases,  the  normal¬ 
ized  fluctuation  level  increases.  Because  the  Reynolds  number  also  increases  as  h 
is  increased,  however,  it  can  not  be  discerned  simply  from  these  data  whether  this 
increase  in  h' /h  is  caused  by  the  change  in  liftoff  height  or  by  the  accompanying 
change  in  Reynolds  number.  This  issue  is  addressed  below. 
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Fig.  4.3  Normalized  rms  level  h' /h  vs.  normalized  liftoff  height  for  the  conditions 
listed  in  Table  4.2  (group  B  nozzles). 

The  data  from  the  different  fuels  do  not  collapse  onto  a  single  curve.  In  general, 
the  rms  levels  for  natural  gas  and  C2H6  are  higher  than  those  of  C2H4  and  C2H2. 
It  is  not  obvious  why  this  should  be  the  case.  The  time  scale  of  the  fluctuations 
decreases  in  the  same  order  as  the  (first  four)  fuels  listed  in  Table  3.3.  The  possibility 
that  the  finite  sampling  time  of  the  measurements  would  cause  lower  measured  rms 
levels  for  the  faster  fuels  was  considered,  and  ruled  out,  as  discussed  in  Appendix  C. 

To  investigate  the  effect  of  Reynolds  number  on  the  flames,  larger  nozzles  of 
group  B  listed  in  Table  4.2  were  used.  The  same  range  of  h  was  investigated,  except 
that  for  acetylene,  the  maximum  attainable  value  was  1.8  because  of  limitations  on 
the  delivery  pressure  of  that  gas.  These  Reynolds  numbers  are  roughly  twice  as 
large  as  those  of  the  group  A  nozzles. 

The  relative  fluctuation  levels  h' /h  are  not  significantly  different  from  those  of 
the  smaller  nozzles  of  group  A,  sis  shown  in  Fig.  4.3.  The  same  general  behavior  is 
evident;  as  the  liftoff  height  increases,  the  relative  fluctuation  level  increases. 
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FlG.  4.4  Normalized  rms  level  b!  fh  vs.  normalized  liftoff  height.  Data  from  group 
A  and  B  nozzles  interpolated  to  Re/,  =  22,000. 

To  isolate  the  effects  of  h  from  those  of  Re*,  interpolation  was  performed  on 
the  data  from  the  group  A  and  B  nozzles  to  achieve  a  constant  Reynolds  number  of 
22,000.  Before  doing  this  interpolation,  the  data  within  each  group  of  nozzles  were 
interpolated  to  four  values  of  h,  namely  1.5,  1.8,  2.2  and  2.45.  Then,  at  each  value 
of  h,  the  interpolation  in  Reynolds  number  was  performed.  The  resulting  plot  is 
shown  in  Fig.  4.4,  where  the  horizontal  axis  has  been  changed  to  include  only  the 
range  of  h  studied.  Clearly,  there  is  an  increase  in  h'/h  with  h,  even  at  a  constant 
Re/,;  as  flames  approach  blowout,  the  relative  fluctuation  level  increases. 


4.2  Temporal  Behavior  Common  to  All  Flames  Studied 

Several  features  of  the  liftoff  height  vs.  time  signals  were  found  to  be  present 
in  all  the  conditions  studied,  and  these  are  illustrated  by  looking  at  results  for 
a  particular  flow  condition.  The  ways  in  which  the  behavior  changes  with  flow 
conditions  are  examined  in  later  sections. 
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The  measurements  of  h  vs.  t  revealed  that  fluctuations  in  h  occur  on  much 
longer  time  scales  than  the  large-scale  time  Tf  =  S(h)/Uc\(h).  This  was  not  expected 
as  it  conflicted  with  ideas  on  the  dynamics  and  structure  of  turbulent  jets.  A  premise 
of  the  strain-rate  model  (Sec.  3.2)  is  that  the  passage  of  large,  turbulent  structures 
controls  the  liftoff  height  fluctuations,  from  which  it  would  be  expected  that  the 
fluctuations  in  liftoff  height  would  occur  on  a  time  scale  approximately  equal  to  Tf . 
A  sample  time  trace  of  h  illustrating  this  long-time  behavior  is  presented  in  Fig.  4.5. 
(Only  1/256  of  the  entire  data  record  is  plotted.)  The  experimental  conditions  are 
the  same  as  those  of  Fig.  4.1.  For  reference,  10  Tf  =  46.7  ms  is  marked  on  the  figure. 
The  observed  fluctuations  clearly  occur  on  times  much  greater  than  Tf.  Although 
small  fluctuations  occur  on  the  jet  time  scale,  larger  fluctuations  occur  on  longer 
time  scales.  Another  observation  is  that  no  periodic  behavior  is  present  in  the  h  vs. 
t  signal. 


FlG.  4.5  Liftoff  height  h  vs.  time.  C2H6,  d  =  3.88  mm,  Us  =  88  m/s,  h  =  151  mm. 

To  quantify  the  time  scale  and  frequency  content  of  the  fluctuations,  the  power 
spectrum  of  the  h  vs.  t  signal  was  computed  with  program  PSD,  developed  by  Prof. 
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Paul  Dimotakis  Then,  all  the  computed  points  within  each  third  of  an  octave 
were  averaged  to  produce  a  smooth  curve  for  plotting.  The  spectral  power  £\(/) 
is  normalized  such  that 

h'2  =  2  f°£n(/) if.  (4.2) 

Jo 

The  spectrum  from  the  data  used  in  Fig.  4.5  is  plotted  in  Fig.  4.6,  where  the  fre¬ 
quency  /  has  been  normalized  using  r$.  At  low  frequencies  the  spectrum  is  fla„;  at 
a  break  frequency  /*,  the  spectrum  begins  its  descent  at  approximately  a  constant 
slope,  somewhat  steeper  than  —2.  There  is  no  peak  in  the  spectrum  that  would 
indicate  periodic  behavior.  Corresponding  to  the  long  time  behavior  noted  above, 
the  break  frequency,  obtained  from  straight  line  fits  to  the  two  linear  portions  of 
the  spectrum,  is  much  less  than  I/7-5,  approximately  0.01/r$. 


1°9 i o ( f  t6) 

FlG.  4.6  Power  spectrum  of  h.  Same  conditions  as  Fig.  4.5. 

As  a  complementary  measure  of  the  behavior  of  h(t),  the  autocorrelation  func¬ 
tion  was  computed  using  a  FFT  algorithm  (Press  et  al.  1988).  This  is  plotted  in 
Fig.  4.7,  where 


Rhh{r)  =  (h(t)  -  h)(h(t  +  t)~  h) 


(4.3) 
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FlG.  4.7  Autocorrelation  of  h.  Same  conditions  as  Fig.  4.5 

has  been  normalized  by  the  mean  square  level  to  force  the  curve  to  go  to  one  at  zero 
lag  time  r.  Again,  10  has  been  marked  on  the  figure,  illustrating  that  fluctuation 
times  are  much  larger  than  .  The  autocorrelation  shows  no  significant  undershoot 
of  zero  or  second  peak,  indicating  that  periodic  behavior  of  h(t)  is  not  present. 
Also  marked  on  the  figure  is  Tj /2,  which  is  defined  as  the  time  lag  r  for  which  the 
autocorrelation  function  falls  to  a  value  of  0.5.  (Values  of  Ty/ 2  have  been  corrected 
for  the  effect  of  non-instantaneous  sampling  time  as  described  in  Appendix  C.)  This 
time  is  used  throughout  this  thesis  as  a  representative  time  of  the  fluctuations  for 
a  given  condition.  Later  in  this  chapter  the  measured  values  of  Ty/ 2  for  different 
experimental  conditions  are  compared. 
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4.3  Variation  of  the  Fluctuation  Time  with  Conditions 

This  section  examines  the  dependence  of  the  fluctuation  time  on  the  flow  condi¬ 
tions.  The  measurement  of  h  vs.  t  reported  in  Sec.  4.2  indicates  that  the  correlation 
time  is  much  larger  than  the  local  jet  time  r^,  at  least  for  one  experimental  con¬ 
dition.  Whether  this  behavior  is  found  in  all  flow  conditions  is  the  topic  of  this 
section.  Of  particular  interest  is  whether  Ti/2  is  proportional  to  7-5  for  all  flow  con¬ 
ditions.  By  changing  the  fuel  type,  nozzle  diameter,  and  liftoff  height  (by  changing 
the  exit  velocity),  ts  can  be  varied  and  the  resulting  correlation  time  measured. 
As  the  experimental  conditions  are  changed,  several  non-dimensional  parameters 
change  which  could  potentially  alter  the  fluctuation  time. 

First  is  the  non-dimensional  mean  liftoff  height  h,  defined  as  Cy  (Ist/^cl),  35 
described  in  Sec.  A. 4.  For  most  of  the  conditions  to  be  described, 

h  *  Yath/ds  ,  (4.4) 

which  is  the  normalization  of  h  used  in  Ch.  3.  With  the  above  definition,  there  is  a 
one-to-one  correspondence  between  h  and  r]st,  the  value  of  77  at  which  Y  =  1st-  At 
low  liftoff  heights,  rjst  is  large  and  burning  takes  place  close  to  the  outside  of  the  jet. 
The  dynamics  of  such  a  flame  can  be  quite  different  from  one  with  a  large  liftoft 
height  close  to  blowout,  in  which  r/st  is  small.  In  the  experiments  to  be  discussed 
in  this  section,  several  fuels  were  used,  having  different  values  of  1st-  To  the  extent 
that  reaction  rates  of  these  fuels  axe  highest  at  Y  =  l^t  and  fall  off  rapidly  for 
higher  or  lower  values  of  Y,  h  is  the  correct  non-dimensionalization  of  the  liftoff 
height.  However,  some  of  the  fuels  have  wider  flammability  limits  than  others,  and 
somewhat  different  dependencies  of  reaction  rate  on  Y/Y’al.  Therefore,  h  is  not  a 
perfect  parameter  to  use  as  a  non-dimensional  liftoff  height  when  different  fuels  are 
compared,  but  it  is  a  reasonable  choice. 

Another  parameter  of  interest  is  the  Reynolds  number.  The  Reynolds  number 
in  the  non-reacting  part  of  the  jet,  just  below  the  liftoff  height,  Re/,,  is  estimated 
as  described  in  Sec.  A. 5.  For  most  of  the  conditions  studied.  Re/,  is  nearly  equal 
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to  the  far-field  Reynolds  number  Re oc  =  t/s<4/t'oo-  The  exception  to  this  is  for  the 
86%  H2  /  14%  C2H6  flames  (to  be  discussed  in  Sec.  4.3.5),  for  which  Re/,  is  ~  10% 
to  20%  lower  than  Re ^  for  the  range  of  h  studied. 

In  the  burning  region  of  the  flame,  the  effective  Reynolds  number  differs  sig¬ 
nificantly  from  Re*  because  of  the  temperature  rise  which  lowers  the  density  and 
increases  the  viscosity.  An  approximate  measure  of  this  effect  can  be  obtained 
from  a  simple  one-dimensional  analysis  of  the  flame.  The  momentum  flux  J  will 
be  unchanged  in  the  initial  burning  region.  J  ~  pU^S2,  and  the  Reynolds  number 
Re  ~  (pUc\6)/p.  Therefore  Re  ~  J1/2 p1/2 / p.  Noting  that  p  ~  T-1,  and  talcing  as 
an  approximation  p  ~  T  s  results  in  Re  ~  T~13.  An  appropriate  average  value  of 
T  can  be  assumed  to  be  T0 0  +  0.6A Ty,  (see  Sec.  A. 6,  which  utilizes  measurements  of 
Becker  and  Yamazaki  (1978))  where  A Tf  is  the  adiabatic  flame  temperature  rise  of 
the  fuel,  ~  1,900  K  for  most  hydrocarbons.  This  implies  that  the  effective  Reynolds 
number  in  the  flame  will  be  ~  0.13  Rex.  As  many  approximations  were  made  in 
this  analysis,  including  the  1-D  assumption,  the  factor  of  0.13  is  only  approximate. 

Additionally,  buoyancy  may  affect  the  dynamics  of  lifted  flames.  The  relative 
importance  of  buoyancy  to  momentum  is  characterized  by  a  Richardson  number 
relevant  to  the  dynamics  of  the  flame  base  ( Rih),  which  is  the  ratio  of  a  buoyancy 
force  to  momentum  flux.  The  flow  entering  the  flame  base  has  a  momentum  flux 
approximately  equal  to  the  source  momentum  flux,  Ja  =  (n/4)d^U2 .  There  is  some 
subjectivity  in  determining  the  proper  region  over  which  to  measure  the  buoyancy 
force  B.  The  flow  field  dynamics  are  influenced  by  the  region  of  the  flame  some 
distance  downstream  of  the  flame  base.  It  was  decided  to  use  the  region  from  the 
liftoff  height  to  one  local  jet  width  6(h)  downstream. 

The  details  of  the  definition  and  calculation  method  of  Rih  are  given  in  Sec.  A. 6, 
and  the  resulting  definition  of  =  Rih'  is 


(4.5) 
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where  is  a  constant  equal  to  0.87  that  accounts  for  the  fact  that  it  is  the  frustrated 
cone  above  h  that  provides  the  buoyancy  force,  rather  than  the  entire  jet  cone.  A p 
is  the  average  value  of  p ^  —  p  in  the  region  extending  from  h  to  h  +  8.  has 
a  length  scale  h  to  the  first  power.  It  is  similar  to  the  definition  of  Becker  and 
Yamazaki  (1978)  and  Becker  and  Liang  (1978),  but  with  h  as  the  relevant  length 
scale,  rather  than  the  flame  length  L ,  and  with  consideration  given  to  the  actual 
density  in  the  flame. 

In  addition  to  the  parameters  of  non-dimensional  liftoff  height,  Reynolds  num¬ 
ber,  and  Richardson  number,  other,  more  subtle  parameters  that  could  affect  the 
dynamics  of  the  liftoff  height  result  from  the  different  chemistry  of  the  fuels.  The 
flammability  limits  are  wider  for  some  of  the  fuels  than  the  others.  Also,  the  flame 
temperature  rise  is  not  exactly  the  same  for  the  fuels.  Although  this  variation  of 
AT/  has  little  effect  on  buoyancy,  and  has  been  accounted  for  in  the  definition  of 
£/,,  the  resulting  dilatation  of  the  gas  could  conceivably  change  the  nature  of  the 
flow  at  the  flame  base.  Also,  the  radiative  properties  of  the  fuels  differ  markedly, 
depending  largely  on  the  amount  of  soot  which  is  produced  in  the  flame. 


4.3.1  Effects  of  Liftoff  Height  and  Fuel  Type 

The  fuels  used  in  the  current  experiments  are  listed  in  Table  3.3,  along  with 
their  properties  relevant  to  these  experiments.  They  are  listed  in  order  of  increas¬ 
ing  reaction  rate  or  decreasing  effective  chemical  time  T|0  =  h/YstUs.  The  first 
experiments  to  be  discussed  here  involved  the  first  four  fuels  listed.  Hydrogen  was 
used  only  as  a  mixture  with  the  other  fuels  in  the  table,  and  hence  discussion  of  the 
experiments  involving  hydrogen  is  postponed  until  Sec.  4.3.5. 

This  section  addresses  the  effects  that  stabilization  distance  and  fuel  type  have 
on  the  temporal  behavior  of  the  flames.  The  effects  of  Reynolds  number  are  ad¬ 
dressed  in  Sec.  4.3.2,  followed  by  a  discussion  of  the  effects  of  buoyancy,  as  char¬ 
acterized  by  the  Richardson  number,  in  Secs.  4.3.3  and  4.3.4.  As  will  be  shown, 
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there  is  a  close  correspondence  between  the  fuel  type  and  the  resulting  Richardson 
number  at  the  flame  base. 

From  the  measured  values  of  rXj2  at  the  conditions  listed  in  Table  4.1  (group 
A  nozzles)  and  Table  4.2  (group  B  nozzles),  an  interpolation  procedure  like  that 
used  to  prepare  Fig.  4.4  in  Sec.  4.1  was  used.  For  each  fuel  and  each  nozzle,  values 
of  were  interpolated  to  four  values  of  h ,  1.5,  1.8,  2.2,  and  2.45.  To  produce 

a  graph  in  which  the  Reynolds  number  is  fixed  ( Reh  =  22,000),  interpolation  from 
the  group  A  and  B  nozzles  was  then  done  at  each  value  of  h. 

The  values  of  Tj j2  normalized  by  Tj  axe  plotted  in  Fig.  4.8  for  each  of  the  four 
fuels.  Because  the  maximum  value  of  h  achieved  for  C2H2  with  its  group  B  nozzle 
was  1.8,  only  the  two  points  corresponding  to  h  of  1.5  and  1.8  are  plotted  for  that 
fuel. 
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FlG. 4.8  Variation  of  normalized  correlation  time  Ti/2/t6  with  normalized  liftoff 
height  h,  at  a  fixed  Reynolds  number  of  Re*  =  22,000. 

All  of  the  fuels  have  correlation  times  much  larger  than  the  local  large-scale 
time,  with  normalized  correlation  times  Ti/2/t$  ranging  from  ~  10  to  40.  The 
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normalized  correlation  times  increase  with  fuel  type  in  the  same  order  as  Table  3.3, 
which  corresponds  to  decreasing  rj0.  The  effect  is  most  dramatic  for  C2H2,  which 
has  values  of  t1/2/ts  more  than  twice  those  of  the  other  fuels.  Even  C2H4,  however, 
begins  to  show  an  increase  in  tx /2/t6  over  the  values  of  natural  gas  and  C2H6.  The 
possibility  that  the  changing  Richardson  number  at  the  flame  base  with  fuel  type  is 
the  cause  of  the  variation  of  correlation  times  is  examined  in  Secs.  4.3.3  and  4.3.4. 

The  other  feature  of  the  fluctuations  evident  from  Fig.  4.8  is  the  variation  of 
correlation  time  with  h.  For  natural  gas,  there  is  very  little  change  in  tx with 
h.  C2H6  shows  a  concave  behavior,  where  T! /2/r#  first  decreases  then  increases  as 
blowout  is  approached.  For  C2H4  and  C2H2,  there  is  a  decrease  in  r, with 
liftoff  height,  although  only  the  lowest  two  values  of  h  are  available  for  C2H2. 


4.3.2  Effect  of  Reynolds  Number 

Figure  Fig.  4.9  shows  the  effect  of  Reynolds  number  on  the  flames.  Here,  tx /2/r6 
is  plotted  against  Re*  in  four  separate  plots,  with  one  plot  for  each  of  the  values 
of  h,  1.5,  1.8,  2.2,  and  2.45.  In  each  of  the  four  plots,  results  from  each  of  the  four 
fuels  are  presented. 

In  addition  to  the  group  A  and  B  conditions,  some  additional  measurements 
were  made  at  still  higher  Reynolds  numbers,  corresponding  to  the  conditions  listed 
in  Table  4.3;  these  points  are  included  in  Fig.  4.9.  For  natural  gas  and  C2H6,  data 
were  taken  only  with  the  group  A  and  B  nozzles.  For  C2H4  an  additional  nozzle 
(group  C)  was  used  for  h  =  1.5,  1.8,  and  2.2.  For  C2H2,  all  three  nozzles  were  used 
for  h  =  1.5,  the  group  A  and  B  nozzles  for  h  =  1.8,  and  only  the  group  A  nozzle 
for  h  =  2.2  and  2.45. 

In  general,  there  is  very  little  variation  of  normalized  correlation  time  with 
Reynolds  number  for  all  the  fuels  except  C2H2.  For  C2H2,  at  least  at  h  <  1.8,  there 
is  a  decrease  in  tx  j2  / Tf,  with  Re*,  for  the  range  of  Re&  shown.  It  is  not  obvious 
why  one  fuel  should  exhibit  a  dependence  on  Reynolds  number  while  the  others 
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FlG.  4.9  Variation  of  Ti/2/t6  with  Re/,  and  fuel  type  at  four  values  of  h  =  1.5,  1.8, 
2.2,  and  2.45.  See  Fig.  4.8  for  definitions  of  symbols.  Note  the  different 
vertical  axes  of  the  top  and  bottom  plots. 


Table  4.3:  Flow  Conditions  Studied,  Group  C  Nozzles 


fuel 

d 

(mm) 

us 

(m/s) 

h 

(mm) 

Reh 

rl/2 

(ms) 

c2h4 

3.88 

115 

87 

30,400 

21.2 

148 

108 

38,300 

23.9 

179 

130 

45,800 

27.7 

C2H2 

2.16 

186 

32 

28,000 

8.3 

245 

43 

35,400 

8.4 

do  not.  One  possibility  is  that  C2H2,  which  has  a  higher  flame  temperature  than 
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the  other  fuels  (Table  3.3),  has  a  lower  effective  Reynolds  number  at  a  given  value 
of  the  cold- jet  Reynolds  number  Re*.  Then,  while  the  range  of  flame  Reynolds 
numbers  is  high  enough  for  Reynolds  number  independence  in  the  other  fuels,  it  is 
low  enough  to  still  have  a  dependence  for  C2H2.  It  would  be  somewhat  surprising, 
however,  if  the  Reynolds  number  in  the  burning  region  were  more  important  than 
that  upstream  of  the  flame  base  (Re/,)  in  determining  the  dynamics  of  the  flame 
base. 


4.3.3  Variation  of  Correlation  Time  with  Richardson  Number 


The  correlation  times  of  different  fuels  can  be  compared  by  calculating  the  value 
of  £/,,  and  plotting  the  corresponding  value  of  T] /2/t6-  Faster  fuels,  with  smaller 
effective  chemical  times,  i\0  —  h/YstUs,  will  have  lower  Richardson  numbers.  This 
can  be  seen  by  combining  the  definition  of  tiq  with  the  definition  of  the  Richardson 
number  (Eq.  4.5)  to  get 


(4.6) 


Here,  C\:3  is  a  constant,  and  g'  =  g  Ap/poo  is  approximately  constant  for  the  various 
fuels.  Then,  if  h  «  Ysth/ds  is  fixed,  and  the  Reynolds  number  is  fixed  by  keeping 
Usds  constant,  £/,  will  be  proportional  to  n'./2- 


Figure  4.10  shows  the  results  of  all  three  groups  of  nozzles,  A,  B,  and  C.  For 
each  fuel/nozzle  pair,  four  points  axe  plotted,  corresponding  to  the  four  values  of 
h.  Beside  each  group  of  points,  the  corresponding  fuel  is  printed.  For  the  group 
A  conditions  (open  symbols),  there  axe  data  points  at  each  of  the  four  values  of  h 
for  each  of  the  four  fuels.  For  the  group  B  conditions  (filled  symbols),  there  are 
four  data  points  for  each  fuel  except  C2H2  for  which  only  the  lower  two  values  of 
h  are  plotted.  The  four  points  from  the  group  C  nozzles  (one  for  C2H2  and  three 
for  C2H4)  axe  plotted  with  a  -I-  inside  the  symbol.  These  points  make  the  graph 
somewhat  cluttered  for  C2H4  where  the  points  lie  close  to  those  for  the  group  B 
nozzles. 
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Fig.  4.10  Variation  of  normalized  correlation  time  with  Richardson  number 

£/,,  for  the  conditions  listed  in  Tables  4.1,  4.2,  and  4.3.  Open  symbols: 
group  A,  filled  symbols:  group  B,  symbols  with  a  +  inside:  group  C. 

There  is  a  reasonably  good  collapse  of  correlation  time  with  Richardson  num¬ 
ber.  Because  the  values  of  £/,  do  not  overlap  for  the  different  fuels,  however,  the 
results  are  not  as  robust  as  they  would  otherwise  be.  Also,  the  acetylene  results 
for  the  two  different  nozzle  diameters  do  not  collapse  very  well.  It  could  be  argued, 
however,  that  the  imperfect  collapse  for  C2H2  is  due  to  slight  errors  on  the  estima¬ 
tion  of  £h  for  the  different  nozzles.  Alternatively,  this  may  be  a  Reynolds  number 
effect,  as  discussed  previously. 

The  variation  of  the  correlation  time  with  £*  is  quite  extreme  for  the  C2H2 
flames  while  for  the  other  fuels,  T\j2  approximately  scales  with  Tf,  (although  the 
upward  trend  is  apparent  for  C2H4).  It  is  surprising  that  at  such  low  Richardson 
numbers  as  the  C2H2  flames  have,  there  would  be  such  a  large  effect  of  further 
reduction  in  Although  extrapolating  the  curve  to  £*  — ►  0  appears  to  show 

*7"  1  /2 / T6  —*  00,  it  is  reasonable  to  expect  that  that  rx /2 / would  actually  go  to 
some  constant  value  as  £*  — ♦  0.  The  Richardson  numbers  of  all  the  flames  studied 
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are  relatively  low  by  conventional  standards.  Becker  and  Yamazaki  (1978)  found 
that  for  £x  <  1,  the  mass  flux  and  momentum  flux  are  relatively  independent  of 
£r  in  hydrocarbon  flames.  Because  they  used  the  approximation  A p  ~  p0 0  in  their 
definition  of  £,  a  more  appropriate  value  to  use  in  comparing  with  the  current 
definition  of  which  assumes  A p  ~  p^  —  p{Too  +  0.6  AT/)  would  be  ~  0.93. 
Similarly,  Gilbrech  (1991)  found  Richardson  number  independence  of  flame  length 
of  low-temperature  flames  for  <fr  <  1.  In  his  definition  of  £,  A p  «  p^-  piT^  +  ATf), 
and  so  a  more  appropriate  value  to  use  in  comparing  with  the  current  definition 
would  be  0.61/3  =  0.84.  Because  the  Richardson  numbers  in  the  current  experiments 
axe  below  that  in  which  momentum  dominance  is  believed  to  hold,  the  dependence  of 
ti/2/t(  on  £h  is  surprising.  However,  because  the  dynamics  of  the  flames  are  largely 
affected  by  the  outer  part  of  the  jet,  where  velocities  are  lower,  such  a  dependence 
on  £h  is  conceivable. 


Fig.  4.11  Variation  of  normalized  correlation  time  Ti/2/Ti  with  Richardson  number 
at  a  fixed  Reynolds  number  of  Re*  =  22,000. 


One  possible  reason  for  acetylene’s  longer  normalized  correlation  times  is  that 
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its  radiation,  which  is  considerably  greater  than  that  of  the  other  fuels  because  of 
greater  soot  production,  might  somehow  affect  the  dynamics  of  the  liftoff  height. 
This  possibility  is  discussed  in  Sec.  4.5. 

As  indicated  in  Fig.  4.9,  there  is  some  variation  of  ri/2  with  Reynolds  number, 
in  particular  for  C2H2.  In  Fig.  4.10,  however,  the  different  points  are  not  at  the  same 
Reynolds  number,  making  it  somewhat  difficult  to  distinguish  the  effects  of  from 
those  of  Reft.  Figure  Fig. 4.11  shows  the  data  after  interpolating  to  Re*  =  22,000, 
as  was  done  in  Fig.  4.8. 

The  data  fall  nearly  onto  a  single  curve.  There  is  a  slight  variation  with  h , 
however.  As  mentioned  above,  it  is  conceivable  that  the  non-dimensional  liftoff 
height  would  have  an  effect  on  the  flame  dynamics,  resulting  in  separate  curves  of 
ri/2 lTf>  vs-  th  f°r  the  different  values  of  h.  The  variation  of  Tj /2 / with  h  is  not 
very  large  for  the  fuels  other  than  acetylene,  however. 


4.3.4  Effect  of  Nozzle  Orientation 

To  help  discern  whether  buoyancy  is  responsible  for  the  variation  of  correlation 
time  with  fuel  type,  measurements  of  h  vs.  t  were  repeated  with  the  nozzles  oriented 
horizontally.  If  buoyancy  is  an  important  parameter,  then  having  the  buoyancy 
force  normal  to  the  flame  propagation  direction  would  be  expected  to  alter  the 
liftoff  height  dynamics,  although  it  would  not  have  the  same  effect  a s  eliminating 
buoyancy  altogether. 

The  flow  conditions  for  the  horizontal-flame  experiments  are  listed  in  Table  4.4. 
The  same  conditions  were  studied  as  for  the  (vertical)  nozzles  in  group  A,  except 
that  C2H6  was  not  used.  Also,  slightly  higher  velocities  were  needed  for  natural 
gas  and  C2H4  to  achieve  the  same  liftoff  heights,  as  shown  by  the  plot  of  h  vs.  Us 
(Fig.  4.12).  This  is  understandable  in  that,  in  the  horizontal  flames,  there  is  no  force 
pushing  the  hot  products  downstream  and  inhibiting  the  ignition  of  fresh  reactants 
as  there  is  with  the  vertical  flames.  This  variation  with  nozzle  orientation  would  be 
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Table  4.4:  Flow  Conditions  Studied,  Horizontal  Nozzles 


fuel 

d 

(mm) 

Us 

(m/s) 

h 

(mm) 

Reh 

t1/2 

(ms) 

n.g. 

4.61 

39 

92 

9,200 

52.7 

46 

108 

10,800 

58.2 

57 

133 

13,200 

74.7 

63 

149 

14,500 

81.1 

c2h4 

2.16 

86 

49 

12,600 

21.5 

103 

60 

14,900 

25.7 

124 

76 

17,600 

28.6 

136 

87 

19,100 

33.6 

c2h2 

1.02 

176 

20 

11,900 

10.9 

197 

25 

13,100 

10.8 

220 

30 

14,400 

10.7 

235 

34 

15,200 

11.5 

Us  (m/s) 


FIG.  4.12  Effect  of  nozzle  orientation  on  mean  liftoff  heights.  Open  symbols:  verti¬ 
cal  nozzles,  filled  symbols:  horizontal  nozzles. 
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most  extreme  for  the  more  buoyant  flames.  However,  even  for  natural  gas,  which 
has  the  largest  Richardson  numbers,  the  velocities  axe  only  ~  15%  larger  for  the 
horizontal  nozzles. 

The  correlation  times  are  plotted  in  Fig.  4.13  for  the  horizontal  and  vertical 
nozzles.  Within  the  resolution  of  the  measurements,  there  is  no  change  in  the 
dynamics  of  the  liftoff  height  with  nozzle  orientation.  This  is  true  for  all  three 
of  the  fuels  studied.  The  largest  effect  would  be  expected  for  natural  gas  flames, 
since  the  Richardson  number  is  the  highest  for  them,  but  neither  they  nor  the  other 
flames  appear  to  be  affected  by  the  nozzle  orientation. 


FlG.  4.13  Effect  of  nozzle  orientation  on  scaled  correlation  time  Ti/2/t{  vs.  Richard¬ 
son  number  £/,.  Open  symbols:  vertical  nozzles,  filled  symbols:  horizontal 
nozzles. 

The  result  that  the  direction  of  the  gravitational  vector  relative  to  the  nozzle 
orientation  has  no  measurable  affect  on  the  liftoff  height  dynamics  makes  the  col¬ 
lapse  of  Ti/2/ts  with  £/,,  a  measure  of  the  importance  of  buoyancy,  more  difficult 
to  explain.  It  is  still  possible  that  both  horizontal  and  vertical  nozzles  could  have 
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comparable  effects  of  buoyant  behavior  on  the  liftoff  height  dynamics.  Other  expla¬ 
nations  for  the  variation  of  r1/2/7'6  with  £/,  are  examined  in  the  following  sections. 

4.3.5  Results  for  Fuel  Mixtures 

To  gain  additional  control  of  the  flow  conditions,  experiments  wcr^  conducted 
using  mixtures  of  fuels.  By  appropriately  choosing  the  proportions  of  the  fuels, 
the  effective  chemical  time  of  the  fuel  mixture  could  be  controlled.  The  conditions 
studied  :  re  listed  in  Table  4.5.  A  mixture  of  50%  C2H2  /  50%  C2H6  (by  volume) 
was  used  primarily  to  compare  with  C2H4.  Mixtures  or  H2  with  each  of  C2H6, 
C2H4,  and  C2H2  were  used  to  obtain  effective  chemical  times  in  the  neighborhood 
of  that  of  C2H2. 

Only  one  nozzle  diameter  was  used  for  the  C2H2/C2H6  mixture,  the  same 
nozzle  as  was  used  for  C2H4  in  the  group  A  conditions.  For  the  H2/C2H6  mixture, 
three  different  nozzle  diameters  were  used  to  be  able  to  interpolate  to  Re*  =  22,000 
at  each  of  the  four  non-dimensional  liftoff  heights.  Two  nozzles  were  used  for  the 
H2/C2H4  experiments,  again  allowing  interpolation  to  Re/,  =  22,000.  Only  a  single 
nozzle  diameter  was  used  for  the  H2/C2H2  mixture. 

To  test  whether  the  collapse  of  rx /2  / T&  vs.  £/,  would  hold  for  fuels  other  than  the 
four  fuels  previously  reported,  the  mixture  of  50%  C2H2  /  50%  C2H6  was  used.  This 
resulted  in  the  same  average  molecular  weight  and  1st  as  C2H4.  The  measurements 
of  h  revealed  that  the  mixture’s  effective  chemical  time  is  slightly  lower  than  that 
of  C2H4,  yielding  slightly  lower  liftoff  heights  at  each  exit  velocity. 

Measurements  of  the  correlation  times  for  the  C2H2/C2H6  mixture  are  shown 
in  Fig.  4.14  along  with  the  group  A  results  and  the  H2/C2H2  results.  The  groups 
of  points  corresponding  to  each  fuel  (or  mixture)  are  labeled,  and  the  symbols 
alternate  between  filled  and  open  in  the  order  of  £/,.  The  50%  C2H2  /  50%  C2H6 
and  C2H4  results  nearly  lie  on  top  of  each  other,  although  the  former  have  slightly 
lower  values  of  £/,  at  each  value  of  h.  The  fact  that  the  50%  C2H2  /  50%  C2H6 
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Table  4.5:  Flow  Conditions  Studied,  Fuel  Mixtures 

fuel 

d 

(mm) 

Us 

(m/s) 

h 

(mm) 

Reh 

rl/2 

(ms) 

50%  C2H2/ 

50%  C2H6 

2.16 

89 

50 

13,300 

19.3 

104 

60 

15,300 

19.4 

122 

72 

17,600 

21.1 

136 

85 

19,500 

26.3 

86%  H2/ 

14%  C2H6 

1.55 

357 

42 

14,600 

10.9 

393 

46 

16,200 

10.8 

2.16 

264 

36 

15,200 

12.9 

320 

43 

18,100 

11.5 

375 

51 

21,300 

11.0 

442 

62 

25,300 

11.3 

3.12 

438 

59 

36,000 

10.3 

509 

68 

42,100 

50%  H2/ 

50%  C2H4 

2.16 

153 

44 

15,600 

172 

51 

17,400 

16.6 

■a 

63 

20,300 

16.9 

|  ..  | 

75 

22,600 

21.8 

3.12 

189 

59 

28,000 

15.5 

248 

76 

36,400 

17.2 

295 

92 

42,900 

19.7 

311 

102 

45,200 

21.4 

50%  H2/ 

50%  C2H2 

1.02 

290 

21 

13,300 

8.5 

265 

17 

12,300 

9.6 

346 

29 

15,600 

8.1 

points  fall  on  the  same  curve  as  the  pure  fuels  provides  support  that  £&  can  collapse 
results  of  different  fuels.  It  should  be  noted,  however,  that  this  is  not  a  particularly 
strong  test  of  the  collapse  with  because  the  mixture  of  50%  C2H2  /  50%  C2H6 
has  the  same  average  chemical  composition,  molecular  weight,  and  yst  as  C2H4. 
Discussion  of  the  H2/C2H2  mixture  results  (left-most  points)  is  postponed  until  all 
other  mixture  results  have  been  described. 


Otiier  fuel  mixtures  were  used  to  achieve  Richardson  numbers  near  to  those  of 
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FlG.  4.14  Variation  of  normalized  correlation  time  Tlf2/Ts  with  Richardson  number 
,  for  the  conditions  listed  in  Table  4.1  and  for  the  50%  C2H2  /  50%  C2H6 
and  50%  H2  /  50%  C2H2  mixtures  listed  in  Table  4.5.  Filled  symbols: 
n.g.,  C2H4,  and  C2H2. 

acetylene  flames.  This  was  done  largely  to  see  whether  the  longer  non-dimensional 
correlation  times  encountered  with  C2H2  were  peculiar  to  that  fuel.  Experiments 
previously  conducted,  and  described  in  Sec.  4.5,  showed  that  the  visible  radiation 
from  C2H2  flames  is  highly  correlated  with  fluctuations  in  h.  It  was  thought  that 
some  feedback  mechanism  from  the  flame  radiation  (and  associated  heat  loss)  to  the 
flame  base  might  be  the  cause  of  the  longer  correlation  times.  Mixing  H2  with  C2Hg 
or  C2H4  provided  a  good  way  to  achieve  flames  with  low  T|0  and  low  £*,,  without 
intense  radiation  from  soot.  Pure  hydrogen  could  not  be  used  because  its  visible 
radiation  is  insufficient  for  detection  with  the  photodiode  array. 

First,  a  mixture  of  86%  H2  /  14%.  C2Hg  was  used  The  effective  chemical  time 
of  the  mixture  was  found  to  be  slightly  larger  than  that  of  C2H2,  but  sufficiently 
close  for  the  present  purposes.  No  radiation  from  soot  was  observed  for  these  flames. 
Following  these  experiments,  studies  of  a  mixture  of  50%  H2  /  50%  C2H4  were  made. 
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Very  little  radiation  from  soot  was  observed  for  these  flames.  The  effective  chemical 
time  of  this  mixture  was  found  to  be  slightly  larger  than  the  86%  H2  /  14%  C2H6 
mixture.  By  using  more  than  one  nozzle  for  each  of  the  two  mixtures,  interpolation 
to  a  constant  Reynolds  number  of  22,000  was  made  possible.  Results  are  shown 
in  Fig.  4.15.  This  figure  is  similar  to  Fig.  4.11,  but  with  additional  data  points  for 
the  mixtures  of  H2  with  C2H6  and  C2H4.  Although  neither  of  the  fuel  mixtures 
achieved  quite  as  low  a  Richardson  number  as  pure  C2H2,  it  can  still  be  seen  that 
there  is  a  substantial  increase  in  /2 / as  £/,  becomes  small.  The  mixture  results 
fall  approximately  on  the  same  curve  as  the  pure-fuel  results,  although  the  50%  H2 
/  50%  C2H4  points  are  slightly  to  the  right  of  the  86%  H2  /  14%  C2H6  points. 


Fig.  4.15  Variation  of  normalized  correlation  time  T\j2lTb  with  Richardson  number 
£h,  comparing  C2H2  flames  with  mixtures  of  H2  with  C2H6  and  C2H4. 
Reh  =  22,000.  Filled  symbols:  n.g.,  C2H4,  86%  H2  /  14%  C2H6. 


The  basic  conclusion  from  Fig.  4.15  is  that  the  large  correlation  times  for  C2H2 
are  not  the  result  of  the  intense  radiation  from  such  flames.  The  H2/C2H6  and 
H2/C2H4  flames,  for  which  little  radiation  from  soot  is  present,  also  have  large 
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(normalized)  correlation  times. 

The  last  fuel  mixture  in  Table  4.5  is  50%  H2  /  50%  C2H2,  which  was  used 
to  achieve  Richardson  numbers  even  lower  than  for  pure  C2H2.  It  was  not  known 
whether  the  trend  toward  higher  Tj/2/t£  would  continue  as  £/,  was  decreased.  The 
results  are  plotted  in  Fig.  4.14,  some  aspects  of  which  has  already  been  discussed 
above.  The  left-most  point  at  each  of  the  four  liftoff  heights  is  from  the  H2/C2H2 
mixture.  It  is  evident  from  Fig.  4.14  that  Ti/2/ts  continues  to  increase  as  is 
decreased  beyond  the  range  achievable  with  C2H2.  If  the  curve  levels  out  at  low 
£h ,  it  must  occur  at  lower  values  than  those  yet  reached  in  these  experiments.  It  is 
also  worth  noting  that  the  Reynolds  numbers  of  the  C2H2  flames  and  the  50%  H2 
/  50%  C2H2  flames  are  comparable,  indicating  that  the  difference  in  r^/rg  is  not 
a  result  of  different  Reynolds  numbers. 

4.4  Simultaneous  Measurements  of  Liftoff  Height  and  Flame  Length 

Fluid  mechanical  time  scales  at  the  flame  tip,  where  the  length  scales  are  larger 
and  the  velocities  are  lower,  are  considerably  longer  than  the  time  scales  at  the  flame 
base  Because  the  measured  fluctuation  times  of  the  liftoff  height  were  found 

to  be  considerably  larger  than  r$(h),  there  was  some  speculation  that  the  behavior 
at  the  flame  tip  might  be  influencing  the  motion  of  the  flame  base.  To  investigate 
this  possibility,  measurements  of  both  the  liftoff  height  and  flame  length  were  made 
simultaneously. 

The  same  photodiode  array  was  used  for  both  measurements  by  positioning 
the  array  such  that  both  the  flame  base  and  the  flame  tip  were  in  view.  While  there 
is  a  sharp  transition  in  light  intensity  at  the  flame  base,  the  intensity  decays  much 
more  gradually  at  the  flame  tip.  This  results  in  some  ambiguity  in  determining  the 
flame  length  T,  but  simply  setting  a  threshold  in  intensity  proved  to  be  adequate. 

Two  flow  conditions  were  studied  using  natural  gas  with  the  4.61  mm  nozzle, 
and  one  condition  using  C2H6  with  the  3.88  mm  nozzle.  The  relevant  parameters 
are  listed  in  Table  4.6. 


Table  4.6:  Conditions  of  h  and  L  Measurements 

fuel 

d 

(mm) 

Us 

(m/s) 

h 

(mm) 

h' 

(mm) 

h 

L 

(mm) 

V 

(mm) 

n.g. 

4.61 

39 

111 

6.1 

1.80 

850 

66 

53 

153 

10.2 

2.48 

900 

64 

c2h6 

3.88 

80 

153 

9.2 

2.29 

1000 

72 

Portions  of  the  two  time  traces  for  the  higher  velocity  natural  gas  condition 
are  plotted  in  Fig.  4.16.  The  liftoff  height  values  have  been  multiplied  by  4  to 
allow  both  h  and  L  to  be  seen  clearly  on  one  plot.  In  the  flame  length  signal, 
ramp  structures  are  evident.  This  is  caused  by  burning  structures  propagating 
downstream  and  burning  out.  Once  the  structure  burns  out,  the  flame  length 
abruptly  drops  back  to  the  next-most  upstream  burning  region.  This  behavior  has 
been  observed  previously  in  liquid  flames  by  Dahm  and  Dimotakis  (1987)  and  in 
hydrocarbon  flames  by  Mungal  et  al.  (1991). 


FlG.  4.16  Liftoff  height  h  and  flame  length  L  vs.  time.  Natural  gas,  d  =  4.61  mm, 
Us  =  53  m/s. 
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Looking  at  the  time  traces,  there  is  evidently  very  little  correlation  between  the 
h  and  L  signals.  This  is  supported  by  Figs.  4.17  and  4.18,  which  show  correlations 
for  the  two  natural  gas  cases  listed  in  Table  4.6.  In  each  of  those  two  figures,  three 
curves  are  plotted.  The  autocorrelation  of  h  and  that  of  L  axe  shown,  along  with 
the  cross  correlation  of  h  and  L,  where  a  positive  time  lag  is  defined  for  L  lagging  h. 
The  autocorrelation  curves  have  again  been  normalized  by  the  mean  square  value, 
forcing  the  curves  to  go  through  one  at  zero  time  lag.  The  cross  correlation  has 
been  normalized  in  a  similar  manner, 


Ph  l(t) 


(h(t)-h)(L(t  +  T)-L) 
( h'L ') 


(4.7) 
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FlG.  4.17  Correlation  of  liftoff  height  h  and  flame  length  L  vs.  time.  Natural  gas, 
d  =  4.61  mm,  Us  =  39  m/s. 


The  flame  length  autocorrelation  shows  some  periodic  behavior,  corresponding 
to  the  burnout  of  large  structures.  For  reference,  the  calculated  time  scales  of  the 
flow  for  the  first  case  in  Table  4.6  are:  r$(/i)  =  6.4  ms,  T(,(L)  =  80  ms,  and  the 
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FlG.  4.18  Correlation  of  liftoff  height  h  and  flame  length  L  vs.  time.  Natural  gas, 
d  =  4.61  mm,  Us  =  53  m/s. 

convection  time  rc  from  h  to  L  is  approximately  120  ms.  For  the  second  case  the 
times  axe:  Ts(h)  =  9.0  ms,  r$(Z)  =  76  ms,  and  rc  =  120  ms.  For  estimating 
T(,{L)  and  rc,  experimental  results  of  Becker  and  Yamazaki  (1978)  were  used  in 
the  velocity  estimates  in  the  flame,  as  described  in  Sec.  A. 7.  The  convection  time 
is  defined  here  as  the  time  for  fluid  to  move  from  h  to  Z  traveling  at  the  mean 
centerline  velocity.  Because  there  is  considerable  uncertainty  in  the  estimates  of 
r$(Z)  and  rc,  they  are  intended  only  as  an  approximate  reference  time  to  compare 
to  the  correlations.  As  seen  in  Figs.  4.17  and  4.18,  the  time  scale  of  the  flame  length 
fluctuations  is  considerable  smaller  than  that  of  the  liftoff  height,  even  though  the 
large-scale  time  at  the  flame  tip  is  considerable  larger  than  at  h.  The  measurements 
of  the  two  quantities  are  somewhat  difficult  to  compare,  however,  because  the  flame 
length  shows  some  periodic  behavior  and  its  autocorrelation  has  such  a  different 
shape  from  that  of  the  liftoff  height. 


In  both  flames,  there  is  only  a  modest  correlation  between  h  and  Z,  which  is 
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negative,  and  peaks  at  a  time  in  the  vicinity  of  the  convection  time  from  h  to  L. 
The  most  important  observation,  however,  is  that  the  two  signals  are  not  strongly 
correlated,  indicating  that  the  flame  tip  behavior  is  not  responsible  for  the  observed 
long  correlation  times  in  the  liftoff  height.  This  is  the  result  that  was  expected. 
Considering  that  a  jet  is  nearly  a  parabolic  flow,  the  dynamics  far  downstream 
should  not  significantly  influence  the  upstream  regions.  Because  L  >  4h,  feedback 
from  the  flame  tip  to  the  liftoff  height  would  be  expected  to  be  weak. 


t  (ms ) 

FlG. 4.19  Correlation  of  liftoff  height  h  and  flame  length  L  vs.  time.  C2H6,  d  = 
3.88  mm,  Us  =  80  m/s. 


Results  similar  to  those  presented  for  natural  gas  were  obtained  with  C2H6 
using  a  3.88  mm  nozzle  as  shown  in  Fig.  4.19.  The  results  from  the  slightly  higher 
Reynolds  number  C2H6  flame  (Re^  =  20,000  for  the  C2H6  flame  vs.  12,000  for  the 
second  natural  gas  flame)  provide  further  evidence  that  the  dynamics  at  the  flame 
tip  are  not  responsible  for  the  behavior  of  h(t). 
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4.5  Simultaneous  Measurements  of  Liftoff  Height  and  Radiation 

As  shown  in  Figs.  4.10,  4.11,  and  4.15,  the  measured  values  of  axe  well 

correlated  with  £/,.  However,  the  fact  that  the  correlation  times  are  unchanged  by 
nozzle  orientation  (Fig.  4.13)  suggested  that  there  might  be  an  explanation  for  the 
collapse  other  than  the  influence  of  buoyancy. 

For  the  first  four  fuels  listed  in  Table  3.3,  the  radiation  from  soot  increases  in  the 
same  order  as  r |Q.  That  is,  C2H2  radiates  more  than  C2H4  which  radiates  more  than 
either  C2H6  or  natural  gas,  neither  of  which  produce  substantial  amounts  of  soot 
under  the  conditions  investigated.  Furthermore,  for  a  given  fuel,  e.g.,  C2H2,  flames 
at  low  non-dimensional  liftoff  heights  have  lower  Richardson  numbers  ( cf.  Eq.  4.6  on 
p.  56),  and  were  also  found  to  have  greater  radiation  from  soot.  Therefore,  for  the 
first  four  fuels  listed  in  Table  3.3,  there  is  a  one-to-one  correspondence  between  the 
Richardson  number  of  a  flame  and  the  fraction  of  heat  radiated  away  from  the  flame. 
It  was  conjectured  that  there  might  be  some  feedback  mechanism  whereby  the 
radiation  in  the  downstream  portion  of  the  flame,  where  fluid  mechanical  time  scales 
axe  larger  than  at  the  liftoff  height,  could  couple  with  the  liftoff  height  fluctuations 
and  cause  the  correlation  times  of  h(t)  to  increase.  This  would  then  explain  the 
collapse  of  ri/2/x«  with  shown  in  Fig.  4.11. 

The  measurements  reported  in  this  section  were  performed  prior  to  those  of 
Sec.  4.3.5  in  which  mixtures  of  H2  with  other  fuels  were  used  to  produce  flames 
with  low  Richardson  numbers  comparable  to  the  C2H2  flames,  but  without  the 
radiation  from  soot.  The  values  of  found  for  the  mixtures  were  comparable 

to  those  of  C2H2,  making  it  very  unlikely  that  the  large  correlation  times  of  C2H2 
result  from  the  radiation.  Nevertheless,  some  interesting  results  were  obtained  on 
the  simultaneous  measurement  of  liftoff  height  and  radiation,  which  are  reported  in 
this  section. 

It  was  noticed  in  visual  observations  of  the  C2H2  flames  that  when  the  flame 
base  dropped  to  a  low  position,  the  radiation  generally  increased,  presumably  be¬ 
cause  of  increased  soot  production.  To  verify  this,  simultaneous  measurements  were 
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made  of  the  liftoff  height  and  flame  radiation  at  a  position  downstream  of  the  flame 
base.  The  radiation  was  measured  with  the  same  photodiode  array  used  to  deter¬ 
mine  the  liftoff  height.  The  photodiode’s  spectral  response  to  radiation  is  not  flat. 
It  peaks  at  ~  700  nm  and  falls  off  rapidly  outside  the  visible  range.  Therefore, 
the  radiation  measurement  was  qualitative,  but  quite  adequate  for  the  purposes  of 
determining  the  degree  to  which  the  flame  radiation  and  liftoff  height  are  correlated. 

The  acetylene  flames  which  were  studied  are  listed  in  Table  4.7,  where  the 
intensity  I  is  in  arbitrary  units,  but  is  consistent  from  one  case  to  the  next.  The 
downstream  distance  at  which  I  was  measured  is  denoted  as  i;.  By  examining 
traces  of  intensity  all  along  the  flame,  it  was  found  that  the  radiation  at  different 
downstream  locations  fluctuates  together.  Therefore,  the  results  presented  are  not 
very  sensitive  to  the  choice  of  xj,  except  possibly  the  magnitude  of  the  light  in¬ 
tensity.  For  all  the  flames,  the  1.02  mm  nozzle  was  used.  For  the  first  three  cases 
listed,  a  neutral  density  filter  2.0  (each  0.3  is  one  stop,  or  a  reduction  in  light  in¬ 
tensity  by  a  factor  of  2)  was  used  to  reduce  the  signal  below  the  saturation  level  of 
the  photodiode  array.  The  fourth  case  had  low  enough  radiation  that  no  filter  was 
required. 

The  measurement  of  I  was  made  at  a  position  along  the  photodiode  array  some 
distance  downstream  from  that  of  /i,  typically  200  pixels.  Because  the  measurement 
of  intensity  at  any  given  pixel  of  the  array  is  the  integration  of  light  over  a  time 
period  A t  ending  at  the  measurement  time,  the  measurements  of  h  and  I  were 
not  exactly  simultaneous.  The  data  were  corrected  for  this  by  assuming  that  the 
measurement  of  /  occurred  at  a  time  lag  from  that  of  h,  equal  to  the  time  needed 
to  scan  from  the  pixel  corresponding  to  h  to  the  pixel  at  which  I  was  measured. 
For  the  acetylene  flames  this  was  a  significant  correction  to  the  data. 

Two  observations  can  be  made  from  Table  4.7.  First,  as  the  liftoff  height 
increases,  the  flame  radiation  decreases,  especially  near  blowout.  Second,  as  the 
liftoff  height  increases,  the  ratio  V  jl  increases.  This  is  because  at  the  high  liftoff 
heights,  the  flame  radiation  is  highly  intermittent,  occurring  in  short,  infrequent 
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Table  4.7:  Conditions  of  h  and  I  Measurements:  C2H2 

Us 

(m/s) 

h 

(mm) 

h' 

(mm) 

h 

XI 

(mm) 

/ 

(arb.) 

(arb.) 

166 

21.4 

0.90 

1.56 

60 

62 

11 

178 

24.9 

1.14 

1.81 

70 

56 

14 

199 

31.0 

1.65 

2.25 

80 

17 

6.7 

213 

35.7 

2.24 

2.60 

90 

0.62 

0.52 

FlG.  4.20  Liftoff  height  h  and  light  intensity  I  vs.  time.  C2H2,  d  =  1.02  mm,  Us  = 
178  m/s,  xj  =  70  mm. 


bursts. 


Portions  of  the  time  signals  of  h  and  I  from  the  second  case  in  Table  4.7  (Us  = 
178  m/s)  are  plotted  in  Fig.  4.20.  It  is  apparent  from  the  time  traces  that  when  the 
liftoff  height  drops,  the  radiation  tends  to  increase,  often  producing  a  sharp  peak 
in  the  I  signal. 

This  is  brought  out  more  clearly  in  a  plot  of  the  cross  correlations  of  h  and  /, 
which  were  computed  for  all  four  cases  in  Table  4.7.  Results  are  shown  in  Fig.  4.21 
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FlG.  4.21  Correlation  of  h  and  light  intensity  I  vs.  time.  Same  conditions  as 
Fig.  4.20. 

for  the  same  case  as  Fig.  4.20.  In  the  upper  right  comer  are  the  autocorrelation 
functions  of  h  and  of  I.  Here,  I  is  shown  to  have  approximately  the  same  correlation 
time  as  the  liftoff  height.  In  the  lower  half  of  the  figure  is  the  cross  correlation,  where 
a  positive  time  lag  is  defined  as  I  lagging  h.  Just  as  with  the  cross  correlation  of  h 
and  L  reported  in  the  previous  section,  ph  /  is  normalized  such  that  the  maximum 
possible  magnitude  is  one.  That  the  minimum  value  is  so  close  to  —1  is  an  indication 
that  the  two  signals  are  highly,  negatively  correlated. 

It  is  difficult  to  discern  from  the  data  by  how  long  the  change  in  liftoff  height 
leads  the  change  in  radiation,  because  the  temporal  resolution  is  insufficient.  The 
peak  in  the  cross  correlation  for  each  case  occurs  at  the  shortest  measured  time  lag, 
t’.e.,  in  the  same  scan  of  the  array,  with  the  time  lag  arising  from  xj  corresponding 
to  a  pixel  which  is  scanned  some  time  after  the  pixel  corresponding  to  h.  The 
convection  time  of  fluid  traveling  at  the  centerline  velocity  from  the  flame  base  to 
x/  was  estimated  as  described  in  Sec.  A. 7  to  be  ~  1.6  ms  for  the  data  shown  in 
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Figs.  4.20  and  4.21.  This  is  slightly  less  than  the  time  lag  for  the  peak  in  the  cross 
correlation  which  was  measured  to  be  1.9  ms.  However,  because  of  experimental 
resolution  difficulties,  with  the  time  between  samples  being  4  ms,  there  is  a  large 
uncertainty  in  the  measured  time  lag. 

A  possible  explanation  for  the  large,  negative  correlation  between  h  and  I  is 
that  when  the  flame  base  drops  to  a  low  value,  the  fuel,  which  is  not  as  well  mixed 
upstream,  has  a  longer  residence  time  in  the  hot  flame,  thereby  generating  more 
soot,  which  greatly  increases  the  radiation.  This  is  a  possible  mechanism  by  which 
the  change  in  liftoff  height  could  cause  the  radiation  to  increase. 

There  are  conceptually  two  possible  ways  in  which  the  flame  radiation  could 
cause  the  liftoff  height  to  change.  In  the  first  scenario,  the  increased  radiation  lowers 
the  temperature  in  the  flame,  which  therefore  changes  the  entrainment,  which  must 
act  to  keep  the  momentum  flux  constant  (assuming  a  low  Richardson  number). 
It  is  unlikely  that  this  could  affect  the  liftoff  height  dynamics,  however,  because 
the  radiation  occurs  downstream  of  the  flame  base,  and  the  largest  affect  on  the 
entrainment  should  also  be  downstream  of  the  flame  base. 

The  second  possibility  is  that  a  fraction  of  the  radiation  from  the  flame  is 
absorbed  by  the  gas  approaching  the  flame,  thereby  raising  its  temperature  and 
increasing  the  reaction  rate.  The  faster  reaction  would  cause  the  flame  to  be  sta¬ 
bilized  further  upstream,  providing  a  feedback  mechanism.  The  major  flaw  in  this 
argument  is  in  the  magnitude  of  the  possible  temperature  rise  in  the  reactants  en¬ 
tering  the  flame.  Modeling  the  radiation  as  a  point  source  of  strength  equal  to 
0.4  times  the  total  heat  release,  located  an  average  distance  of  1.5  h  away  from  the 
approaching  reactants,  talcing  the  velocity  along  the  ray  77  =  0.15,  and  using  an 
extremely  conservative  estimate  for  the  absorptivity  of  the  mixture  to  be  0.01  m_1, 
the  calculated  temperature  rise  is  less  than  0.01  K.  More  details  on  the  calculation 
are  provided  in  Sec.  A. 8.  This  temperature  rise  would  not  be  nearly  sufficient  for  a 
noticeable  change  in  the  chemical  kinetics  which  could  alter  the  liftoff  height. 


The  cross  correlation  functions  of  the  other  three  cases  listed  in  Table  4.7  are 
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similar  to  that  of  the  first  case  (Fig.  4.21).  The  minimum  correlation  value  always 
occurs  at  the  shortest  measurable  time  lag.  The  values  of  this  minimum  correlation 
for  the  four  cases  are:  —0.92,  —0.94,  —0.93,  and  —0.82.  The  correlation  does  drop 
somewhat  for  the  highest  lifted  flame.  Its  radiation  is  also  substantially  less  than 
the  other  cases  as  indicated  in  Table  4.7. 

To  compare  these  results  for  acetylene,  which  radiates  brightly  from  soot  pro¬ 
duction,  with  less  sooting  flames,  similar  experiments  were  conducted  using  natural 
gas,  C2H6,  and  C2H4.  Additionally,  measurements  of  I  were  made  at  more  than 
one  location.  The  conditions  of  the  experiments  involving  fuels  other  than  C2H2 
are  listed  in  Table  4.8. 


Table  4.8:  h  and  I  Measurements:  n.g.,  C2H6,  and  C2H4 

fuel 

d 

(mm) 

us 

(m/s) 

h 

(mm) 

XI 

(mm) 

n.g. 

6.25 

60 

181 

300,  400,  500 

c2h6 

3.88 

81 

150 

240,  340,  440 

C2H4 

2.16 

110 

73 

125,  200 

The  first  condition  listed  in  Table  4.8  involved  a  non-sooting  flame  of  natural 
gas  at  a  non-dimensional  liftoff  height  of  h  =  2.2.  To  determine  how  the  correlation 
of  h  and  I  is  affected  by  the  choice  of  downstream  location  x/,  measurements  of  I 
were  made  at  three  positions  simultaneously. 

The  cross  correlations  of  h  and  I  at  the  three  locations  are  shown  in  Fig.  4.22. 
The  peak  values  in  the  correlations  are  small,  indicating  that  h  does  not  substan¬ 
tially  affect  I  in  non-sooting  flames.  As  xj  is  increased,  the  correlation  diminishes, 
with  maximum  values  of  p(r)  being  0.12,  0.070,  and  0.036  for  x/  =  300,  400,  and 
500  mm,  respectively.  Additionally,  the  sign  of  the  correlation  is  opposite  of  that 
of  the  sooting  C2H2  flames.  This  is  a  result  of  the  completely  different  mechanisms 
of  radiation  of  the  two  fuels.  For  natural  gas,  the  radiation  comes  from  chemi¬ 
luminescence  of  C2  and  CH  radicals.  From  Fig.  4.22  it  appears  that  an  increase 
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FlG.  4.22  Correlation  of  h  and  light  intensity  I  vs.  time.  Natural  gas,  d  —  6.25  mm, 
Us  =  61  m/s,  h  =  181  mm,  x j  =  300,  400,  and  500  mm. 

in  liftoff  height  may  increase  the  total  amount  of  chemical  reaction  at  a  location 
downstream  of  h.  Alternatively,  the  increase  in  iifto^  height  may  change  the  en¬ 
trainment,  and  consequently  the  stoichiometry  at  the  downstream  location.  This 
would  affect  the  chemiluminescence,  which  is  very  sensitive  to  composition  ( e.g ., 
Chomiak  1990).  Because  the  response  of  the  photodiode  array  is  dependent  on  the 
wavelength  of  light,  a  change  in  the  spectral  content  of  the  radiation  would  change 
the  measurement. 

Because  the  correlation  functions  are  not  sharply  peaked,  it  is  difficult  to  iden¬ 
tify  a  time  lag  at  which  the  correlation  is  maximum.  There  does  appear  to  be  a 
slight  shift  to  longer  times  as  the  distance  downstream  xj  is  increased.  If  convection 
times  from  h  to  x/  axe  calculated  as  described  in  Sec.  A. 7,  these  timec  are  rc  =  13, 
25,  and  38  ms  for  the  three  values  of  x/.  The  correlation  curves  peak  at  larger  times 
than  these  in  general.  Also,  the  measured  shift  to  longer  times  with  increasing  x/ 
is  not  as  extreme  as  the  calculated  convection  times. 
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FlG.  4.23  Correlation  of  light  intensity  I  at  different  locations  xj  vs.  time.  Natural 
gas,  d  =  6.25  mm,  Us  —  61  m/s,  h  =  181  mm,  x/  =  300,  400,  and  500  mm. 

It  is  also  interesting  to  look  at  the  correlations  of  I  at  the  different  locations. 
This  gives  some  indication  of  the  convection  speed  of  the  flame  structures.  Cross 
correlations  of  Jj  J2  and  I\  I3  are  shown  in  Fig.  4.23  along  with  the  autocorrelation 
of  Jj,  where  the  subscripts  1,  2,  and  3  correspond  to  xj  =  300,  400,  and  500  mm, 
respectively. 

The  autocorrelation  of  I\  shows  quasi-periodic  behavior,  with  a  period  of 
100  ms.  For  comparison,  the  estimated  large-scale  time  at  x/ti  is  15  ms.  A  buoy¬ 
ancy  time  defined  as  Tb  =  (2  6(x)  Pool  9  Ap)1/2  is  ~  180  ms.  It  is  quite  possible  that 
this  oscillation  is  driven  by  buoyancy. 

By  having  measurements  of  light  intensity  at  several  locations,  it  is  possible  to 
determine  the  convection  times  from  one  location  to  the  next.  The  peaks  in  the 
correlations  occur  at  time  lags  of  38  and  68  ms  for  correlations  of  I\  J2  and  I\  J3, 
respectively.  These  are  longer  than  the  calculated  times  of  12  and  25  ms,  indicating 
that  the  structures  convect  at  speeds  lower  than  the  centerline  value. 
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Fig.  4.24  Correlation  of  h  and  light  intensity  I  vs.  time.  C2H6,  d  =  3.88  mm,  Us  = 
81  m/s,  h  =  150  mm,  x/  =  240,  340,  and  440  mm. 

Measurements  of  a  C2H6  flame,  which  produces  essentially  no  soot,  were  also 
made.  The  correlations  of  h  and  I  measured  at  three  different  locations  are  shown 
in  Fig.  4.24.  The  results  are  qualitatively  similar  to  those  of  natural  gas.  The 
estimated  convection  times  from  h  to  xj  are  9,  16,  and  24  ms.  The  peak  in  the 
correlation  at  xj  =  300  mm  occurs  at  ~  27  ms.  No  peaks  are  discernible  in  the 
correlations  at  the  other  two  locations. 

Finally,  measurements  were  made  with  C2H4,  which  does  have  radiation  from 
soot,  although  not  as  much  as  C2H2.  The  radiation  intensity  was  measured  at  two 
locations,  xj/h  ~  1.7  and  2.7.  The  lower  location  xjt\  is  below  the  region  of  intense 
soot  radiation,  while  x/i2  is  in  the  beginning  of  the  sooting  region.  The  cross  corre¬ 
lations  are  shown  in  Fig.  4.25.  At  the  lower  location,  there  is  very  little  correlation 
between  h  and  /.  At  x/i2,  however,  there  is  a  larger  (negative)  correlation.  The 
correlation  is  not  nearly  as  strong  as  that  of  C2H2,  however.  It  peaks  at  a  positive 
time  lag  of  9.4  ms,  which  is  close  to  the  calculated  convection  time  from  h  to  x/,2 
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FIG.  4.25  Correlation  of  h  and  light  intensity  I  vs.  time.  C2H4,  d  =  2.16  mm,  Us  = 
110  m/s,  h  =  74  mm,  2/  =  125  and  200  mm. 

of  9  dh  3  ms. 

While  it  was  difficult  to  determine  the  time  by  which  the  soot  production  lags 
the  liftoff  height  fluctuations  for  C2H2  because  of  limitations  in  the  resolution  of 
the  measurements,  it  is  clear  that  this  time  is  approximately  the  time  for  fluid  to 
convect  from  h  to  2/  for  C2H4.  In  all  likelihood,  this  is  also  the  case  for  C2H2 
flames  as  well. 


4.6  Discussion 

The  preceeding  results  on  the  temporal  behavior  of  the  liftoff  height,  the  simul¬ 
taneous  measurements  of  h  and  L,  and  the  measurements  of  h  and  I  still  leave  open 
the  question  of  what  determines  the  time  scale  of  the  fluctuations.  Although  the 
collapse  with  is  compelling,  the  lack  of  an  effect  of  nozzle  orientation  leaves  room 
for  speculation.  Additional  information  may  be  found  by  examining  what  happens 
to  the  magnitude  of  the  fluctuations  as  the  time  scale  changes. 
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FlG.  4.26  Relation  between  normalized  fluctuations  and  normalized  correlation  time 
at  a  fixed  Reynolds  number  of  Re*  =  22,000.  Filled  symbols:  n.g.,  C2H4, 
86%  H2  /  14%  C2H6. 

Figure  4.26  shows  Tj /2 / Tf,  vs.  W  jh  for  all  the  fuels  and  fuel  mixtures  for  which 
two  different  nozzles  were  used  (allowing  interpolation  to  Re ^  —  22,000).  This 
figure  shows  that  there  is  not  a  one-to-one  correspondence  between  Ti/2/ts  and 
h!  jh.  Although  there  is  a  decrease  in  h'  jh  and  an  increase  in  t1/2/t(>  for  the  fuels: 
natural  gas,  C2Hg,  and  C2H4,  the  values  of  h'  jh  are  approximately  the  same  for 
C2H4  and  C2H2,  for  which  Tj /2/t«  differs  by  a  factor  of  2.  Further,  for  the  mixtures 
of  H2  with  C2H6  and  C2H4,  values  of  h' jh  are  comparable  to  those  of  natural 
gas,  but  rj/2/r{  is  much  larger.  Any  explanation  of  the  variation  of  rj^/r*  with 
fuel  type  must  not  require  a  corresponding  change  in  the  relative  rms  level  of  the 
fluctuations. 

It  is  possible  that  some  subtle  property  of  the  fuels  is  responsible  for  the  varia¬ 
tion  in  One  candidate  is  the  relative  dependence  of  the  chemical  properties 

on  Y/Yst.  For  instance,  the  laminar  flame  speed  versus  Vyyrst  data  of  Egolfopoulos 
and  Law  (1991)  and  Egolfopoulos,  Zhu  and  Law  (1991)  plotted  in  Fig.  4.27  show  a 
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FlG.  4.27  Variation  of  laminar  flame  speed  with  mass  fraction.  Data  for  C2H6, 
C2H4,  and  C2H2  are  from  Egolfopoulos,  Zhu,  and  Law  (1991);  H2  data 
are  from  Egolfopoulos  and  Law  (1991). 

widening  of  the  values  of  Y/Ysl  for  which  S/Smax  is  large,  in  the  order  of  the  fuels: 
C2H6,  C2H4,  C2H2;  this  is  the  same  order  as  the  increase  in  The  flame 

speed  of  hydrogen  is  also  shown.  Its  curve  is  wider  than  any  of  the  other  fuels,  and 
it  is  shifted  toward  the  rich  side  of  stoichiometric.  The  shapes  of  the  S/Smax  vs. 
Y/Y st  for  mixtures  of  H2  with  C2H4  and  C2H6  have  not  been  measured,  but  they 
would  likely  fall  between  those  of  the  individual  components  of  the  mixture. 


4.7  Comparison  with  other  Measurements 

Birch  and  Hargrave  (1989)  measured  liftoff  heights  of  underexpanded  natural 
gas  flames.  Their  mean  results  have  already  been  discussed  in  Sec.  3.3.2.  They 
recorded  the  flame  fluctuations  using  a  video  camera  and  automatically  searched 
for  the  liftoff  height  along  a  ray  close  to  the  jet  edge,  rather  than  by  looking  at 
the  entire  flame  base  as  was  done  in  the  current  experiments.  They  reported  the 
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autocorrelation  for  one  of  the  flow  conditions:  d  =  50.6  mm,  stagnation  pressure  = 
28.1  bar,  h  =  2670  mm,  h!  =  186  mm.  From  their  autocorrelation  curve,  rXj2  is 
found  to  be  120  ms. 


FIG.  4.28  Data  point  of  Birch  and  Hargrave  (1989)  (+)  plotted  on  graph  of  tx/2/ts 
vs.  £/,.  All  other  points  (open  symbols)  sire  at  Re/,  =  22,000. 

Because  of  the  jet  being  choked,  estimates  of  the  flow  conditions  at  the  liftoff 
height  are  more  difficult.  Calculating  Us  and  <4  by  approximating  the  expansion 
to  atmospheric  pressure  as  an  isentropic  process,  as  was  done  in  Sec.  3.3.2,  Us  = 
810  m/s  and  ds  =  114  mm.  If  the  flow  conditions  are  then  calculated  by  the 
same  methods  used  previously  for  subsonic  jets,  the  values  are  h  =  1.3,  £/,  =  0.23, 
Re/,  =  6.3  x  106,  and  =  5.5  ms.  The  low  value  of  £/,  for  this  natural  gas 
flame  results  from  the  large  Re0 0  (cf.  Eq.  4.6).  The  value  of  tx/2/ts  ~  22  falls 
approximately  on  the  curve  of  ri/2/r«  vs-  th,  as  shown  in  Fig.  4.28.  The  other  data 
points  in  this  figure  are  at  Reoo  =  22,000.  The  collapse  of  rj /2  / with  £/,  appears 
to  hold  even  for  values  of  Re much  greater  than  those  investigated  in  the  current 
experiments.  The  uncertainties  in  the  estimates  of  £/,  and  are  greater  than  those 
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of  the  current  experiments,  however.  Birch  and  Hargrave  also  reported  a  quasi- 
periodic  aspect  to  the  h  vs.  t  signal,  with  a  period  of  ~  1.6  s,  which  is  ~  300  r^. 
They  attributed  this  to  the  passage  of  a  helical  vortex  through  the  flame  base,  which 
forces  the  low  point  of  the  flame  front  to  rotate  around  the  periphery  of  the  jet.  No 
such  behavior  was  found  in  the  current  measurements. 

Another  set  of  measurements  on  the  temporal  behavior  of  lifted  flames  is  that 
of  Chen  et  al.  (1989),  who  measured  the  radial  location  of  the  reaction  zone  using 
a  thin  ceramic  filament  and  a  photodiode  detector  to  record  the  radiation  from 
the  filament.  They  reported  values  of  the  rms  of  the  velocity  of  the  flame  position 
divided  by  the  rms  of  the  flame  position,  which  they  refer  to  as  a  strain  rate.  Using 
CH4  with  a  5  mm  nozzle,  at  an  exit  velocity  of  50  m/s  and  a  co-flow  velocity  of 
0.15  m/s,  their  mean  liftoff  height  was  100  mm.  The  measured  value  of  their  strain 
rate  was  520  s-1  at  an  axial  location  of  120  mm.  Their  co-flow  velocity  was  low 
enough  not  to  significantly  affect  the  jet  spreading  rate  of  centerline  velocity  at  h , 
and  so  r6  can  be  estimated  from  free-jet  scaling  laws  to  be  ~  4.2  ms.  This  gives 
a  value  of  strain  rate  times  r$  of  0.46,  indicating  that  the  radial  fluctuations  of 
the  reaction  zone  downstream  of  the  liftoff  height  occur  on  time  scales  of  order  r*, 
much  smaller  than  the  time  scales  of  fluctuations  in  liftoff  height.  There  is  some 
difficulty  in  comparing  their  measured  strain  rates  with  values  of  rj /2  in  the  current 
experiments,  however. 


4.8  A  Simple  Model  for  the  Fluctuation  Time  Scale 

A  simple  model  against  which  to  compare  the  observed  long-time  behavior  of 
lifted  flames  is  presented  here.  It  is  based  on  the  idea  of  a  turbulent  flame  prop¬ 
agating  upstream  against  the  approaching  gas,  and  gives  no  consideration  to  the 
turbulent  structure.  It  is  not  intended  as  a  plausible  explanation  for  the  experi¬ 
mental  results,  but  rather  as  a  simple  analysis  based  on  the  mean  properties  of  the 
jet. 
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It  is  assumed  that  the  flame  base  lies  on  the  mean  stoichiometric  contour  and 
is  perturbed  about  its  mean  liftoff  distance  h,  but  remains  on  the  contour.  The 
local  velocity  at  the  perturbed  location  h  =  h  +  h'  is  given  by 

—  d?/ 

u(h)  =  u(h)  +  —  h'  ,  (4.8) 

dx 

where  h'  is  the  perturbation  and  the  derivative  is  taken  along  the  stoichiometric 
contour.  The  propagation  speed  of  the  flame  S#  will  also  be  changed  and  is  given 

by 

St(h)  =  St(h )  +  ^  h'  ,  (4.9) 


where,  at  the  equilibrium  position,  u(h)  =  St(h).  Then  h(t)  is  described  by 


(4.10) 


Assuming  that  ^  =  0  (constant  propagation  speed  along  the  stoichiometric  con¬ 
tour)  the  time  constant  of  the  exponential  return  to  the  equilibrium  position  pre¬ 
dicted  by  this  model  is 

=  1/(||)  ■  (4.11) 

This  derivative  is  evaluated  from 


a  =  CuUs  —  fu{t]) 

x 


where 


/„(•»)  =  e-c--’ 


(3.2a) 


Taking  the  derivative  along  Y  =  yst, 


is  evaluated  from  the  equation  for  mean  concentration, 


(4.12) 


Y  =  CY  —  fy(ri) 

x 


(3.3) 


where 


At  Ysth/ds  =  2,  rm  =  18  Tf, ,  indicating  that  this  model  does  predict  fluctuation 
times  considerably  larger  than  r«.  The  strong  dependence  of  rm  on  Ysth/ds  is  not 
realized  in  the  experiments.  However,  the  trend  is  correct;  as  Ysih/da  increases, 
ri/2 !t6  decreases,  at  least  for  C2H4  and  C2H2. 
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CHAPTER  5 

Imaging  of  Fuel  and  Flame  Structures 


Experiments  were  conducted  in  which  two-dimensional  images  of  both  fuel 
concentration  and  reaction  zones  were  obtained.  Methane  was  used  as  the  fuel, 
and  the  reaction  zones  were  identified  by  the  presence  of  CH,  which  is  a  reaction 
intermediate.  Methane  was  imaged  by  Raman  scattering  and  CH  by  fluorescence. 
These  experiments  were  conducted  at  the  Combustion  Research  Facility  of  Sandia 
National  Labs,  in  collaboration  with  Dr.  Robert  Schefer  and  Dr.  Mehdi  Namazian. 
The  goal  of  the  experiments  was  to  gain  insight  into  the  interaction  of  the  turbulent 
structure  of  the  jet  and  the  flame  zones,  thereby  providing  additional  understanding 
of  lifted  flames,  supplementing  the  global  measurements  of  liftoff  height  described 
in  previous  chapters. 

In  the  flow  condition  described  in  this  chapter,  the  fuel  was  premixed  with  air, 
with  an  initial  mass  fraction  of  fuel  in  the  jet  fluid  of  0.62.  The  nozzle  diameter 
was  5.4  mm,  Us  was  23.7  m/s,  and  the  mean  liftoff  height  was  found  to  be  86  mm, 
making  h  ~  1.6.  The  Reynolds  number  was  Re*  =  6,700. 


5.1  Apparatus  and  Instrumentation 

As  with  the  experiments  discussed  previously,  the  nozzle  was  designed  to  have 
turbulent  pipe  flow  at  the  exit;  its  length  to  diameter  ratio  was  200.  The  jet  emerged 
from  a  solid  cylindrical  body  50  mm  in  diameter,  a s  illustrated  in  Fig.  5.1.  The  two 
cameras  used  to  record  the  images  were  located  95  mm  from  the  jet  centerline,  and 
shields  protecting  the  cylindrical  mirrors  of  the  multipass  cell  were  70  mm  from  the 
centerline. 
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Fig.  5.1  Experimental  apparatus  used  in  imaging  of  CH4  and  CH  (top  view). 

The  imaging  technique  used  was  that  described  by  Namazian  et  al.  ( 1988).  The 
DIANA  facility  dye  laser  at  the  Sandia  Combustion  Research  Facility  was  modified 
to  produce  a  beam  with  two  lines  at  431.5  and  444  mm.  Using  a  multipass  cell, 
the  beam  was  formed  into  a  0.3  mm  thick  sheet  passing  through  the  jet  centerline. 
The  431.5  nm  beam  caused  CH  fluorescence  at  489  nm,  and  coincidentally  Raman 
scattering  from  CH4  also  at  489  nm.  To  separate  the  CH4  and  CH  contributions, 
the  444  nm  line  was  used  to  produce  a  Raman  signal  from  CH4  at  510  nm.  The 
signals  were  recorded  with  two  intensified  vidicon  detectors  located  on  opposite 
sides  of  the  laser  sheet,  one  camera  with  a  10  nm  filter  centered  at  510  nm  and 
the  other  camera  with  a  10  nm  filter  centered  at  490  nm.  The  cameras  were  gated 
to  stay  on  for  4  ps.  To  calibrate  the  images  of  both  cameras,  the  test  section  was 
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filled  with  methane  using  a  100  mm  orifice  and  the  scattering  signals  recorded.  By 
normalizing  the  measurement  of  each  camera  with  the  signal  for  pure  CH4  at  room 
temperature,  and  then  taking  the  difference  of  the  recorded  signals  from  the  two 
cameras,  CH4  and  CH  were  measured  separately.  Five  hundred  shots  were  taken  of 
the  flame,  and  an  additional  100  shots  of  the  non-reacting  jet  were  taken. 

The  normalized  CH4  concentration,  denoted  [CH4],  is  scaled  such  that  a  value 
of  1  represents  pure  fuel  at  ambient  temperature;  it  can  be  converted  to  number 
density  by  multiplying  by  2.5  x  1019  molecules/cm3.  The  relationship  of  the  CH 
signal  to  CH  number  density  depends  on  temperature  and  composition.  Because 
some  regions  very  near  the  flame  base  may  be  farther  from  chemical  equilibrium 
than  others,  the  CH  signal  is  used  as  a  marker  of  flame  position,  but  not  interpreted 
quantitatively  in  terms  of  number  density.  The  CH  signal,  normalized  by  the  signal 
for  pure  CH4  as  described  above,  is  denoted  [CH]. 

In  order  to  increase  the  resolution  of  the  images,  viewing  was  restricted  to  one 
side  of  the  flame  using  53  pixels  in  the  horizontal  (y)  direction  and  110  pixels  in  the 
vertical  (x)  direction.  This  resulted  in  pixel  spacings  corresponding  to  ~  0.5  mm  at 
the  image  plane.  In  the  process  of  aligning  the  two  images  before  subtracting  them 
(to  obtain  the  [CH]  measurement),  the  data  were  interpolated  and  stored  as  arrays 
of  90  by  90  elements. 


5.2  Observations  on  Flame  Structure 

The  simultaneous  measurements  of  fuel  concentration  and  reaction  zones  re¬ 
vealed  many  interesting  features.  One  observation  that  can  be  made  is  that  there 
is  no  single,  characteristic  reaction  zone  structure  of  lifted  flames.  Several  shots 
were  selected  for  presentation  to  illustrate  the  range  of  structure  types  seen  in  the 
images. 

Figure  5.2  shows  contour  plots  of  [CH4]  and  [CH]  for  one  of  the  shots.  The 
contour  levels  for  [CH4],  starting  from  the  outside,  are  0.05,  0.1,  0.15,  0.2,  and  0.3. 
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(a)  (b)  (c) 


FlG.  5.2  Instantaneous  shot  illustrating  curved  reaction  zones.  Contour  plots  of  (a) 
CH4  number  density,  (b)  [CH]  signal,  and  (c)  the  superposition  of  (a)  and 
(b).  The  contour  levels  for  [CH4]  are  0.05,  0.1,  0.15,  0.2,  and  0.3.  The 
levels  for  [CH]  are  0.3,  0.34,  0.38,  0.42,  0.46,  0.50,  0.54,  0.58,  0.62,  and 
0.66. 

For  reference,  the  stoichiometric  level  is  0.095  at  room  temperature.  There  are  ten 
contours  for  [CH],  beginning  at  0.3,  which  was  found  to  be  a  good  choice  for  a 
threshold  to  distinguish  between  burning  and  non-burning  regions.  Labeling  the 
individual  contours  was  omitted  to  keep  the  figures  from  being  too  cluttered. 

Figure  5.2  illustrates  one  of  the  more  common  types  of  flame  structure  in  which 
the  reaction  zone  lies  near  the  outer  edge  of  the  jet,  and  approximately  follows  the 
contour  of  [CH4]  =  0.5  along  a  curved  path,  showing  evidence  of  the  large  structure 
of  the  turbulence.  The  interface  is  sharp,  suggesting  that  the  combustion  is  taking 
place  in  a  diffusion-flame  like  manner,  with  essentially  pure  air  on  the  outer  side 
and  a  fuel-rich  mixture  on  the  inner  side  of  the  flame.  At  some  locations,  the  CH4 
concentration  on  the  outer  side  of  the  flame  is  not  zero,  however. 

In  other  instantaneous  images,  less  evidence  of  the  large  turbulent  structure  is 
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FlG.  5.3  Instantaneous  shot  illustrating  relatively  straight  reaction  zones,  (a)  [CH4], 
(b)  [CH4]  and  [CH],  (c)  cuts  across  flame;  solid  line  [CH4],  vertical  axis  is 
0, . . . ,  0.2;  dashed  line  [CH],  vertical  axis  is  0, . . . ,  1.  For  contour  levels  see 
caption  of  Fig.  5.2. 

apparent.  Figure  5.3  illustrates  this  structure  type,  where  the  reaction  zone  is  rather 
straight  with  small-scale  convolutions,  again  following  the  [CH4]  =  0.5  contour 
reasonably  well.  As  with  Fig.  5.2,  the  combustion  resembles  that  of  a  diffusion 
flame.  Four  radial  cuts  have  been  plotted  in  (c)  to  highlight  the  variation  of  fuel 
concentration  moving  through  the  reaction  zone.  There  is  some  unbumed  fuel  on 
the  outer  part  of  the  reaction  zone,  even  at  the  most  downstream  location  plotted. 

This  structure  differs  from  that  of  a  pure  diffusion  flame  in  which  no  fuel  is 
present  on  the  oxidizer  side  of  the  reaction  zone.  Computational  results  of  Puri 
et  al.  (1987)  are  shown  in  Fig.  5.4  for  a  counterflow  diffusion  flame  of  CH4,  with  a 
strain  rate  of  55  s-1.  The  concentrations  of  CH4  and  CH  are  plotted  as  normalized 
number  density  to  allow  direct  comparison  with  the  present  experimental  results,  at 
least  for  CH4.  The  strain  rate  in  the  calculation  is  lower  than  that  in  the  turbulent 
flames  studied  here,  where  d\  ~  175  s-1  (using  the  result  a\  ri0  ~  7.5,  with  rj0  ~  43 
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Fig.  5.4  Structure  of  the  reaction  zone  in  a  counterflow  diffusion  flame,  from  Puri 
et  al.  (1987),  for  a  strain  rate  of  55  s-1. 

ms,  as  discussed  in  Sec.  3.2).  The  strain  rate  will  affect  the  profiles,  in  particular 
the  width  of  the  reaction  zone.  However,  the  main  point  to  be  made  from  Fig.  5.4 
is  that,  in  a  strained  diffusion  flame,  [CH4]  drops  to  low  values  before  the  region  of 
high  [CH],  and  is  essentially  zero  on  the  oxidizer  side  of  the  reaction  zone. 

Not  all  of  the  images  suggest  that  combustion  is  occurring  at  the  interface 
between  fuel-rich  fluid  and  pure  air.  In  some  shots  there  is  a  large  region  of  burning 
at  the  base  of  the  flame,  suggestive  of  premixed  combustion  propagating  against  the 
oncoming  flow.  This  type  of  structure  is  found  in  Fig.  5.5,  where  the  burning  region 
at  the  flame  base  is  seen  to  occur  over  a  range  of  CH4  concentrations  from  0.05  to 
0.1.  Above  the  flame  base  the  flame  zone  is  thinner,  appearing  more  diffusion- like. 

In  other  images  such  as  Fig.  5.6  there  axe  disconnected  burning  structures.  In 
the  three-dimensional  flow  field,  there  might  be  a  connection  between  the  different 
regions  seen  in  the  image.  These  isolated  flame  structures  can  be  at  different  axial 
distances  as  in  Fig.  5.6,  or  partially  overlapping,  with  one  structure  closer  to  the 
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FlG. 5.5  Instantaneous  shot  resembling  premixed  combustion, 
see  caption  of  Fig.  5.2. 

centerline  than  the  other.  There  is  less  of  a  correlation  between 
and  the  reaction  zones  in  this  image. 


y  (mm) 


For  contour  levels 
the  [CH4]  contours 


5.3  Statistical  Measures 

Several  statistical  measures  of  the  structure  of  the  fuel  and  reaction  zones  were 
calculated  from  the  ensemble  of  the  images.  Many  of  these  required  knowledge  of 
the  liftoff  height,  which  was  found  from  the  [CH]  field.  From  each  image  both  the 
liftoff  height  h  and  the  radial  stabilization  location  y*  were  found  by  searching  for 
the  most  upstream  occurrence  of  [CH]  above  a  threshold  of  0.3.  To  discriminate 
from  noise  in  the  data,  only  groups  of  at  least  five  adjacent  array  elements  above  the 
threshold  qualified  for  finding  the  liftoff  height.  Here,  adjacent  array  elements  are 
defined  as  any  of  the  eight  elements  surrounding  a  given  element,  h  is  then  defined 
as  the  most  upstream  distance  within  the  group,  and  y*  as  the  radial  location  of 
maximum  [CH]  in  the  row  at  x  =  h. 
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FlG.  5.6  Instantaneous  shot  illustrating  separated  reaction  zones.  For  contour  levels 
see  caption  of  Fig.  5.2. 


As  a  quantitative  measure  of  the  flame  structure,  the  probability  of  having 
burning  (defined  as  [CH]  above  the  threshold)  in  the  vicinity  of  the  liftoff  position 
was  calculated.  For  each  image,  after  finding  the  liftoff  position  (both  axial  h  and 
radial  y/, ),  the  presence  or  absence  of  burning  was  recorded  versus  displacement 
(Ax,  Ay)  from  the  liftoff  position.  All  the  images  were  then  averaged  to  produce 
the  probability.  By  definition,  the  probability  of  having  burning  at  Ax  =  Ay  =  0 
is  1.  The  results  are  shown  in  Fig.  5.7  as  a  contour  plot,  where  contours  levels  are 
shown  starting  at  0.9  (inner-most  curve)  and  decreasing  to  0.1.  To  keep  the  plot 
from  being  too  cluttered,  only  the  levels  0.2  and  0.1  are  labeled. 

That  the  contours  are  more  elongated  in  the  vertical  direction  indicates  that  the 
reaction  zone  at  the  flame  base  is  most  likely  to  be  aligned  approximately  vertically. 
The  fact  that  the  probability  falls  off  rapidly  with  distance  results  from  the  reaction 
zones  being  relatively  thin,  and  having  a  wide  range  of  orientations.  In  this  averaged 
representation,  there  is  no  one  common  flame  shape  which  appears. 
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Fig.  5.7  Contours  of  probability  of  finding  reaction  ([CH]  >  0.3)  at  locations  dis¬ 
placed  from  the  liftoff  location  by  Aar  in  the  vertical  direction  and  Ay  in 
the  horizontal  direction.  Nine  contours  are  shown  from  0.9  to  0.1;  only  the 
lowest  two  contour  levels  are  labeled,  1:  0.1,  2:  0.2. 


As  with  the  temporal  measurements  discussed  in  Ch.  4,  probability  density 
functions  can  be  computed  from  the  data  in  the  imaging  experiments.  The  pdf 
of  h  from  the  imaging  data  is  shown  in  Fig.  5.8.  The  pdf  is  broader  than  that  of 
the  measurements  using  the  photodiode  array  (Fig.  4.1),  having  a  normalized  rms 
of  h! /h  =  0.11.  This  higher  rms  and  corresponding  wider  pdf  is  probably  a  result 
of  the  fact  that  in  the  imaging  experiments,  only  a  thin  slice  through  one  side  of 
the  flame  is  used  in  determining  h,  rather  than  the  entire  flame  base.  It  is  clear 
from  the  images  that  local  extinction  occurs  downstream  of  the  flame  base,  and  this 
probably  occurs  at  the  social  location  of  the  flame  base  as  well.  This  would  cause 
the  measured  liftoff  height  in  the  imaged  plane  to  fluctuate  more  than  the  average 
over  the  flame  base. 
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The  pdf  of  the  radial  flame  base  position  is  shown  in  Fig.  5.9,  where  y/,  has 
been  normalized  by  yst,  the  radial  position  of  mean  stoichiometric  concentration 
at  x  =  h,  as  determined  from  the  shots  of  a  non- reacting  jet  at  the  same  flow 
conditions.  The  peak  in  the  pdf  is  slightly  to  the  outside  (fuel  lean  side)  of  the 
stoichiometric  position  with  a  mean  value  of  yh/ys t  of  1.17.  As  with  the  axial 
fluctuations,  the  radial  fluctuations  are  large,  with  y'h/yh  =  0.23  and  y'h/h  =  0.033. 
Although  the  flame  position  never  reaches  the  centerline,  the  fluctuations  are  of  the 
order  of  the  jet  radial  extent. 


FiG.  5.10  Probability  density  function  of  [CH4]  at  the  flame  base  (x  =  h,  y  =  y/,), 
and  at  a  location  just  upstream  of  the  flame  base  (x  =  0.98  h,  y  =  y/,). 


The  pdf  of  methane  number  density  at  the  flame  base  is  plotted  in  Fig.  5.10, 
along  with  that  at  0.98  h,  where  temperatures  axe  lower.  The  mean  value  of  [CH4] 
at  h  is  0.054,  and  at  0.98  h  the  mean  is  only  slightly  higher  at  0.058.  Both  are 
substantially  less  than  the  stoichiometric  value  at  ambient  temperature  of  0.095. 
The  number  density  at  h  is  lowered  by  the  higher  temperatures  in  the  reaction 
zone.  It  is  quite  possible  that  heating  of  the  reactants  even  at  0.98  h,  which  is 
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~  1.4  mm  below  ft,  lowers  the  density  somewhat.  The  pdf  at  0.98  ft  is  somewhat 
broader  than  at  ft.  This  is  probably  a  result  of  the  fuel  concentration  contours  being 
convoluted,  so  that  moving  upstream  even  1.4  mm  can  result  in  being  in  a  region 
of  either  higher  or  lower  concentration.  Even  the  pdf  at  ft  is  quite  wide,  with  the 
rms  being  45%  of  the  mean. 

Noise  in  the  measurements  will  contribute  to  the  measured  rms,  and  so  its 
contribution  must  be  estimated.  For  this  estimate,  both  the  background  rms  level 
n'0  (at  [CH4]  =  0)  and  the  rms  with  the  test  section  filled  with  methane  n\  (at 
[CH4]  =  1)  were  used.  The  equation 

n'2{x.y)  =  n'^{x,y)  +  [CH4]2(x,y)  (n'^(x,y)  -  n'o(x,y)) 


was  then  used  to  estimate  the  noise  at  each  location  (a-,  y),  using  the  mean  value 
of  [CH4].  To  account  for  the  fluctuation  of  the  radial  flame  position,  the  noise  was 
integrated  over  its  pdf.  Pyh.  i.e.. 


n'2  =  r 

Jo 


Pyh{yk)n'2(h.yh)dyh  . 


The  integrated  noise  n'2  was  found  to  be  0.014.  This  is  a  significant  contribution  to 
the  measured  rms  of  [CH4]  of  0.024.  Accepting  this  estimate  of  the  noise,  the  rms 
of  the  number  density  is  (0.0242  —  0.0142)1/2  ~  0.020,  making  the  number  density 
pdf  0.020/0.024  =  0.83  times  as  wide  as  the  one  plotted  in  Fig.  5.10.  While  there  is 
significant  noise  in  the  data,  after  accounting  for  it.  there  is  still  a  large  variation 
of  [CH4]  at  the  flame  base. 


The  large  variation  is  a  result  of  the  wide  range  of  conditions  present  at  the 
flame  base  at  different  times.  It  is  not  possible  to  determine  whether  the  mass 
fraction  of  CH4  or  the  temperature  at  the  flame  base  varies  more,  and  therefore  is 
the  principal  cause  of  the  variation  in  number  density  of  CH4. 


From  looking  at  the  images,  it  appears  that  the  reaction  zones  follow  contours 
of  [CH4]  reasonably  well.  To  get  some  quantitative  measure  of  this,  statistics  on 
the  [CH4]  level  in  the  center  of  the  reaction  zone  at  each  downstream  location  x 
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FlG.  5.11  Mean  [CH4]  number  density  at  the  location  of  peak  [CH]  vs.  distance 
downstream  of  h. 

above  h  were  computed.  Here,  the  center  of  the  reaction  zone  is  defined  as  the 
radial  location  of  maximum  [CH],  subject  to  the  level  being  above  the  threshold  for 
finding  h.  The  mean  and  rms  values  of  [CH4]  vs.  ( x  —  h)  are  shown  in  Fig.  5.11. 
The  mean  [CH4]  level  in  the  reaction  zone  is  0.54  at  the  base  of  the  flame,  decreases 
over  the  next  15  mm  to  0.4,  then  remains  at  0.4  for  at  least  the  next  15  mm.  It  is 
not  surprising  that  the  mean  level  is  below  the  stoichiometric  level  of  0.95,  because 
of  the  high  temperatures  in  the  reaction  zone.  In  fact,  it  is  rather  surprising  that 
the  number  density  is  as  high  as  it  is  even  downstream  of  the  liftoff  height. 
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CHAPTER  6 

Conclusions  and  Summary 


6.1  Mean  Liftoff  Heights 

Natural  gas,  C2H6,  and  C2H4  exhibit  a  nearly  linear  relationship  between  exit 
velocity  and  liftoff  height,  as  had  been  found  by  Kalghatgi  (1984)  for  several  fuels 
including  CH4  and  C2H4.  There  is,  however,  a  slight  decrease  in  the  required 
velocity  needed  to  achieve  a  given  liftoff  height  as  the  nozzle  diameter  is  increased. 
In  non-dimensional  terms,  Ysth/ds  is  nearly  proportional  to  Y^rxUs/ds,  although 
there  is  a  slight  decrease  in  YgtrxUs/ds  at  a  fixed  value  of  Ysth/ds  as  the  Reynolds 
number  Re^  is  increased.  For  acetylene  the  relationship  between  Us  and  h  exhibits 
a  large,  negative  extrapolated  intercept  to  the  h  axis,  at  least  for  a  1.0  mm  nozzle,  for 
which  Reoo  is  less  than  15,000.  With  larger  nozzles  (and  larger  Reynolds  numbers), 
the  extrapolated  (non-dimensional)  intercept  is  reduced. 

A  strain-rate  model  based  on  far -field  behavior  of  turbulent  jets  provides  a 
plausible  explanation  for  the  linear  h  vs.  Us  relationship.  It  also  accounts  for  the 
variation  in  the  slope  of  the  line  with  fuel  type. 

Experiments  in  which  air  was  premixed  with  the  fuel  prior  to  exiting  the  nozzle 
showed  that  the  slope  of  the  ( Y9i/Y0)h/da  vs.  (Tst/T’0)2C/s/ds  line  increases  as  Y0 
is  decreased.  However,  the  intercept  on  the  ( Yai/Y0)h/ds  axis  decreases  as  Y0  is 
decreased,  creating  a  fan  pattern  in  the  plot.  This  is  in  contrast  with  expectations 
based  on  far-field  seeding  of  turbulent  jets,  which  predicts  a  collapse  of  the  data  at 
any  value  of  Y0.  The  discrepancy  is  attributed  to  near- field  effects,  which  are  more 
prominent  when  Y0  is  low  because  such  flames  are  restricted  to  lower  values  of  h/d^ 
before  blowout. 
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6.2  Fluctuations  of  the  Liftoff  Height 

The  spatial  scales  of  fluctuations  in  h  are  of  the  order  of  the  local  large  scale  of 
the  jet.  Values  of  h' /h  fall  between  0.05  and  0.08.  There  is  an  increase  in  h' /h  as 
flames  approach  blowout,  or  h  becomes  large.  There  is  a  modest  variation  of  h'/h 
with  fuel  type,  the  reason  for  which  is  not  known. 

Time-resolved  measurements  of  the  liftoff  height  revealed  that  fluctuations  oc¬ 
cur  on  time  scales  considerably  longer  than  the  local  jet  time  r$.  Using  different 
fuels  and  flow  conditions  to  achieve  different  r$,  it  was  found  that  the  normalized 
correlation  time  t^/t*  is  not  a  constant.  Rather,  for  all  the  conditions  inves¬ 
tigated,  there  is  a  good  collapse  with  the  local  Richardson  number  in  which 
ti/2/t6  increases  as  decreases.  Both  Ts  and  £/,  depend  largely  on  the  effective 
fuel  chemical  time  t\0.  The  fuels  (and  fuel  mixtures)  studied,  in  order  of  increasing 
Ti/2/n,  are:  natural  gas,  C2H6,  C2H4,  50%  C2H2  /  50%  C2H6,  50%  H2  /  50% 
C2H4,  86%  H2  /  14%  C2H6,  C2H2,  and  50%  H2  /  50%  C2H2.  Experiments  in  which 
the  nozzles  were  oriented  horizontally  showed  no  difference  in  the  temporal  behavior 
of  h,  suggesting  that  it  may  not  be  the  influence  of  buoyancy  which  is  responsible 
for  the  variation  of  tx/2/t6. 

Feedback  from  the  flame  tip  is  not  responsible  for  the  long  correlation  times 
nor  for  the  variation  of  Ti/2/t6  with  fuel  type.  This  was  confirmed  by  simultane¬ 
ously  measuring  h(t)  and  the  flame  length  L(t)  and  finding  that  they  are  nearly 
uncorrelated  in  time. 

There  is  a  monotonic  relationship  between  the  measured  values  of  tx  /2  /  r#  and 
the  visible  radiation  from  soot  of  the  above  listed  fuels,  other  than  the  mixtures 
with  hydrogen.  In  particular,  C2H2  radiates  brilliantly  and  has  substantially  higher 
values  of  r,/2/T$  than  the  other  fuels.  To  investigate  whether  the  radiation  and 
liftoff  height  are  correlated  in  time,  simultaneous  measurements  were  made  of  h 
and  the  visible  light  intensity  I  at  a  downstream  location  in  the  flames.  For  C2H2 
flames,  radiation  from  the  flame  is  strongly,  negatively  correlated  with  h.  h  and  /  are 
moderately  correlated  for  C2H4,  which  produces  less  soot,  and  nearly  uncorrelated 
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for  C2H6  and  natural  ga s,  which  produce  essentially  no  soot  under  the  conditions 
studied. 

This  correlation  between  h  and  /  is  not  responsible  for  the  variation  in  Ti/2/ts 
with  fuel  type,  however.  Measurements  with  fast,  non-sooting  mixtures  of  H2  with 
C2H6  and  C2H4,  with  values  of  r*  and  close  to  those  of  C2H2,  revealed  comparable 
values  of  Ti/2/t6. 

There  is  not  a-one-to  one  correspondence  between  Ti/2/ts  and  h'/h.  This 
eliminates  any  explanation  of  the  variation  of  /2 / 77,  with  fuel  type  which  would 
require  a  corresponding  change  in  h'/h. 

It  is  possible  that  the  variation  in  ri/2/77  with  fuel  type  results  from  details  of 
the  chemical  properties  of  the  fuels,  rather  than  £/,.  In  particular,  the  dependence 
of  the  reaction  rate  on  composition  may  change  the  structure  at  the  flame  base  and 
hence  the  flame  stability.  There  is  a  correspondence  between  Ti/2/t6  measured  in 
these  experiments  and  the  width  of  the  curve  of  laminar  flame  speed  vs.  Y /  Yst . 


6.3  Imaging  of  Fuel  and  Reaction  Zones 

Planar  images  of  methane  number  density  and  CH  (to  mark  the  reaction  zones) 
were  taken  of  a  flame  lifted  to  h  =  1.6  at  a  Reynolds  number  of  Re/,  =  6,000.  In 
some  of  the  images  the  large-scale  structure  of  the  turbulence  is  quite  evident,  with 
large,  curved  reaction  zones,  located  near  the  outer  part  of  the  jet.  Sometimes  there 
there  are  multiple,  separated  reaction  zones  suggesting  local  extinction.  In  other 
images  the  reaction  zones  are  rather  straight,  but  still  located  near  the  outer  edge  of 
the  jet,  suggesting  diffusion-like  combustion.  Not  all  of  the  images  appear  to  have 
combustion  taking  place  in  a  diffusion-like  maimer,  however.  There  is  sometimes  a 
large  region  at  the  base  of  the  reaction  zone  which  appears  to  be  premixed  flame 
propagation,  followed  by  more  diffusion-like  reaction  zones  downstream.  Even  in 
the  diffusion-like  structures,  the  concentration  of  CH4  was  not  zero  on  the  air  side 
of  the  reaction  zone. 
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Fluctuations  in  both  the  axial  (h)  and  radial  (y/,)  flame  stabilization  location 
were  found  to  be  large,  on  the  order  of  the  local  jet  width  6,  with  h' /h  =  0.011  and 
y'h/h  =  0.033.  The  mean  value  of  yn  was  found  to  be  ~  1.2  yst,  where  yst  is  the 
location  of  mean  stoichiometric  concentration  in  the  non-reacting  jet.  At  the  flame 
base,  the  mean  value  of  [CH4],  the  number  density  normalized  by  that  of  pure  fuel 
at  ambient  temperature,  was  found  to  be  0.054.  Large  variations  in  [CH4]  at  the 
flame  base  were  observed  in  the  individual  shots. 


6.4  Summary 

Measurements  of  mean  liftoff  height  confirmed  the  nearly  linear  dependence  of 
h  on  Us,  with  a  slight  increase  in  h  with  Reynolds  number.  This  linear  dependence 
is  in  agreement  with  a  strain-rate  model  based  on  the  far-field  scaling  of  turbulent 
jets.  When  increased  amounts  of  air  are  added  to  the  nozzle  fluid,  the  slope  of  the  h 
vs.  Us  line  increases  faster  than  predicted  by  far-field  jet  scaling.  The  spatial  scales 
of  the  fluctuations  of  h  were  found  to  be  of  the  order  of  the  jet  large  scale,  while  the 
time  scales  were  found  to  be  much  larger  than  the  local  large-scale  time.  The  non- 
dimensional  fluctuation  time  varies  with  fuel  type,  and  collapses  with  the  Richardson 
number  at  the  flame  base,  £*,.  There  is  no  change  in  the  measured  fluctuation  time 
scale  when  the  nozzles  are  oriented  horizontally,  suggesting  that  the  fluctuation  time 
scale  may  depend  on  some  parameter  other  than  £/, .  The  variation  is  not  caused  by 
feedback  from  the  flame  tip  nor  by  feedback  from  radiation.  It  may  be  due  to  details 
of  the  chemistry  of  the  fuels.  Experiments  in  which  fuel  concentration  and  reaction 
zones  were  imaged  showed  a  large  range  of  structure  types.  Large  variations  in  the 
axial  and  radial  flame  stabilization  location  and  fuel  number  density  at  the  flame 
base  were  observed. 
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APPENDIX  A 

Calculation  of  Fluid  Flow  Quantities 


A.l  Momentum  Velocity  and  Momentum  Diameter 

From  the  calculated  mean  velocity  exiting  the  nozzle  U0,  obtained  as  described 
in  Sec.  2.1.1,  the  equivalent  source  velocity  Ua  was  calculated,  along  with  the  equiv¬ 
alent  source  diameter  ds.  Us  and  ds  are  defined  in  terms  of  the  momentum  flux  JQ 
and  mass  flux  m0  exiting  the  nozzle  by  the  solution  of  the  equations 


rh0  =  poo  ~  d\  Us 
Jo  = 


(A.l) 


To  calculate  Us  and  ds,  the  mean  velocity  profile  for  turbulent  flow  in  a  pipe 
as  reported  by  Nikuradse  (1932)  was  used, 


-  =  (»--)’• 
“cl  \  To) 


(A-2) 


where  uci  is  the  centerline  velocity,  r  is  the  radial  coordinate,  and  r0  is  the  noz¬ 
zle  radius.  The  exponent  n  depends  on  the  pipe  flow  Reynolds  number  Re0  = 
( p0  U0d)/p0,  as  reported  by  Nikuradse. 


By  integrating  the  velocity  profiles,  U9  and  ds  are  foirnd  to  be 

U9  (n+l)(n/2  +  l)2 
U0  ~  (2n  +  l) 


(A-3) 


<k  =  (hSn  (  <2n  +  1)  V/2  (A  4) 

d  \(n  +  l)(n/2  + 1)2/  ■ 

To  a  reasonable  approximation  for  the  Reynolds  numbers  investigated  (10,000  < 
Re0  <  50,000), 

Ua  =  1.02  U0 

(A. 5) 

d,  =  0.99(  pJp^'U. 


(A.4) 
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A.2  Mean  Strain  Rate 

The  mean  strain  rate  in  the  far  field  of  a  turbulent  jet  can  be  obtained  from 
the  mean  velocity  profile, 

«  =  C.V.±-f.(t))  .  (3-2) 

x 

From  the  continuity  equation,  the  radial  velocity  is  obtained, 


1  r^A- 

v  = - /  r  —  dr  . 

r  J0  dx 


(A.6) 


Using  Eq.  3.2,  this  becomes 


/.(■»)  -  i  jT  V  /.«)  •  (A.7) 


Partial  derivatives  of  the  velocity  components  are  then  given  by 

tt  =  c.  ^  r.w 


dv_ 


Us  dg  /  2 


dv_  Uads 


The  strain  rates  are  given  by 


= 


-  {v2  Kin)  +  V  fn(v))  =  Vfc 
(n  f'u(r))  +  ^  • 


(A-8) 


(Trr  — 


1  ( du  dv\ 

°,T  =  5  V  a? +  &  j 

*«  =  2  =  (/.(*)-£  j[V.0»)*>)  ■ 


Using  a  Gaussian  for  the  velocity  function, 


fu(v)  =  e 


-  .-cl*’ 


(A-9) 


(A.10) 
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these  become 


The  principal  strain  rates  in  the  x-r  plane  are 


+  *lr 


(A.ll) 


(A. 12) 


The  most  compressive  strain  rate,  am,  is  the  minimum  of  aQ,  ag,  and  agg.  The 
most  extensive  strain  rate,  a\,  is  the  maximum  of  the  three. 


A.3  Velocity,  Concentration,  Density  just  below  Flame  Base 

The  centerline  velocity  in  the  far  field  of  a  non-buoyant,  non-reacting  jet  is 
given  by 

Uc i  =  CUUS-  .  (A. 13) 

x 

This  equation  relies  on  the  fact  that  as  a:  — ♦  oo,  p  — ►  poo-  Provided  that  h  is 
out  of  the  near  field,  and  far  enough  downstream  that  p  ~  p^,  Eq.  A. 13  provides 
an  adequate  estimate  for  Uc\  at  x  =  h.  This  is  sufficiently  accurate  for  all  of  the 
conditions  presented  here,  except  those  using  a  mixture  of  86%  H2  /  14%  C2H6, 
for  which  the  centerline  density  at  x  =  h  is  considerably  less  than  p0 0.  In  order 
to  obtain  a  more  accurate  estimate  for  Uc\  in  those  flames,  another  method  for 
estimating  Uc\,  as  well  as  Y'ci,  was  employed. 
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The  momentum  flux  J  and  mass  flux  of  fuel  rhf  remain  constant  at  any  axial 
location  x  below  the  flame.  The  time-averaged  values  of  rhf  and  J  axe  given  by 


rhf  =  J  ^pY  u  +  (pY)' u'  2  nr  dr 
J  =  J  ^pu2  +  ~pu'2  +  2uu'  p'  +  p'  u'2  2nrdr  . 


(A.14) 


The  p'  u'2  term  was  assumed  to  be  small,  and  was  neglected  in  the  momentum 
equation;  all  other  terms  were  estimated  from  measurements  in  turbulent  jets.  The 
radial  behavior  of  u  and  Y  were  assumed  to  be  given  by  the  similarity  functions 

fci  = 


=  fM 


=  Mn) 


Ud 

Vy^ 

Yd 

Y'u' 


(A. 15) 


=  fY'(v) 


=  fY'u'(v)  ■ 


Y*Ud 

It  was  assumed  that  the  statistics  of  u  and  Y  are  unaffected  by  the  initial  den¬ 
sity  ratio  Pol  poo-  The  functions  in  Eq.  A. 15  were  estimated  from  measurements  in 
turbulent  jets.  The  concentration  functions  /y  and  fy  were  taken  from  data  of 
Dowling  (1988).  The  mean  velocity  function  /„  was  taken  from  Chen  and  Rodi 
(1980)  and  the  rms  velocity  function  /„-  from  Wygnanski  and  Fiedler  (1969).  The 
cross  correlation  function  of  velocity  and  concentration  / y'u •  was  obtained  from 
measurements  of  So  et  al.  (1990)  who  reported  velocity-density  correlations  in  a  jet 
with  an  initial  density  ratio  of  p0/Poo  =  0.64. 

The  density  as  a  function  of  mass  fraction  is  given  by 

~pZ  =  ^(MWoo/MW.-lJ  +  l  *  (A'16) 

To  obtain  estimates  for  the  terms  in  Eq.  A.14  involving  p' ,  Eq.  A. 16  was  linearized 


about  Y: 


»  ^-Y' 

~  d YY  ' 


(A.17) 
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Equations  A.  14  can  then  be  combined  with  the  similarity  functions  in  Eq.  A.  15 
to  obtain  two  simultaneous  equations  for  Uc\  and  Fch 


1/8  (7)  *  = r"  t:  h  + y“  ^  3F  /* 


(A. 18) 


+  fl)  +  Ya  —  £f.fr: 

Poo  Ul 


rjdrj 


Equations  A. 18  were  solved  for  ycj  and  Uc  1  at  the  location  x  =  h.  When 
p/ Poo  —  1?  these  equations  reduce  to  the  asymptotic  forms 


(A.19) 


x 

with  the  values  of  C„  =  6.17  and  Cy  =  5.92  foimd  from  the  solution  of  Eqs.  A. 18. 
This  value  of  Cu  is  very  close  to  the  value  reported  by  Chen  and  Rodi  (1980)  of 
6.2.  The  value  for  Cy  obtained  from  Eq.  A.  18  is  somewhat  higher  than  previously 
measured  values,  however.  Chen  and  Rodi  recommend  5.0,  and  Dowling  (1988) 
finds  Cy  =  5.0  at  Re  =  5,000  and  4.7  at  Re  =  16,000. 


A. 4  Non-Dimensional  Liftoff  Height 

The  definition  chosen  for  the  non-dimensional  liftoff  height  relies  on  the  con¬ 
centration  field  of  the  jet.  Because  the  fuels  studied  have  maximum  reaction  rates 
at  compositions  near  the  stoichiometric  value  Kst,  it  was  decided  that  selecting  a 
parameter  that  would  be  a  function  of  r^st,  the  value  of  rj  for  which  Y  =  Fst,  would 
be  the  best  choice. 

In  the  fax  field  of  the  jet,  the  mean  concentration  field  is  given  by 

V  =  CyK-Mv)- 

X 


(3.3) 


109 


From  this  equation,  it  is  clear  that  ( Yst/Y0)h/ds  is  a  function  of  T)sl,  and  is  the 
correct  non-dimensional  liftoff  height.  This  is  the  exact  parameter  used  in  Ch.  3  for 
the  non-dimensional  liftoff  height. 

For  all  of  the  conditions  studied  except  those  using  a  mixture  of  86%  H2  /  14% 
C2H6,  this  definition  would  have  sufficed,  because  Eq.  3.3  is  sufficiently  accurate. 
However,  to  be  consistent  with  the  calculation  of  Uc\  and  described  in  Sec.  A. 3, 
a  slightly  different  definition  of  the  non-dimensional  liftoff  height  h  was  used: 

h  =  CY  y-  •  (A.20) 

Using  this  definition,  h.  is  a  function  of  r)st,  and,  as  pc\  — *•  poo ,  h  — *  ( Yst/Y0)h/ds . 


A.5  Reynolds  Number  at  the  Liftoff  Height 

The  Reynolds  number  in  the  cold  jet,  just  below  the  flame  base  (Re^)  was 
determined  for  each  experimental  condition  by  neglecting  my  possible  influence 
from  the  flame  downstream.  In  the  far  field  of  the  jet,  where  sufficient  air  has 
mixed  with  the  jet  fluid  (fuel)  such  that  p  «  poo  and  p  ~  poo,  the  Reynolds  number 
becomes  independent  of  x,  and  is  Re oo  =  Uads/ i/qo.  To  get  a  slightly  better  estimate 
of  Re*,,  the  values  of  Uc \  and  Yc\  calculated  as  described  in  Sec.  A.l  were  used. 

The  local  Reynolds  number  is  proportional  to  Uc\  6/u.  Noting  that  6  oc  x,  the 
precise  definition  for  Re^  is 

He*  -  (£)  U-ir  ■  <A-21> 

where  u  =  p/p  is  evaluated  for  the  mixing  cup  concentration 

Y  =  ~  ,  (A.22) 

m 

where  rhf  is  the  mass  flux  of  fuel  and  m  is  the  total  mass  flux,  both  evaluated 
at  x  =  h.  Dividing  by  the  velocity  decay  constant  C„  ensures  that  as  x  — ♦  00, 
He  h  ^  He  00  * 


no 


A.6  Richardson  Number  at  the  Liftoff  Height 


The  Richardson  number  at  the  flame  base  Rik  is  defined  as  the  ratio  of  a 
buoyancy  force  B  in  the  vicinity  of  the  flame  base,  to  the  momentum  flux  entering 
the  flame  base  J  =  (7r/4)d^f/s2.  The  buoyancy  force  was  calculated  from  the  liftoff 
height  h  to  a  distance  ah.  This  distance  was  taken  as  the  size  of  one  large  structure, 
making  a  equal  to  1  +  Cg,  where  the  jet  spreading  rate,  Cs  =  6/x  «  0.44.  Then, 


B  = 


-  P)dV  , 


(A. 23) 


and  Rift  ~  R/ J.  In  order  to  facilitate  comparison  of  this  definition  of  a  Richardson 
number  relevant  to  the  liftoff  height  dynamics  with  definitions  for  the  flame  tip 
dynamics,  a  constant  factor  was  added  to  this  definition. 


Becker  and  Liang  (1978)  used  Ril  =  (g  L3)/{U^  d%)  as  the  Richardson  number 
for  the  flame  length,  noting  that  p^  —  p  «  pQ 0  in  hydrocarbon  flames,  and  that 
the  volume  ~  L3 .  The  actual  volume  of  a  cone  of  height  L  and  spreading  rate  Cs 
is  Kone  =  (tt/12 )C|  L3  =  0.051  L3.  Therefore,  the  definition  of  Ri/,  in  the  current 
experiments  is 

BJa  =  7T  Cl  a3  Ufdl  '  (A‘24) 

Defining  A p  is  the  spatially  averaged  value  of  A p  in  the  first  flame  structure, 


—  _  h  /05 C<  z(po0  -p)2nrdrdx 
I-hh  J05C(X2irrdrdx 

Eq.  A. 23  becomes 


(A. 25) 


Ri*  =  (a3  -  l)/«3(jA^  hz)HVl<ti)  .  (A.26) 

The  factor  a3  —  1  is  included  because  only  the  region  above  h  is  at  low  density;  the 
factor  of  1/a3  is  included  because  without  it,  this  would  actually  be  a  Richardson 
number  at  i  =  ah,  rather  than  at  x  =  h. 


Ill 


The  remaining  difficulty  in  this  definition  was  to  determine  p  in  the  flame  re¬ 
gion.  For  that,  the  data  of  Becker  and  Yamazaki  (1978)  for  propane  flames  attached 
to  the  nozzle  were  used  in  the  following  manner.  By  measuring  dynamic  pressure 
and  temperature,  they  calculated  the  mixing  cup  density  p  =  m/Q,  where  m  is  the 
mass  flux  and  Q  is  the  volume  flux  at  the  downstream  measuring  station.  They 
found  that  to  a  good  approximation,  p/ Poo  is  a  function  of  the  local  Richardson  num¬ 
ber  and  independent  of  the  distance  along  the  flame  x/L.  In  particular,  they  found 
that  for  £x  <  2,  p/poc  ~  0.20,  where  their  definition  of  is  ( x/ds)g/(Ug  d f^)1/3. 

Although  their  measurements  were  on  attached  flames  (often  using  hydrogen 
stabilization),  their  data  were  used  to  approximate  the  density  in  lifted  flames. 
Account  was  taken  of  the  properties  of  each  fuel  in  applying  their  results  from 
propane  flames.  Because  all  of  the  flames  studied  in  the  current  experiments  had 
ix  <  2,  with  the  distance  ah  used  for  r,  the  result  that  ( p/Poo )  =  0.2  was  used. 
This  was  converted  into  a  mean  temperature  rise  by 


AT 

AT/ 


Too  f  Poo  _  -A 

AT/  \  p  ) 


=  0.6  , 


(A. 27) 


where  AT/  is  the  adiabatic  temperature  rise  for  a  stoichiometric  mixture  of  fuel 
and  air.  Approximating  A p  as  p^  —  p  gives 


A  p  =  poo 


Too  MW  \ 
Too  +  ATMWoJ 


(A. 28) 


The  average  molecular  weight  in  the  flame  MW  was  estimated  by  solving  for  the 
mass  fractions  of  product,  fuel,  and  oxidizer  from  the  equations 

Yp  =  AT/AT/ 

y,t  Tp  "b  ^  /  =  ^hfuei/^hair 


Yp  +  Yf  +  Y0  =  1  . 


(A. 29) 
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A. 7  Convection  Time  Estimates 

Estimates  of  convection  times  in  the  flames  were  made  to  compare  with  the 
simultaneous  measurements  of  h  and  I  at  27,  as  well  as  the  measurements  of  h  and 
L.  For  the  h  and  I  measurements,  estimates  of  the  convection  times  from  ft  to  r/ 
and  from  xjj  to  07,2,  etc.  were  made.  For  the  h  and  L  experiments,  estimates  of  the 
convection  time  from  h  to  L  were  made.  Centerline  velocities  were  estimated  at  the 
locations  h,  xj  and  L\  these  velocities  were  used  in  estimating  both  the  convection 
times  and  the  local  large  scale  times  r* . 

For  the  acetylene  flames  listed  in  Table  4.7,  for  which  the  Richardson  numbers 
at  the  measurement  locations  xj  were  all  low,  it  was  assumed  that  the  centerline 
velocity  is  reasonably  well  approximated  by  the  isothermal  result, 

Uel  =  Cu  Us  —  .  (A. 30) 

x 


The  convection  time  of  fluid  traveling  at  the  centerline  velocity  is  then 

xj  ~  & 

Tc  ~~  2CuUsds  • 


(A- 31) 


Even  at  low  Richardson  numbers,  the  change  in  density  may  affect  the  centerline 
velocity.  However,  it  was  felt  that  this  method  would  provide  a  more  accurate 
estimate  of  Uc\  for  the  acetylene  flames,  than  the  method  described  below,  which 
was  used  for  the  other  fuels,  listed  in  Table  4.8. 


For  the  fuels  other  that  acetylene,  two  different  estimates  for  the  velocity  at  xj 
were  made.  Because  most  of  the  locations  were  not  too  far  downstream  of  h ,  the 
first  method  was  simply  to  use  Eq.  A. 30.  Some  of  the  locations  xj  were  far  enough 
downstream  of  the  flame  base  that  buoyancy  would  affect  the  velocities,  and  so  an¬ 
other  method  for  estimating  velocities  was  taken  from  Becker  and  Yamazaki  (1978). 
They  investigated  attached  flames,  but  their  results  were  used  as  an  approximation. 
Their  Eq.  38  can  be  written 

IT2,  =  [72  (ll)  —  (2.6  -I-  0.51  £)3  , 

V  1  )  Pc\ 


(A. 32) 
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where  they  have  defined  £  as  (g  x3/U*  d^)1/3.  To  obtain  the  centerline  density  pc\, 
their  Fig.  9  was  used,  approximating  the  data  as 

jci  =  0.2  -f  0.072  ,  (A.33) 

ds 

for  (Yst  x)/ds  <  6.  Here  Tc\  =  ( Tc\  —  Too)/(Tf  —  Too).  Then  pc\  was  found  from 


Poo  _  &Tf  fci  4-  Tqq 
Pci  Too 

to  give  an  approximation  for  the  centerline  velocity.  This  second  method  always 
gave  estimates  of  velocities  larger  than  the  isothermal  result,  Eq.  A. 30.  The  average 
of  the  two  methods  was  taken  as  the  estimate  of  Uc\(xj).  The  reason  for  not  simply 
using  the  more  complex  calculation,  which  accounts  for  buoyancy,  was  that  the 
equations  used  were  obtained  form  measurements  •>»;  attached  flames,  and  it  was 
felt  that  the  true  value  of  Uc\  would  lie  between  the  two  estimates. 


For  estimating  the  velocity  at  x  —  Z,  a  method  similar  the  last  one  was  used. 
First,  the  calculated  flame  length  Lc a  was  found  from  Becker  and  Yamazaki’s  Equa¬ 
tions  (40)  and  (58) 

(  0.202,  if  £  <  1 1 

=  \  >  ,  (A. 35) 

l  0.18  +  0.022£,  if  £  >  1  J 

where  rp  =  (d,  0/Lca.  MW,)  with  0  =  (MWm  T//MW,  T^)1/2.  Then  Eq.  A.32  was 
used  to  find  Uc\  at  Lca.  However,  Eq.  A.33  could  not  be  used  because  it  is  not 
applicable  above  (V^t  x)/ds  =  6.  Instead,  to  estimate  pc\  at  Lca,  their  Fig.  12  was 
used.  Three  quantities  were  obtained  from  the  graph:  the  maximum  normalized 
centerline  concentration  Td.m,  the  location  of  this  maximum  location  xm/Lc a  =  xm, 
and  x  for  which  Tcj  =  0.8,  x,8.  Then  linear  extrapolation  was  used  to  get  Tc\  at  Tc»- 


Because  the  measured  flame  length  L  was  in  general  less  than  the  calculated 
flame  length  £ca,  linear  interpolation  was  used  to  find  Uc\  at  Z,  using  the  velocities 
at  h  and  LCh  for  the  interpolation. 
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The  velocity  estimates  at  locations  other  them  h  have  rather  large  uncertainties, 
probably  ±30%,  and  are  intended  only  as  a  rough  measure  against  which  to  compare 
the  correlation  times  measured  in  the  experiments. 

Given  these  velocity  estimates,  local  passage  times  were  estimated  as 

0  AA  x 

n(x)  =  jff- T  ,  (A.36) 

Ud(x) 

and  the  convection  time  from  one  x\  location  to  another  x2  is  estimated  as 


x2  —  Xi 

1/2  (Ud(Xl)  +  Ucl(x2))  ' 


(A.37) 


A.8  Bound  on  Temperature  Rise  of  Reactants  from  Radiation 

A  very  crude  estimate  for  the  maximum  temperature  rise  in  the  reactants 
approaching  the  flame  base  caused  by  absorption  of  thermal  radiation  is  presented 
here.  The  result  is  quoted  in  Sec.  4.5  as  a  verification  that  this  temperature  rise  could 
not  be  a  coupling  mechanism  between  the  liftoff  height  and  downstream  radiation 
of  the  flame. 

The  following  nomenclature  is  used  in  this  section:  q  is  the  total  heat  release 
rate  of  the  flame;  fr  is  the  fraction  of  q  that  is  radiated;  AH  is  the  heat  release 
per  unit  mass  of  fuel,  assumed  to  be  50  x  106  J/kg;  a  is  the  absorptivity  of  the 
mixture  of  fuel  and  air,  integrated  over  the  spectrum  of  the  radiation;  and  rc  is  the 
convection  time  from  the  nozzle  exit  to  h,  integrated  along  the  ray  q  —  0.15. 


Approximating  the  radiation  as  a  point  source,  located  an  average  distance  of 
1.5  h  from  the  reactants,  the  temperature  rise  during  the  convection  from  the  nozzle 


exit  to  h  is 


AT  = 


frqa 

Poo  Cp  (1.5  h) 


2Tc  ’ 


(A.38) 


where  Cp  of  the  reactants  is  taken  to  be  1000  J/kg-K.  Noting  that 
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and  calculating  the  convection  time  rc  by 


fh  dx  _  0.08  h 2 

Jo  u(x,Ti)  fu(v)Usds 


the  temperature  rise  can  be  found  to  be 


AT  = 


0.08  AH 

1 6/.(i*)/r  Cp 


q  ds 


h 2 

(1.5  h)2  ' 


(A. 40) 


(A.41) 


Taking  /r  =  0.4.  and,  as  a  very  conservative  estimate  aa  =  0.01  m_1  results  in 
a  temperature  rise  of  only  0.002  K.  This  calculation  has  a  tremendous  number  of 
simplifications,  but  is  adequate  to  show  that  the  temperature  rise  associated  with 
absorption  of  radiation  is  quite  small. 
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APPENDIX  B 

Tabulated  Mean  Liftoff  Heights 


B.l  Measurements  Without  Air  Premixing 


Table  B.l:  h  Measurements  for  n.g.,  Without  Air  Premixing 


d 

(mm) 

us 

(m/s) 

h 

(mm) 

ua 

(m/s) 

h 

(mm) 

Us 

(m/s) 

h 

(mm) 

4.61 

33 

89 

44 

124 

56 

164 

39 

106 

50 

143 

6.25 

42 

133 

51 

159 

71 

213 

45 

141 

58 

175 

74 

229 

48 

152 

64 

195 

7.73 

39 

137 

59 

187 

81 

244 

49 

166 

70 

216 

90 

273 

Table  B.2:  h  Measurements  for  C2H6,  Without  Air  Premixing 

d 

(mm) 

ua 

(m/s) 

h 

(mm) 

(m/s) 

h 

(mm) 

Us 

(m/s) 

h 

(mm) 

3.12 

46 

74 

58 

98 

71 

124 

53 

88 

65 

110 

3.88 

48 

89 

66 

119 

87 

162 

53 

98 

73 

132 

93 

180 

59 

108 

81 

147 

4.61 

58 

115 

82 

154 

111 

210 

67 

130 

91  ! 

167 

74 

140 

101 

187 

117 


Table  B.3:  h  Measurements  for  C2H4,  Without  Air  Premixing 


d 

(mm) 


2.16 


Us 

(m/s) 

h 

(mm) 

Us 

(m/s) 

h 

(mm) 

73 

44 

92 

59 

83 

52 

102 

67 

95 

66 

134 

93 

116 

80 

149 

102 

111 

82 

145 

103 

125 

92 

168 

119 

Us 

(m/s) 


114 


126 


165 


181 


200 


h 

(mm) 


77 


91 


117 


134 


143 


Table  B.4:  h  Measurements  for  C2H2,  Without  Air  Premixing 


d 

(mm) 


1.02 


Us 

(m/s) 

h 

(mm) 

Us 

(m/s) 

h 

(mm) 

173 

20 

204 

28 

182 

23 

217 

31 

195 

26 

234 

36 

226 

32 

260 

38 

239 

34 

266 

40 

U% 

(m/s) 


240 


h 

(mm) 


37 
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B.2  Measurements  with  Air  Premixing 
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APPENDIX  C 

Correcting  Autocorrelation  for  Finite  Sample  Time 


One  of  the  possible  causes  of  error  in  the  current  measurements  of  h  vs.  time 
is  the  non-zero  sampling  time  of  the  photodiode  array.  Averaging  over  a  length  of 
time,  rather  than  talcing  an  instantaneous  measurement,  can  cause  changes  in  the 
measured  autocorrelation  function.  Thus,  both  the  rms  level  h'  and  the  correlation 
time  Tij 2  can  be  affected.  This  appendix  describes  the  method  for  estimating  and 
correcting  the  associated  errors. 

The  linear  photodiode  array  used  in  these  experiments  has  512  pixels  which 
are  scanned  successively  by  the  controlling  electronics  and  output  onto  a  serial  line 
which  is  then  input  to  the  A/D  board  of  the  computer.  After  scanning  the  last 
pixel,  there  is  a  delay  of  4  clock  cycles  (the  time  needed  to  read  4  pixels)  before  the 
next  scan  of  the  array,  starting  with  the  first  pixel.  Each  pixel  integrates  the  light 
incident  on  it  during  the  entire  time  between  successive  scans,  making  its  output 
equal  to  the  average  light  intensity  during  the  previous  array  sample  time  At  (ms 
per  scan).  Thus  each  pixel  intensity  is  an  average  over  a  slightly  different  time 
interval.  In  the  current  measurements  of  h,  however,  an  approximate  simplification 
can  be  made.  Because  h  moves  at  most  a  small  fraction  of  the  512  pixels  during 
any  time  step,  the  average  of  each  of  the  pixels  in  the  near  vicinity  of  h  is  taken 
over  approximately  the  same  time  interval. 

Two  steps  remain  in  finding  an  approximate  correction  to  the  measurements  of 
h.  In  Sec.  C.l,  it  is  argued  that  this  averaging  of  light  intensity  on  the  linear  array 
has  the  effect  of  averaging  the  liftoff  height  over  the  same  time  interval.  In  Sec.  C.2 
a  technique  for  correcting  the  measured  values  of  h'  and  Tj/2  is  presented. 
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C.l  Effect  of  Integration  of  Light  on  Measurement  of  h 

To  approximate  the  effect  of  the  integration  of  light  over  a  time  interval  on  the 
measurement  of  liftoff  height  h,  the  light  intensity  profile  is  linearized  about  the 
threshold  for  determining  h ,  It.  This  linearized  profile  is  given  by 

I(x,t)  =  It  4-  ^  (x-h(t))  ,  (C.l) 

where  ^  is  the  slope  of  the  I  vs.  x  curve,  assumed  to  be  constant  over  the  time 
interval  At.  Talcing  a  time  average  over  At, 

/(*)  =  It  +  ^  (x-h(t))  .  (C.2) 

Therefore,  at  x  =  h,  I  =  It,  indicating  that  averaging  the  light  intensity  /  over  a 
time  interval  At  has  the  effect  of  averaging  the  liftoff  height  h(t)  over  the  same  time 
interval,  subject  to  the  above  assumptions. 


C.2  Effect  of  Sampling  Time  on  the  Measured  Autocorrelation 

If  the  true  process  to  be  measured  is  designated  asi(f)  and  the  measurement 
occurs  over  an  interval  At,  then  the  measured  process  w(t)  can  be  written  as 

i  #>t+At/2 

w(t)  —  ~rr  /  x(t)  dt  .  (C.3) 

If  the  integration  in  Eq.  C.3  is  approximated  by  the  trapezoidal  rule,  w(t )  can  be 
approximated  by  z(t)  as 

w(t)  «  z(t)  =  l/2  [ar(#  —  At/2)  +  x{t  +  At/2)]  .  (C.4) 

The  autocorrelation  of  the  process  z(t)  can  now  be  written  in  terms  of  the  autocor¬ 
relation  of  x(t)  as 

Riz(t)  =  l/4[2 Rzx(t)  +  Rxz{t  +  At)  +  Rzx(t  -  At)]  .  (C.5) 
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Because  Rxx(At )  =  Rxx(—At),  for  r  =  0,  Eq.  C.5  reduces  to 

Rzz(0)  =  V2i?xx(0)+  'hRxx{At)  .  (C.6) 

Given  the  values  of  the  measured  autocorrelation  function  Rzz(iAt),  for  0  < 
i  <  N,  it  is  desired  to  find  the  corresponding  values  of  the  autocorrelation  of  the 
true  process,  Rxx(i  At).  In  order  to  do  this,  one  additional  equation  is  needed  at 
the  end  of  the  interval, 

Rxx{NAt)  =  Rz:(NAt)  .  (C.7) 

This  is  equivalent  to  extrapolating  the  values  of  RZZ(N  —  1)  and  RZZ{N)  to  RZZ(N  + 
1).  Equation  Eq.  C.7  completes  the  set  of  N  1  equations 

l/2Rxx(0)  +  l/2Rxx(l)  =  Rz:(  0) 

\/4RIX(i -1)  + 1/2  Rzx(i)  + 1/4  Rxx(i  +  l)  =  Rzz(i)  i  =  1,2, . . . ,  N  -  1 

R*r(N)  =  RZZ(N)  . 

(C.8) 

From  the  computed  autocorrelation  of  liftoff  height  an  attempt  was  made  to 
invert  the  above  set  of  equations  and  obtain  the  autocorrelation  in  the  limit  of  zero 
sampling  time.  However,  this  set  of  equations  is  extremely  sensitive  to  the  values  of 
Rzz(i),  in  particular  near  i  —  0.  Solutions  to  the  equations  typically  gave  oscillatory 
behavior;  values  of  Rzz(i)  for  even  values  i  were  higher  than  for  odd  values  of  i. 
An  alternative  method  of  correcting  the  autocorrelation  measurements  w as  then 
chosen. 

From  the  autocorrelation  of  the  data  set  with  the  best  relative  resolution 
(largest  r,^/ At),  Eqs.  C.8  were  used  to  compute  the  autocorrelation  which  would 
be  obtained  by  averaging  each  pair  of  adjacent  data  points.  From  the  new  auto¬ 
correlation,  different  values  of  h!  and  tx/ 2  were  found.  This  process  was  repeated 
until  tx/2/ At  was  less  than  that  of  any  case  studied.  From  the  autocorrelations  at 
the  different  effective  sampling  times  At,  a  least  squares  fit  was  preformed  on  the 
equation 

=  1  ,  _  (c.9) 

rl/2(0)  (Tl/2(At)/A<)C* 
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From  the  least  squares  fit,  the  values  of  C\  =  0.391  and  C2  =  1.042  were  obtained. 
This  equation  was  then  used  to  solve  for  Ti/2(0),  the  limit  of  rXj2  as  At  — ♦  0. 

Similarly  for  the  rms  levels, 


h’m  =  . ,  c3 

h'(  0)  (r1/2(Af)/Af)^ 


(C.10) 


with  C3  =  —0.113  and  C\  =  1.004  was  used  to  estimate  the  values  of  h'  in  the  limit 
At  -4  0. 


These  equations  resulted  in  at  most  a  5%  increase  in  the  value  of  h'  and  a  15% 
decrease  in  the  value  of  r1y2.  The  accuracy  of  this  correction  technique  was  checked 
in  two  ways.  First,  the  scheme  described  to  determine  the  values  of  C\, . . . ,  C4  was 
applied  to  another  autocorrelation  function  and  the  corrected  values  of  tx/2  and  h! 
compared  at  the  various  resolutions.  For  1.7  <  Ti/2(At)/At  <  6.0,  the  values  of 
h'  agreed  with  one  another  to  within  0.5%  and  the  values  of  rj/2  to  within  1%. 

For  the  second  test,  data  were  taken  at  two  different  sampling  times,  At,  but 
under  the  same  experimental  conditions.  The  correction  equations  (Eqs.  C.8  and 
C.9)  were  applied  to  h'  and  Tj/2  and  compared.  For  tx /2( At)/At  =  2.8,  the  (cor¬ 
rected)  value  of  Tj/2/t6  was  1%  higher,  and  the  value  of  h' /h  was  0.5%  higher  than 
at  Ti/2(At)/At  =  5.9. 

Having  confidence  that  the  correction  technique  was  accurate,  it  was  applied 
to  all  the  data  to  adjust  the  values  of  h '  and  Ti/2. 


124 


APPENDIX  D 

Accuracy  Estimates  and  Experimental  Considerations 


In  this  appendix  estimates  are  given  for  the  accuracy  of  various  measured  and 
calculated  quantities.  This  includes  the  measurement  of  flow  velocity  and  the  cal¬ 
culated  autocorrelation  time  rj/ 2,  as  well  as  the  large  scale  time  at  the  liftoff  height 
r$,  which  is  used  to  normalize  Ti/2- 

Along  with  identifying  the  size  of  errors  in  the  measurements,  some  analysis 
of  the  sensitivity  of  the  measurements  to  the  exact  experimental  conditions  and 
measurement  technique  is  presented.  This  includes  the  arrangement  of  screens 
around  the  ^cune  for  protecting  it  from  room  disturbances  and  the  choice  of  what 
fraction  of  the  flame  width  to  image  when  finding  the  liftoff  height. 


D.l  Measurement  of  Flow  Velocities 

Velocities  of  gases  exiting  the  nozzles  were  calculated  from  measurements  with 
laminar  flow  elements  (LFE’s),  across  which  the  pressure  drop  was  measured  with 
a  Barocel  pressure  transducer  as  described  in  Sec.  2.1.1.  The  exit  velocity  was  then 
calculated  by  assuming  that  the  gas  reached  atmospheric  pressure  and  temperature. 
The  sources  of  error  in  the  measurements  axe  1)  error  in  calibration  of  the  LFE,  2) 
error  in  measurement  of  pressure  drop  A p,  3)  error  in  knowledge  of  the  gas  viscosity, 
used  to  convert  A p  to  flow  rate,  and  4)  the  gas  not  reaching  a  uniform  temperature 
of  Too  at  the  nozzle  exit  (applicable  at  the  higher  Mach  numbers). 
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The  LFE’s  have  reported  accuracies  of  ±0.5%.  By  fitting  the  points  to  a 
parabola  of  A p  vs.  flow  rate,  this  accuracy  should  be  achievable  at  any  flow  rate  in 
the  usable  range.  In  the  current  experiments,  the  repeatability  was  approximately 
±0.5%,  and  the  accuracy  approximately  ±2%,  as  checked  by  comparing  the  two 
flow  meters.  The  pressure  transducer  is  accurate  to  ±0.1%  of  the  reading,  and  so 
it  should  not  contribute  substantially  to  the  measurement  error.  Viscosities  were 
taken  from  Reid  ei  al.  (1978).  In  comparing  them  to  the  CRC  Handbook  (Weast 
1974),  the  values  were  found  to  be  within  1%  for  all  gases  except  acetylene,  which 
differed  by  2%.  Because  the  flow-rate  is  approximately  proportional  to  Ap//i,  the 
corresponding  uncertainty  in  estimating  U0,  is  ~  2%  for  C2H2  and  ~  1%  for  the 
other  fuels. 

For  the  higher  Mach  number  flows,  the  largest  possible  error  in  calculating 
velocities  stems  from  the  assumption  that  the  exit  temperature  equals  Ta c.  Heat 
transfer  calculations  indicate  that  this  should  be  a  reasonable  assumption.  However, 
to  estimate  the  largest  possible  error,  calculations  were  also  done  assuming  adiabatic 
expansion  to  atmospheric  pressure  and  compared  with  the  isothermal  calculations. 
The  adiabatic  velocity  estimates  are  always  lower,  being  0.5,  2.0,  2.5,  and  6.5% 
lower  than  the  isothermal  velocities  for  the  highest  velocity  cases  of  the  four  fuels: 
natural  gas,  C2H6,  C2H4,  and  C2H2,  respectively.  The  calculation  of  ds  is  always 
larger  for  the  adiabatic  calculation,  by  a  percent  approximately  half  that  for  the 
velocity  calculation.  Therefore,  the  calculation  of  Us  ds,  on  which  both  the  Reynolds 
number  (Re0 Q)  and  large-scale  time  (t^)  calculations  rely,  would  be  only  0.25,  1.0, 
1.3,  and  3.3  %  lower  in  the  extreme  case  of  adiabatic  flow  to  the  nozzle  exit. 
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D.2  Resolution  Limits  Caused  by  Pixel  Spacing 

The  photodiode  array  used  in  the  measurement  of  h  has  512  pixels,  spaced  at 
50  /zm  intervals.  Each  pixel  has  a  width  of  2.5  mm.  For  each  flow  condition,  the 
placement  of  the  array  was  chosen  such  that  the  width  of  the  flame  image  just  filled 
the  width  of  the  array.  This  determined  the  distance  at  the  flame  corresponding  to 
the  spacing  between  two  pixel  centers,  Ah.  Specifically,  Ah/h  =  0.008. 

This  resolution  can  affect  the  accuracy  of  both  the  mean  liftoff  height  k  and 
the  other  quantities  derived  from  the  measurement  of  h  vs.  t.  The  time  averaging 
will  allow  h  to  be  measured  to  accuracies  better  than  Ah,  and  so  A h  is  an  upper 
bound  on  the  error  of  the  measurement  of  h. 


An  approximate  analysis  of  the  effect  of  the  pixel  spacing  on  measurement  of 
the  rms  fluctuations  h!  can  be  made  by  assuming  that  the  error  associated  with 
binning  the  true  liftoff  height  ht  to  an  integer  number  of  pixels  for  the  measured 
liftoff  height  hm,  e  =  hm  —  ht,  is  uncorrelated  with  ht.  The  rms  of  the  error  e  is 
then  e'  =  Ah/\/l2.  Using  the  experimental  result  that  h'  jh  as  0.07  results  in  the 
estimate  for  the  error  in  determining  the  rms  of  h  as 


h' 


e'  Ah  h 


Ah  h  h'  003 


(D.l) 


D.3  Errors  Caused  by  Integrating  Measurement  over  Time 

In  Appendix  C  the  technique  for  reducing  the  error  of  the  measurement  of  rj /2 
caused  by  the  non-zero  sample  time  of  the  linear  array  is  discussed.  Because  several 
approximations  were  made  in  deriving  the  correction  technique,  errors  will  remain. 

To  estimate  the  remaining  errors,  two  experiments  were  conducted  with  iden¬ 
tical  flow  conditions,  but  with  different  temporal  resolutions,  changed  by  adjusting 
the  clocking  rate  of  the  array.  As  acetylene  was  found  to  have  the  shortest  (dimen¬ 
sionally)  correlation  times  of  all  the  pure  fuels,  conditions  were  selected  to  match 
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the  resolution  of  the  acetylene  experiments.  For  this,  ethylene  was  used,  with  mea¬ 
surements  first  made  at  At  =  4.4  ms,  and  then  at  a  time  step  of  10.0  ms.  This 
longer  time  was  selected  to  match  the  resolution  of  the  C2H2  group  A  experiments 
where  Tx/2/At  %  2.3< 

The  autocorrelation  functions  were  computed  and  corrected  for  sampling  time 
as  described  in  Appendix  C.  There  is  very  little  difference  between  the  estimates  of 
T\ I2I Tb  f°r  the  two  cases  with  the  values  being  14.8  and  15.0  for  the  fast  and  slow 
sampling,  respectively. 

The  measured  rms  is  also  affected  by  the  sample  time,  and  was  corrected  by 
the  same  technique  as  the  autocorrelation.  Values  of  k' /h  are  0.0597  and  0.0600 
for  the  feist  and  slow  sampling,  respectively.  This  is  verification  that  the  reported 
values  of  h' /h  are  reasonably  accurate,  even  for  the  acetylene  flames. 

An  additional  test  of  the  temporal  resolution  effects  was  made  by  modifying 
the  data  acquisition  system,  allowing  just  256  of  the  512  pixels  of  the  linear  array  to 
be  sampled.  This  enabled  a  sampling  time  of  2.2  ms,  compared  with  the  standard 
sampling  time  of  4.4  ms.  The  results  from  a  mixture  of  86%  H2  /  14%  C2H6  again 
indicated  close  agreement  in  the  measured  autocorrelation.  Namely,  as  Tx/2/At  was 
increased  from  2.5  to  5,  the  measured  t1^2/t6  decreased  only  6%. 


D.4  Accuracy  of  Estimate  of 

T|/2  is  defined  as  the  time  at  which  the  autocorrelation  function  of  liftoff  height 
reaches  a  value  of  0.5.  To  estimate  the  accuracy  to  which  this  time  scale  is  known, 
the  following  procedure  was  used.  The  data  record  of  h  vs.  t  was  broken  up  into  eight 
pieces.  The  autocorrelation  was  then  calculated  for  each  piece,  and  the  resulting 
value  of  Tx/2  determined.  From  these  eight  measurements  of  Ti/2,  the  standard 
deviation  was  calculated.  To  estimate  the  accuracy  to  which  the  value  of  rj  j2  is 
known  when  using  the  entire  data  record,  it  was  assumed  that  the  eight  individual 
pieces  were  statistically  uncorrelated.  Under  this  assumption  the  standard  deviation 
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of  the  measure  of  t1  /2  for  the  entire  data  record  would  equal  the  standard  deviation 
for  the  individual  measurements  divided  by  y/8. 

As  an  example,  using  the  seventh  condition  listed  in  Table  4.2  (ethane,  d  = 
3.88  mm,  h  «  2.2)  the  standard  deviation  of  the  eight  individual  measurements  was 
found  to  be  4.5%  of  the  mean  value.  Dividing  this  by  >/8>  gives  an  accuracy  of  1.5% 
on  the  estimate  of  rXj2-  With  the  exception  of  the  eighth  case  in  Table  4.2,  during 
which  the  gas  supply  was  depleted,  the  accuracies  of  all  the  measurements  of  txj 2 
axe  estimated  to  be  better  than  ±3%. 


D.5  Calculation  of  Local  Large-Scale  Time 

The  large-scale  time  at  the  liftoff  height  7g(h)  is  defined  as  8/Uc\.  The  velocity 
was  estimated  as  described  in  Sec.  A. 3,  with  S  defined  as 

6  =  C6x,  (D.2) 

where  Cg  =  0.44.  For  most  of  the  cases  studied,  pc\(h)  %  Poo,  and  the  velocity 
estimate  reduces  to 

Uci  =  CUUS^-,  (D.3) 

ds 

with  Cu  =6.17. 

One  source  of  error  in  this  calculation  stems  from  the  fact  that  no  jet  virtual 
origin  x0  was  used  in  the  above  equations,  e.g.,  replacing  x  with  (x  —  x0).  If  both  XJC\ 
and  8  were  to  have  the  same  virtual  origin,  then  the  fractional  error  in  estimating 
rg  would  be  ~  2 x0/h.  Talcing  x0  —  2d  and  h  =  20  d  (the  lowest  liftoff  heights 
studied),  results  in  a  fractional  error  of  ~  20%.  For  the  highest  liftoff  heights  the 
error  would  be  reduced  to  ~  10%.  This  is  possibly  the  largest  source  of  error  in  the 
presented  data. 
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D.6  Fraction  of  Jet  Width  Imaged 

The  choice  of  what  fraction  of  the  flame  width  to  image  onto  the  photodiode 
array  could  potentially  affect  the  results  of  both  the  mean  liftoff  height  and  the 
temporal  behavior. 

Chapter  3  describes  measurements  of  h  vs.  t  that  were  made  for  relatively  short 
times,  where  only  the  mean  values  were  reported.  In  those  measurements,  for  the 
flames  near  blowout  (specifically  for  flames  with  Yrsth/ds  >  2.2),  the  outer  parts  of 
the  flames  were  chopped  off  of  the  image.  At  the  highest  liftoff  height  (Ysth/ds  ss 
2.6),  approximately  0.84  6  was  imaged.  To  see  what  effect  that  would  have  on  the 
mean  measurements,  two  experiments  were  conducted  which  were  identical,  except 
for  the  distance  of  the  array  from  the  flame.  This  changed  the  fraction  of  the  flame 
width’s  image  on  the  array,  as  well  as  the  resolution  of  the  measurement.  For  the 
first  case,  0.656  was  imaged,  and  for  the  second  case,  the  entire  jet  width  S  was 
imaged. 

The  fuel  was  natural  gas,  d  =  6.25  mm,  Us  =  61  m/s.  The  measured  value  of 
h  for  the  first  case  was  182.6  mm,  and  for  the  second  case  182.4  mm,  essentially 
identical  results.  This  indicates  that  the  measure  of  h  is  not  very  sensitive  to  the 
fraction  of  the  flame  width  imaged,  at  least  down  to  w/6  =  0.65,  where  w  is  the 
width  imaged. 

The  same  measurements  are  also  able  to  give  an  indication  of  how  sensitive  the 
temporal  measurements  are  to  the  fraction  of  the  flame  width  imaged.  The  auto¬ 
correlations  from  the  two  cases  are  plotted  in  Fig.D.l,  along  with  a  third  curve  to 
be  discussed  in  Sec.  D.7.  There  is  no  significant  change  in  the  autocorrelation  when 
the  fraction  of  the  flame  imaged  is  reduced  to  0.65.  Whether  a  further  reduction, 
say  to  just  a  thin  strip  along  the  centerline,  would  change  the  time  scales  measured 
has  not  been  investigated. 
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FlG.  D.l  Autocorrelation  of  liftoff  height.  Natural  gas,  d  =  6.25  mm,  Us  =  61  m/s. 

Solid  line:  fraction  of  width  imaged  w/6  —  0.65,  dashed  line:  w/6  =  1, 
dot-dash:  same  as  dashed,  except  a  larger  screen  diameter  was  used. 

D.7  Effect  of  Diameter  of  Screen  Surrounding  Flame 

To  see  how  the  time  scales  of  the  h  fluctuations  are  affected  by  the  surrounding 
cylindrical  screen,  measurements  were  made  in  which  a  screen  with  a  larger  diameter 
was  used,  to  compare  with  results  of  the  standard,  smaller  screen. 

These  results  from  the  larger  82  cm  screen  are  compared  to  those  of  the  68  cm 
screen  in  Fig.  D.l.  As  there  is  no  measurable  difference  in  correlation  times,  the 
reported  results  are  relatively  insensitive  to  screen  diameter  used,  as  long  as  the 
screen  is  large  enough.  The  possibility  that  the  large  correlation  times  measured 
scales  in  some  way  with  the  screen  diameter  can  therefore  be  ruled  out. 
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APPENDIX  E 

Photographs  of  Lifted  Flames 


Photographs  of  natural  gas,  ethylene  and  acetylene  flames  are  shown  in  this 
appendix.  All  photographs  were  taken  with  a  Nikon  35  mm  camera  using  T-max 
3200  speed  film  pushed  to  6400  during  development.  The  range  of  liftoff  heights  in 
the  photographs  of  each  fuel  corresponds  roughly  to  the  range  1.8  <  h  <  2.6. 

Figure  E.l  shows  natural  gas  flames  using  a  nozzle  diameter  of  4.61  mm,  taken 
with  an  exposure  time  of  1/1000  s,  at  f/1.4.  Flames  at  four  different  flow  velocities 
are  shown,  with  two  photographs  at  each  velocity,  e.g.,  (a)  and  (a'),  except  at  the 
highest  velocity  for  which  only  one  photograph  was  taken.  From  these  photographs 
it  is  clear  that  the  flame  base  is  not  a  flat  front.  Rather  it  can  be  quite  irregular. 
This  is  particularly  evident  in  parts  (c)  and  (d)  of  Fig.  E.l,  which  are  flames  rather 
close  to  blowout.  In  some  of  the  photographs,  tongues  of  flame  protruding  down 
toward  the  nozzle  can  be  seen.  As  an  aid  to  judging  the  degree  to  which  these 
photographs  represent  instantaneous  images,  the  ratio  of  the  exposure  time  to  Tf, 
ranges  from  ~  0.13  to  0.10  in  going  from  (a)  to  (d). 

Photographs  of  C2H4  flames  are  shown  in  Fig.  E.2,  taken  with  an  exposure 
time  of  1/4000  s,  at  f/2.0.  The  bright  regions  in  the  tops  of  the  photos  of  (a)  and 
(a')  are  radiation  from  soot,  and  the  visible  light  is  yellow,  wheres  the  flame  base 
appears  blue.  The  soot  level  decreases  as  the  blowout  velocity  is  approached.  The 
flame  base  is  again  irregular,  although  it  has  a  slightly  different  appearance  from 
the  natural  gas  flames.  The  undulations  are  smaller,  giving  the  base  a  rough  shape, 
but  with  less  large-scale  deformation.  For  these  images,  the  ratio  of  the  exposure 
time  to  Tf s  ranges  from  ~  0.19  to  0.12  in  going  from  (a)  to  (d),  indicating  only 
slightly  poorer  resolution  than  the  natural  gas  photos  of  Fig.  E.l. 
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Finally,  photographs  of  acetylene  flames  are  shown  in  Fig.  E.3,  taken  with  an 

exposure  time  of  1/4000  s,  at  f/8.0.  The  bright  regions  downstream  of  the  initial 

flame  base  are  again  caused  by  soot.  These  soot  regions  axe  reduced  as  the  blowout 

► 

condition  is  approached,  and  in  (d)  and  (d')  the  soot  is  not  visible  at  all.  The  flame 
base  is  very  bright  (blue)  and  the  region  just  downstream  is  dim  (not  even  visible 
in  the  photo  of  (d) ).  The  exposure  time  divided  by  r*  ranges  from  ~  0.91  to  0.62 
for  these  images,  indicating  that  there  is  some  blurring. 
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